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ABSTRACT 


Satellite operations scheduling, like many other real life-scheduling undertakings, is a 
complex problem. Although NP-complete, it is more challenging than, say, job shop or open 
shop, due to the dynamic nature of the constraints involved. It is difficult to derive here the 
objective function(s). Many qualitative decisions are also involved. Therefore, this problem 
cannot be solved by conventional optimisation methods and the availability of many textbook 
models appears to be of little direct utility in satellite operations scheduling. 

In this sUidy, we consider the problem of optimally resolving satellite visibility cla.shes over a 
ground station. Such opportunities develop about 250 to 270 times each week in ISRO's LEO 
satellites network consisting of 8-10 ground stations, all individually tracking and 
communicating with the satellites. We propose a non-linear objective function for scheduling 
support of TfC type of operations. It also covers the issue of dynamic priority. Then three 
meta-heuristic approaches, viz. Greedy (local) Search, Tabu Search and Simulated Annealing 
are used to tackle the clash resolution problem. Each of these methods is able to find solutions 
for clashing satellite visibilities that are optimal or near to optimal. To use local search and yet 
avoiding getting trapped in a local optima, wc used a new method of generating diversely 
dispersed initial solutions, one in each disjoint search region. This enabled reaching near-to- 
global optima in all runs. The study concluded with the comparison of the three above 
approaches and the GA approach used earlier to solve the clash resolution problem. Sample 
program codes in JAVA are attached. 
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Chapter 1 

The Challenge of Scheduling Satellite Operations 

1,1 Introduction to Satellites and their Motion 

A satellite is an object which revolves around another object. For example, the moon 
is a satellite of the earth and the earth is a satellite of the sun. Man-made or artificial 
satellites are those which orbit round the Earth in their fixed orbits. The shape of this 
orbit may be elliptical or close to circular. 

India is currently having two series of satellites 

• Indian National Satellite (INSAT) Series 

• Indian Remote Sensing (IRS) Series 

Satellites of the INSAT series including the latest launch of GSATl by GSLVl are 
gco-slatioiuuy, moving at a height ol' around 36000 kin, and therefore having a time- 
period of 24 hrs. I'hese satellites are used for communication, broadcasting and 
meteorological purposes. 

Satellites of IRS series are low earth orbit (LEO) satellites moving at a height of 
700 to 900 km. 'I'heir time period is around 100 minutes. IRS satellites are used for a 
large number of remote sensing applications in areas such as agriculture and crops 
management, irrigation planning, Ibrest and bio-resources hianagement, ocean 
study, rural development, urban planning, etc. I'he present constellation of satellites 
of IRS series consists of lA, IB, 1C, ID, P2, P3, P4, and SROSS. 


Most Indian satellites are fabricated at ISRO Satellite Centre (ISAC), Bangalore while 
the launch vehicles are fabricated at I'iruvananthpuram. As stated above. Low Earth 
Orbit (LEO) satellite orbits at 700-900 kilometres above the earth, resulting in very 
low latencies. While this seems to be a simple solution, it is not that simple to 
implement LEO technology. As one brings the satellite closer to earth, the number of 
satellites needed to cover all areas on the earth increase and high-speed 
communication handoffs between satellites become complex. A LEO satellite may be 
visible for half an hour before it needs to hand-off the signal to the next incoming 



satcllilc. I'liis means that the antenna at the customer's premises has to be sufficiently 
inlcliii'.c-nl to ilelcel and swileh fioin the iiassiiip, salelliles. 

1.1.1 Scheduling the Launch [13] 

'fhe launching of Indian satellites into required orbit is done from Shriharikota. 
l.aunching hundreds of multimillion-dollar satellites is a challenge. Teledesic had 

initial plans of a constellation of 840 satellites but that has 
now dropped to 288. It plans to launch all its satellites within 
a two-year time frame; Celestri has simileur plans. A 
staggering 1,700 satellites will be launched in the next 
decade— more than 10 times the 150 commercial satellites 
now in orbit. This creates big capacity problems for 
launching these satellites. The existing services just do not 
have the capacity. In 1998 alone, we saw more than 60 launch carrying more than 100 
satellites into orbit. To exacerbate things further, the launch reliabilities haven't 
improved since the dawn of space age and 1 in 10 launches fail. This makes the 
insurance premium on a satellite skyrocket to the point where it begins to rival the 
cost of the satellite. 

The main purpose of IRS series of satellites is nation’s resource monitoring and its 
management. Also, satellites generate revenues by taking images of various 
locations on the globe as requested by global customers. 

Although the path of a satellite is a fixed orbit, by virtue of spiq motion of earth, its 
ground trace is spiral shaped. Also, the spiral traversed each day is not exactly the 
same as that of previous day. The spiral path covered on a particular day is repeated 
after a fixed period of time, which may be different for different satellites. This 
period of lime is called ‘repent cycle’ of the spacecraft. (IRS-IM -> 2 day.s, IRS-1 C 
-> 24 days, etc.) In other words, repeat cycle of a satellite is a set of different paths 
the satellite can follow. 

As explained above, a satellite covers and serves different points of globe at 
different points of its motion. Therefore, it is very important to monitor the health 
parameters of various sub-systems of the satellite continuously and also keep track 
of its orbit. This is achieved by means of special “chain” or ground station 
configurations set up at geographically distributed locations that can establish radio- 




link with the spacecraft and thus can communicate (send and receive signals) with it. 
1 hcsc configurations are commonly called ground-statiom. 
t iround-slalions are exclusively of two types- 


• Telemetry, Tracking and Commanding (TTC) stations, and 

• Payload (PL) stations. 

A set of two or more stations is called ground network. 

1.1.2 Advancements in Micro Satellites: [15] 

Micro satellites provide a number of unique advantages not found in larger satellites, 
but also suffer from a number of constraints. The most significant advantage is the 
cost of launching which is around $6 000 per kg compared to $40 000 per kg for 
larger satellites. Given the cost advantage and the miniaturisation of electronics and 
mechanics to build in the required functionality, the constraints amount to power 
availability and opportunity for launching. 

The major limitations on micro satellites arises because of following constraints 

1 . Orbits and launchers 

2. Power 

3. Construction 

4. Architectures 

1.2 Operations on Satellites 

Satellite operations may be classified as follows: 

• 7TC operations: These operations are required for the health monitoring of the 
satellite subsystems. 

• Payload operations: Payload operations are used in imaging the ground from 
LEO satellites. 

• Command Operations: These are of very small duration and are operations 
performed to send signals to the space craft for resetting its various devices, 
.swilching ON/Ol'F the cameras, commanding the timers, etc. 



All TTC stations are equipped with almost identical systems of telemetry, reception, 
tracking and commanding. The I fC network also provides dedicated services for 
payload operations and spacecraft health keeping throughout the mission life. Payload 
stations receive data from spacecrafts, as scheduled. 

1.3 What is a “clash”? 

In MSS software documentation [1 7] the definition of a visibility clash given by 
ISTRAC is “Over a ground station if two or more spacecrafts pass simultaneously or 
one after the other such that the time difference between first spacecraft LOS and 
second spacecraft .A.OS is less than the station reconfiguration time, then the two 
spacecrafts are said to have a visibility clash”. 
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Fig 1.1 Visibility Clash 





Also when a single satellite is visible from two or more ground stations, its important 
to decide which station will support that particular satellite. This situation is called a 

station clash. 



The other clash situation arises when two or more satellites controlled by the same 
controller require support at the same time. The controller can’t give commands to 
two satellites at the same time and hence the schedule must decide which satellite is to 
be supported. The situation is called a controller clash. 



1 .4 Constraints that affect task scheduling. 


The problem we wish to tackle in this investigation is the optimum allotment of tasks 
to open ‘windows’. Here an open window develops due to the interaction of a satellite 
and a ground station as the satellite ‘passes over’ that station. This is called ‘visibility 
window’ as the satellite is visible to the ground station during this time interval. The 
scheduling of operations in an open window is a difficult task. The complexity arises 
due to the nature and types of constraints and interrelationships among different 
operations. 

It is difficult to find a direct match of this problem to conventional machine 
scheduling problems described in the literature as flow, job or open shop. Various 
types of special constraints exist in satellite support scheduling. 

1 .4. 1 Constraints induced due to the nature of satellite operations. 

• Some operations should be scheduled in requested visibility. These include the 
requested payload operations. 

• Some operations require other operations as prerequisites and/or post- 
requisites. 

• A special category of operations is ‘routine’ operations. These operations are 
performed on routine basis. However, these operations have preferred and 
feasible passes defined for them. 

• The next category is the ‘special request’ operations. These operations come 
very rarely but when present, these must be done in the requested pass. This 
may force some other operations to be postponed or preponed. 

1 .4.2 Constraints induced due to ground stations: 

• An operation can only be performed over a ground station if the station has the 

required capability to perform that specific operation. 

• A station can support a satellite only if it has capability to support that specific 

satellite. 

• Each station has its own reconfiguration time that is a big constraint. 



1 .4.3 Satellite induced constraints: 

• Satellites have their own duties and priorities defined. For example, IRS-1 C is 
a payload satellite while SROSS is meant for scientific operations. 

1 .4.4 Controller caused constraints: 

• In performing a control operation, all control signals are given by a human 
controller. A single controller may control more than one satellite. However, 
the controller can’t give commands simultaneously to two different satellites, 
and hence the previously described controller clash situation arises. 

1 .4.5 Constraints induced due to visibility scenarios: 

• A scenario of visibilities can be a complex set with many different satellites 
visible at the same time at a given ground station. Other possibility is that 
same satellite is simultaneously visible from two different ground stations. 
These are complex clash situations that have to be resolved. 

1.5 Thesis Objective: 

1. To develop an objective function for scheduling support of TTC type of 
operations that will tackle the issue of dynamic priority by taking into 
consideration various factors discussed with ISRO. 

2. To check the performance of different metaheuristic algorithms such as 
Tabu Search, Simulated Annealing and Greedy Search on the problem of 
clash resolution of satellite visibilities over a station in the satellite 
scheduling system. 

3. To implement the algorithms on sample problem and compare the results on 
the basis of following parameters. 

• Final solution obtained. 

• Time taken 

• Consistency of algorithm. 

4. Compare these results obtained with the results already obtained by applying 
GA to the same situation. 



5. I'o produce a weekly schedule of clash resolved support using the 
best optimisation method. 

1.2Thesis Organization; 

Chapter I provides an introduction to the concept of satellites and its motion. This is 
followed by the problem overview and the literature review in Chapter 2. Chapter 3 
deals with the development of the objective function to help resolve visibility clashes 
for TTC type of operations and the mathematical formulation of the problem given by 
Kumar and Bagchi |3j. Chapter 4 gives an alternate approach proposed in this study to 
solve the problem using the Greedy search heuristic, its implementation, and the 
results. Chapter 5 discusses a method to solve the problem using Tabu Search, its 
implementation, and the results. In Chapter 6 we have presented the Simulated 
Annealing (SA) approach for the same problem. Also, we have optimised the critical 
parameters in SA. In the same chapter, we discuss the implementation and the results 
obtained by using SA. In Chapter 7 we discuss the approach to solve the same 
problem using GA as given by Kumar and Bagchi [3]. In Chapter 8 we compare the 
first three methods. Then we compare the results and computation time of these three 
algorithms with GA. The JAVA programs for these algorithms are given in 
Appendices. 



Chapter 2 

Problem Overview and Literature Review 


2.1 Problem overview: 

A satellite-scheduling problem is the problem of assigning start times, end times and 
resources to a set of constrained tasks, such that multiple mission objectives are 
satisfied, i'he mission tasks and resources are subject to numerous physical, geometric 
and operational constraints. Satellite mission planning problem differs from 
traditional scheduling problem in several ways [6]. Perhaps the most obvious 
difference is the fact that some of the resources (i.e. satellites) are orbiting the earth. 
This places an additional set of constraints on when a task can be executed. Also, 
satellite mission planning problems may involve task start time' preferences, 
consumable resources, periodic tasks, variable duration tasks, and tasks that can be 
pre-empted. In addition these problems are usually over constrained. 

Kumar and Bagchi [3] have described the problem addressed here in terms of 
scheduling terminology as, 

• The opportunity window is open for a specified period within which the 
resources must be allocated so as to maximise the total profit contributed. 

• Products are perishable, have no stock/salvage value, and consumed as 
produced. 

• Each product has its own marketing opportunity. 

• Market windows of products may clash, forcing production of only one at 
one time. 

• If a product is produced, a specified minimum quantity of it must be 
produced. 

• Once started, production can’t be interrupted. 

• For each product there is a maximum quantity that must be produced. 

• A finite reconfiguration time is required to switch the machine from one 
task to another. 

• The machine itself is available for a specified period. 

• An arbitrary profit function determines how much profit each product 
contributes. 



A simple version of this problem resembles parallel machine scheduling. The goal is 
to find an optimal set of start times, such that the capacities, availability and time limit 
constraints are met and as objective function is optimised. 

A visibility clash occurs when one or more satellites visible at a ground station have 
cither overlapping visibilities or their respective AOS (Acquisition of signal) and LOS 
(Loss of signal) difference is less than the station reconfiguration time. The problem is 
to find the optimum support schedule that maximises the returns (or the total value 
generated) from supports. A clashing situation is difficult as the patterns in which 
satellites moves are enormous. The solution search space is large and discontinuous. 

1 he problem is NP-complcte. [3J 
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2.2 Literature Review 


In India, there is not much work going on over the satellite-scheduling problem. Some 
researchers such as Kumar and Bagchi [3] have worked on this problem. In ISRO, 
some work has been done for satellite scheduling and clash resolution. They have 
developed software called MSS (multisalellite scheduling software) which generates 
the operations schedule and allocates 'FTC network resources for different satellites 
supported at Spacecraft Control Centre (SCC). Satellite specific constraints, ground 
station configurations, satellite priorities, mission requirements and priorities of 
certain payload and special operations, as well as visibility conflicts are taken into 
consideralion wliilc generating operations schedule. 1 17] 

The clash resolution logic given by JSTRAC is as follows [17J 

The CLASH module considers the sorted visibility records of all the spacecrafts for 
one day and these records are divided orbit-wise into different sets. The records in a 
set arc also arranged station-wise and for each station the records that do not have any 
visibility clash are determined. Such records are called confirmed records i.e., the 
support for the spacecraft is confirmed. The weights of all the confirmed records are 
computed. These weights have a significant role to play at the time of clash 
resolution. The Confirmed Record Weightage (CREW) is a function of maximum 
elevation, no. of operations, station priority and operation priority. 

CREW = f (Maximum elevation , no. of operations, station priority, operation 
priority) 

for each ground station among the available clash records, clash weights are derived 
based on the spacecraft priority. The clash weightage (CW) is a function of maximum 
elevation, spacecraft priority, no. of operations, operations importance, exclusive pass, 
minimum operations met, next pass visibility, previous pass visibility, last pass of the 
day, CREW . 

CW = g (Maximum elevation. Spacecraft priority. No. of operations, operations 
importance, exclusive pass, minimum operations met. next pass visibility, previous 
pass visibility, last pass of the day, CREW). 



The formulation for deriving CREW and CW has been mathematically modeled b)< 
taking in to account the function variables and spacecraft control parameters. Having 
derived clash weightages for all the clash' records, the combination is selected for 
which the clash weightage is maximum, i.c. maximize the weightage function (WF ) 
defined below: 


m ti 


WF= SE 

'^1 7-1 


CWij 
2* j 


m ~ no. of spacecrafts. 

n= no. of stations for spacecraft 


Pickup the i and j for which the above weightage function is maximum 
This procedure identifies one optimum combination of records for which the visibility 
clash is resolved, the remaining records in the set are tagged as ' NO SUPPORT'. The 
same process is carried out for the entire data considered for scheduling. 

Kumar and Bagchi [3] have worked on this problem and have evolved with a GA 
((ienetic Algorithm) based solution to tlie visibility clash resolution problem. In the 
paper presented {OR-2()0() XXXIIl Annual Convention of the OR Society of India, 
Ahmedahad, Operations management conference 4, Madras), they have described a 
new method for optimising ground station support to LEO (Low Earth Orbiting) 
satellites. The methodology invokes Darwin’s principles of natural evolution to 
handle objective functions that are arbitrarily shaped and constraints that do not 
possess tractable structure. They have indicated up to 60% improvement over 
manually developed support plans. I'heir method can also help find optimal location 
for ground stations and support capability deployment in widely diverse scenarios. 

Wolfe and Sorensen [19J have described three different approaches to assigning tasks 
to earth observing satellites. A fast and simple priority dispatch method is described 
and shown to produce acceptable schedules most of the times. A look ahead algorithm 
is then introduced that outperforms the dispatcher by about 12% with only a small 
increase in run time, fhe third method they used is CiA. The genetic approach uses 
two additional binary variables, one to allow the dispatcher to occasionally skip the 
job in the queue and another to allow the dispatcher to occasionally allocate the worst 



posilion to the job. I’hese variables are included in the recombination step in natural 
way. They implemented the algorithms on “window con.strained packing problem” 
which is similar to the station clash resolution problem. 

The Cicncric Resource, livent and Activity Scheduler (GREAS) Application 
I 'ramcwork (A1- ) is a robust collection of software components for building satellite 
mission planning applications. [9J The GREAS Scheduling Framework utilizes a 
constraint programming approach for solving these satellite mission-planning 
|)roblcms. Gonstraint programming uses information contained in the problem to 
“prune” the search space, rapidly identifying feasible solutions. It is ideally suited for 
opcratit)nal problems, which require fast, feasible answers. Constraint programming 
methodology allows the natural expression of very complex relationships, including 
logical expressions. The Scheduling Framework builds upon the constraint 
programming capabilities of commercial software ILOG Solver and ILOG Scheduler 
to meet the challenges of solving satellite mission planning problems. 

The Space Telescope Science Institute [10] has developed a general constraint based 
scheduling framework for scheduling operations of NASA’s Hubble Space felescopc. 
The approach is to treat the problem as a constrained optimisation problem. It uses a 
hcuiTstic repair-based scheduling technique called multistart stochastic repair. 1'his 
technique first makes a trial assignment (or initial guess) of activities to times, based 
on heuristics. Such a schedule will generally have temporal or other constraint 
violations, as well as resource capacity overloads. Then they apply heuristic repair 
techniques to try to eliminate constraint violations, until either a pre-established level 
of effort has been expended or there are no conflicts left. The third step is to eliminate 
conflicts by removing any activities with constraint violations, or by relaxing 
constraints, until a feasible schedule remains. The heuristics employed in are 
stochastic. 

Agnese and Brousse [IJ apply depth first, branch and bound methods and greedy 
search to set task allocation to satellite visibility windows. The exact methods such as 
depth first, branch and bound etc. have capability of providing optimal solutions at the 
cost of computational time. 



Bottcher, Jahn, Lutz and Werner [4] have worked on basic design problems of 
LLO/lC'O-based networks. In first part of the paper titled “Analysis of basic system 
parameters of communication networks based on low earth orbit satellites”, the 
estimates for the necessary number of satellites, orbits and number of communication 
channels per satellite are derived. Since the latter is a crucial quantity, the number of 
communication links and channels per link are derived with a more elaborate model in 
the second part of the paper. This includes the radio links from the satellites to mobile 
users and to gateways, as well as intersatellite links and terrestrial lines. They also 
introduce a formal model for LEO/ICO- based networks and propose a method for the 
evaluation of link capacities, given the network topology and the traffic requirements. 
As an example, two constellations are investigated in detail. One of these 
constellations is the IRIDIUM system proposed by Motorola. 


Kennedy, Gregory and William [11] have developed a prototype expert system for 
scheduling satellite supports for satellites in Global Positioning System. I'he 
prototype used Texas Instruments’ Personal Consultant Plus software and 
incorporated a number of unique search routines and graphics. This is described in 
detail below. 

Scheduling communications contact between satellites and tracking stations ( satellite 
support) is a complex manual process that requires specialized knowledge and 
expertise. A satellite support occurs when a line-of-sight communication link is 
established between the satellite and a ground remote tracking station (RTS) to 
transfer electronics data down. Each RTS can communicate with one satellite at a 
time, rhe length of time available for a satellite support (support window) varies from 
a few minutes for a low-earth orbit satellite (200 KM) to continuously for a 
geosynchronous satellite (36,000 KM) visible at an RTS. In geosynchronous orbit, the 
satellite’s location remains fixed with respect to the surface of the earth. 

The planner analyst (PAs) determines the number of support required for each 
satellite and the preferred schedule. The range schedulers must schedule the RTSs to 
support 50 to 100 satellites. The scheduling objectives are to minimise the total 
number of supports required to accomplish each satellites mission and to provide R'fS 
schedulers as much as possible. 



The Air Force Satellite Control Facility at Onizuka Air Force Station, California, 
rcciucsled to the authors to examine the feasibility of using expert systems (IvSs) to 
assist the planner analysts. They used ES technology to develop a prototype advisor 
.system to help planner analysis schedule .satellite supports. They selected the CHobal 
Positioning System (GPS) satellite program for prototype development for two major 
reasons: first, each of the 20 satellites in the GPS system requires separate satellite 
support scheduling, and second, the GPS program if of medium scheduling difficulty 
because the periods GPS satellites are “visible” to an RTS are shorter than for 
geosynchronous satellites and longer than for low-earth-orbit satellites. 

Their objective was to demonstrate that, using ES technology for scheduling satellite 
supports was feasible and cost effective. Based on the success of this prototype, 
several ES projects have been initiated by the Air Force and supporting contractors. 

fhe satellite support scheduling Process 

Scheduling satellite supports is a three way process: First, satellite users send requests 
to the Mission Control Complexes (MCC) for a satellite to perform particular task at a 
certain time. I'hcre arc several MCCs, each responsible lor the satellites performing a 
parliciilar mission. Each .salellile is assigned a unique operations number, for example 
5111. 

Second, PAs prepare seven-day schedule of the supports required for each satellite 
(referred to as Program Action Plans or PAPs). For each satellite, PAs prepare cover 
sheets that outline the support requirements for the seven-day period. Next, PAs 
prepare satellite visibility charts (llg), which indicate all times each .satellite is visible 
at each RTS for the first 24- hour period. Finally, they generate PAPs which include 
the following information: satellite vehicle operations number, planning period, date 
of each support, start/stop times of the support window, contact duration time, visible 
R'l'Ss throughout the support window, type of support, and special support 
requirements (for example, unique equipment). 

Third, each satellite’s PAP is transmitted to the Range Control Centre, an office where 
schedulers build seven-day and 14 day schedules assigning each satellite support to a 
vi.sihle RTS. 
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Figure 1: A sample planner analyst's visibility chart which indicates the times that satellite 
#5111 is visible at each of the seven RTSs for a 24-hour PAP period. Each satellite is assigned a 
unique tape color pattern and rolls of tape with that color pattern are used by the PA 
throughout the support scheduling process to show the duration of the support windows. The 
exact start and end times of visibility are handwritten at each end of the tape strips. 


Why Use Expert System Technology? 

Expert systems are appropriate for three reasons: (I) expert knowledge is important in 
developing PAI’s, (2) the scheduling process is very rule oriented, and (3) automating 
the PAPing process should save money. 

Only certified PAs prepare the PAPs. Each PA is trained on the job for three to six 
months and must pa.ss qualifying examinations to become certifted. Because of the 
complexity of the process and the importance of the satellite supports, each PAP is 
checked by a second certified PA for aceuracy and completeness. 

In preparing PAPs, PAs use several documents and also rely on their extensive 
experience and expertise. They refer to computer printouts of ephemeris data, which 
pni\iik' IMS x'isihililics niul orbit f-’.conuMrics for each satellite vehicle. Various 
documents (test operations instructions and memograms) provide the rules for each 
type of support. These heuristic methods and the scheduling rules can be modelled 
using ES technology. 

Finally, managers at the Air Force Satellite Control Facility believe that there is a 
clear economic incentive to automate the process. Currently the process takes 20 
hours per week per satellite, and 20 satellites are planned for the GPS system. An ES 
could reduce the 400 hours per week by an order of magnitude. 
















re 2: The prototype system design uses a hierarchical frame knowledge representation 
s .cture. Each frame represents a sublask to be performed and contains the procedural and 
heuristic rules to perform each task. When a frame is evaluated it also has access to knowledge 
contained in parent frames. The knowledge representation graph is searched in a depth first 
manner [Nilsson 1980). 


K n () w! edge . A cci u i si lion 

Acquiring knowledge about the existing planning process was critical to their work 
developing a prototype system. They asked the Air Force Satellite Control Facility to 
assign to the project a PA expert who was very interested in it. After extensive 
planning, the knowledge engineers made two one-week site visits. Afterwards, the 
expert reviewed transcripts of the knowledge they had acquired during our visits. 
They continued to consult with the expert by phone each week throughout the 
research. 

Prototype System Design 

Based on an extensive trade-off analysis of ES shells, they selected LISP-based 
Personal Consultants Plus (PC Plus). PC Plus [Texas Instruments 1986] provides for 
forwiird and backward chaining, frames and rules, explanations, user-defined 
lunctioiis. and the graphics interfaces necessary to develop the PAP I'orms. 












( iMiiurc 2) sliDws ihe knowledge reprcsentalioii scheme. For each of the boxes 
ill lii’uie, (licrc is li.imc lliat has unit|iie rules, [latciil frames, goals, user 
prompts, and parameters, riie ES begins the consultations in the schedule frame and 
leiiuests the user to identify the satellite vehicle and the PAP period date. The ES 
then computes the remaining dates. 

'I’he visibilities frame is evaluated next. The ES transfers the appropriate visibility 
data from an external data disk into the knowledge base. Next the display-vis frame is 
evaluated. The ES asks the user whether he or she wants to display the satellite 
visibility charts for any day during the PAP week. 

The support-requirements frame is then evaluated to determine the required support 
(for example, satellite state of health (SOH)). Next, each applicable support frame is 
evaluated. Each of the 18 support frames contains the unique parameters and rules 
necessary to determine such PAP information as support start time and duration. The 
prototype implements SOH supports only but is designed to accommodate all 18 
support types. The other 17 support type obtains data on specific satellite subsystems, 
for example, the battery, the thermal control system, and the rate gyro system. 

The ES then evaluates the master-schedules frame, which contains the scheduling 
parameters and rules. Visibility search routines identify the RTSs that are visible at a 
particular time. If the satellite is not visible from any RTS at the stated time, these 
routines search for the next best time that the satellite is visible from one or more 
RTSs. Each frame evaluation produces the PAP schedule for one support. 

When all supporters are scheduled,, the display-PAP frame displays the completed 
PAP form. At the expert’s request, all graphic displays were designed to mirror the 
current manually generated PAP documents. 

User Intcilacc 

The user interface is menu-driven. To generate a PAP, the user must answer 10 to 15 
questions by the expert system. A sample question is 
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Figure 3: A sample visibility chart developed by the ES shows the same data that was manually 
generated in Figure 1. The chart shows the visibility times on June 16, 1986 for #5111 at each of 
the seven RTSs. 
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I if'.iiro 4: A s.MnpIc riofp.uii Action I’l.in form clcvclopcd by the expert xystern. The expert 
system scheduled the same remote tracking stations as the planner analyst. 








Enter the operations number of tlie vehicle to be PAPed: 

5111 

5112 
..9763 

Ecstilts 

Figure 3 is a sample ES visibility chart for satellite #5111 on June 16, 19S6. This 
chart shows the same visibility data as the manually developed chart in, figure 1. The 
I'iS saves the PA 30 niinutcs per siitellite per week. 

A sample of a completed PAP is sliown in figure 4. Given the last scheduled support 
lime from the user (June 15 at 2200 hours), the ES determined that a support was 
required once every six hours and that the best time window for the next support 
occurs at 032()-()400 (I'igure 4). Tlien the l-iS invoked the visibility search routines to 
determine if vehicle #5111 was visible from any RTS during this time window (see 
the visibility chart in figure 3). It was visible from Pogo, Ross, Cook and Hula(figure 
4). Allowing the maximum time permitted between supports (six hours) , the next 
.support should occure at 0920-1000 hours on June 16, 1986. This choice minimizes 
the total number of .supports required. However, the visibility .search routines 
determined # 5111 was not visible at any RTS at that time. As a result, the search 
routines found the support could be scheduled and still meet the six-hour constraint 
(Figure 4) 


Conclusion: 


Because of this research and prototype effort, the Air Force and its supporting 
aerospace contractors have increased their efforts to apply expert systems to 
scheduling .satellite siipporls. U.siiig a similar design, IBM Federal Systems Division 
is developing a full scale prototype ES called REGE ( for Request Generator) for the 
GPS PAPing process. REGE currently schedules five support types. By late 1988, 
IBM plans to have an operational version of REGE that will schedule 12 to 13 support 
types. When fully implemented for all 18 support types, REGE will reduce the PAP 
generation time by an order of magnitude. 

The success of this prototype demonstrates that ES technology can successfully 
automate the PAPing process. ESs are ideal for problems for which the solutions are 



obtained using experts’ heuristics and rules. However, for problems as complex as 
GPS PAPing process, extensive prototype programs are required to develop 
operational expert systems. Knowledge enginecr.'ng offers an appropriate 
methodology and the required software development tools. 

After going through the literature available in satellite operations scheduling area, in 
the next cliaptcr we present the development of an evaluation function for 'll’S 
operations scheduling. 



Chapter 3 

Derivation of the Objective Function to Help Resolve Visibility 
Clashes. 

A “visibility clash” occurs when one or more LEO satellites visible at a ground station 
have either overlapping visibilities or their respective LOS (Loss of Signal) and AOS 
(Acquisition of Signal) difference is less than the station reconfiguration time. The 
problem is to find the optimum support schedule, that maximizes the returns (or total 
value generated) from supports provided. This chapter proposes two separate 
objective 1‘unctions to help resolve visibility clashes optimally. 

.A clashing situation is difficult to precisely specify, as the patterns in which satellites 
moves arc enormous. A clashing situation as envisaged by ISTKAC (ISRO 
Telemetry, Tracking and Command Network) can have any of possible combinations 
ol' overlaps as shown in the following figures. 

3.1 Possible Visibility Clash Topologies: 

Fig 3.1 
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ho 3.2 NO l A 1 ION REFRESEN l ATlON OF A CLASH SITUATION 


Total visibility of all three satellites > 







3.2 The Notations used (Fig 3.2) in the Mathematical Model of a typical 
clash resolution problem (Kumar, Bagchi 2000) 

J - Difference of first clashing AOS and last LOS. 

I = Total number of visibility clasbes 

/ = 0 1 . 

P= Value function, to be maximiTied (optimised) 

«,= Start of visibility (AOS) of saLcllite i. 
hr End of visibility (LOS) of satellite /. 

.v,“ Start of support of satellite i. 

c^/= End of support ofsatellite i. 

min - Minimum time necessary for a support. 

max = Maximum time allowable for support. 

/•= Reconfiguration time. 

v/= Value contributed to P per unit time by supporting satellite /. 

.Y,= Support duration for visibility ofsatellite /. 

/, = Reconfiguration time is added to end of support of previous supported satellite. /, 
is a binary variable which reflects whether satellite / is supported or not. If r, =l it 
implies that satellite i was supported while /, =0 implies that the satellite was not 
supported. Thus, 

I, is defines as binary variable sucli that: 
if x, 0 ; tn j-O 


if x,>0 : 


ti \ l~ ^ 



3.3 An Evaluation Function for TTC Support Value Generation 

U is known that for ll'C operations the return per unit support time has direct 
correlation to the duration of support. From discussions held with ISRO experts we 
derived following points. 

• fl'C operations require a minimum support of 8 minutes (480 seconds),whenever 
support is provided. 

• Almost all operations can be done in this 8-minute span, 'fhus the returns obtained 
in supporting a satellite for minimum support time are high. 

• A support of more then 8 minutes is not necessarv', but if some additional 
visibility time is available, support may be extended beyond 8 minutes. Thus, 
more then 8 minutes support is desired but not required. 

I'hus we conclude that return function should be designed in such a way there is no 
return if a satellite is supported for a duration of less then 8 minutes but the return per 
unit time is high if support is exactly 8 minutes. However, the return per unit time 
diminishes after 8 minutes. We may explain it as follows. If a satellite is supported for 
10 minutes, the returns obtained are more for the “8”’ minute” as compared to the 
minute”, and returns in the “9‘'’ minute” are more then returns in the “10^'^ minute”. 

Considering these factors, an exponentially diminishing return function was designed, 
'fhis function meets all requirements listed above while incorporating a diminishing 
rate of return. A simple exponential function can be represented as 

(3.3.1) 

where 1 is time measured from origin. I'his exponentially diminishing function is 
shown in Figure 3.3. 

At t=0, f-X , thus X is the ordinate of the exponential curve at t=0. 


12 



Time ^ 


Fig. 3.3 


To mould this diminishing function to our requirement, we introduce one more term, 
B. 


(3.3.2) 

where B is a factor which will determine the actual decay rate of this exponential 
curve as time t increases. 


3.3.1 Determining value of B. 

In this problem, as we already indicated, the minimum support needed is 8 minutes. 
Similarly there is an upper limit also (16 minutes), after which the returns per unit 
time (or the value generated from support) becomes negligible. Based on the 
discussions with ISRO, we assumed that the function value might be assumed to 
reduce to 1% of its initial value (A.) when support time becomes 16 minutes. That is if 

( '180 socoiul.s 

/=l%of?. 

Therefore the value of B may be found as follows. 

Att-0, f=X 

Alt -480, /-l%of?t 

Hence 

which gives 

B=0.00958/A. (3.3.3) 


At an instant t, the total value of support provided up to t is the area covered between 
the exponential rate curve and the abscissa. We may determine this total value by 
integrating the exponential function with the integration going from 0 to t. 

f 

■fotal Value- Jic- ^"'dt , 

0 

which gives 

Total Value of Support = (1- (3.3.4) 


But in real life situation, the return function is not exactly exponential rather it has a 
discontinuous nature. This is shown in Figure 3.4 



Total returns (V) 


Total Return ={A+Exp} 
Fotal Return =0 


for t>480; 
for t<480; 


....(3.3.5) 


Where 
A=480*X 
Exp=(l- c^"V» 

B=0.00958/?l 

t=(Support Duration-480) 

The term A ensures high returns of the total support provided when a satellite gets a 
support of 480 seconds. 


Total Return as a function of / is shown in Figure 3.5 
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3.3,2 The Derivation of Lambda: 


The next challenge in the model is to find the value of X which in general would be 
different for every visibility. X is the rate of return when a satellite is supported for 8 
minutes ( minimum support for TTC operations). In other words, X is the height of the 
exponential rale curve at zero time for the pure exponential function. 

The utility of exponential decay logic depends on the derivation of X. As a result of 
discussions with ISRO experts, we concluded that X for a support (a combination of 
satellite, station and operation) depends on several factors as listed in Table 1. These 
factors each have their own impact on X. The relative impact of these factors is 
expressible as shown in terms of objectively stated priority and penalty based on the 
knowledge of LEO satellite technology. 


TABLE 3.1 FACTORS CONTRIBUTING TO SATELLITE SUPPORT 


S.No 

FACTOR 

VALUE 

PENALTY 

1 

Maximum elevation 

2 


2 

Satellite priority 

10 


3 

Exclusive pass (orbit over a single station) 

8 


4 

PL settings (Timer settings) 


10 

5 

Speeial TTC operations (must be done in that 
pass) 

10 


6 

Service (visibility gap) 


5 

7 

Minimum operations in the satellite 

6 


00 

Exactly ‘n’ ascending operations per satellite 


2 

9 

Exactly ‘m’ descending operations per satellite 


2 

10 

Controller constraint 


2 

11 

Ground station constraint 

2 



Combining these factors in an equation form gives X. The potential method of 
combining these factors reflects their nature and their effect on satellite support. 
Factors not considered in X may be used to guide the actual heuristic search for task 
allocation. 


A.=Max (FactorExclusive, FactorCyclic, FactorTimer, FactorSupportGap) 
+FactorSateUite+F actorEIevation. (3.3.6) 


FactorExclusive, FactorCyclic, FactorTimer and FactorSatellite are set on a scale of 1 
to 10. Other factors depend on dynamic parameters. FactorSupportGap depends on 
the time gap between current visibility and previous supported visibility for the same 














satellite. I'uetor elevation is step iunction assumed for illustration to be rising from 0 
(for elevations up to 5°) to 10 (for elevations >= 60°). 


Note that the above expression for X has two parts. The first part is the maximum of 
lour factors while the other part added to it. It is noticeable that the llrst part has 
terms that correspond to the actual satellite operation . FactorSatellite and 
FactorElevation correspond to satellite concerned and the quality of the visibility 
window and are independent of the actual operations. Further, the first part has factors 
that are very important in support decisions; if any one factor among these is high, the 
visibility should be su|iporled. I'his requirement is ensured by taking maximum of 
Factor Exclusive, FactorCyclic, FactorTimer and FactorSupportGap. 

3.3.2. 1 Effect of Parameters of A 

FactorExclusive: This factor assumes values of either 0 or 10. If a visibility is 
exclusive in that specific orbit, this factor becomes 10. Otherwisti it is 0. Thus if a 
visibility is exclusive, this factor increases the value of X so that the likelihood of 
support for this visibility rises. 

FactorCvclic: This factor also takes values of 0 or 10. If a cyclic operation falls in this 
slot of visibility, this factor becomes 10 otherwise it is 0. Thus if a satellite requires a 
cyclic operation in this particular visibility, this factor becomes 10 and thus increases 
the value of X and hence the likelihood of providing support for this visibility 
becomes high. 

FactorTimer: This factor also takes values of 0 or 10. If a Timer operation ( required 
by a subsequent payload (PL) operation) falls in this slot of visibility, this factor 
becomes 10, otherwise it is 0. 'fhus by giving ‘Timer Operation’ a value of 10 we 
ensure that the subsequent PL operation is serviced as timer setting is prerequisite for 
a PL operation on certain satellites. 

f^actorSupportGap: Generally speaking, a satellite has to be supported throughout the 
day. An important consideration in satellite support scheduling is to keep the gap 
between supports small. factorSupporlGap ensures that a satellite will have an 
increased likelihood of support after a certain time gap. The value of this factor 
depends on the time gap between eurrent visibility and previous supported visibility 
Ibr the same satellite. As this gap increases, FactorSupportGap also builds up. After 



some gap this will build to such a high value that it will exceed any other factor in 
value. Thus such a a satellite will receive a support. 

FactorSatellite: This factor is a satellite-dependent factor and is independent of 
operations and visibility. This factor is additive term in X. This factor indicates that 
when a conflict develops, more important satellite (as indicated by network 
management) should have a higher probability of getting a support. 

FaetorF.levation: The Factor Elevation has a step function. This factor also is additive 
in X. This indicates that a higher elevation visibility is more suitable for doing an 
operation as compared to low elevation visibility. 

3.4 Evaluation function for Payload Operations Support 

For Fayload operations, the return per unit time is constant; hence the evaluation 
runclinn becomes linear. U is (hen simply the product of the support time (Xj) and 
returns per unit lime (Vj), beyond any minimum support required. 

F l!overalli{Vi*Xi} (3.4.1) 

T he function is shown in fig. 3.6 




3.5 The Mathematical Model for TTC Support 


Maximize, 

P=Zoveralli{Ai+((l-i''"'')/B)} (3.5.1) 

Subject To: 

i) Start of support of satellite must be at AOSi or later and it should be equal to or 

less then LOSi. 

bi>Si>ai (3.5.2) 

ii) luul of support of satellite must be at LOSj and it should he equal to or more 


than A O.S). 

ai<ei<bi (3.5.3) 

iii) Reconfiguration allowance. 

Si> { Maximum (c^y^ tk)]T r (3.5.4) 

k (I I) I 

iv) Duration of support of satellite i. 

Xi.ei-Si (3.5.5) 


v) Constraint for minimum time of support. This may be either 0 or greater than min 

Xi=0 or Xi>min (3.5.6) 

vi) Constraint for maximum time of support. This should be less than max. 

Xi<max (3.5.7) 

vii) Non negativity constraint 


Si, ei > 0 ; ; 



'I’he objective function 3.5.1 and the presence of constraint 3.5.6 makes the problem 
non-linear. 

In the above formulation, for convenience and without loss of generality we made the 
assumption that reconfiguration timings are constant and have same value ’r’ for all 
set of ground stations and satellites. 

The above objective functions, (3.3.5) and (3.4.1) will be used in the subsequent 
chapters to actually resolve visibility clashes optimally. (3.3.5) models the objective 
for TTC visibility support while (3.4.1) models the objective for a PL support. 



Chapter 4. 

Greedy search Approach for Optimal Resolution of Visibility 
Clash 

A Greedy Search or Neighbourhood Search is a local search algorithm, which is based 
on iterative improvement. The application of a greedy search algorithm presupposes 
the definition of configurations, a cost function and a generation mechanism, i.e. a 
simple prescription to generate a transition from a configuration to another one by 
small perturbation. “Configuration” means a set of values of decision variables. The 
major disadvantage of this algorithm is that it ends up in a local optimum solution. 
Therefore, in combinatorial optimisation problems, some variation of greedy search is 
used to skip out from local optima. The major advantages and disadvantages of this 
algorithm are discussed in more detail in section 4.1.1 of this chapter. 

In this chapter we refer closely to the mathematical formulation of a sample global 
oplimisalion problem encountered in I.IT) satellite support scheduling. 

4.1 Approach using Greedy search 
4.1.1 Greedy Search: 

1 he simplest type of search in oplimi/.alion is loettl, greedy .search, in which we start 
w ith a solution X. At each step of the search, we choose the neighbor of X, X{n)i 
such that X(n)i is lower than the cost of any other neighbor and is also strictly 
smaller (or larger, if a maximum rather than minimum is desired) than any previously- 
encountered cost. The algorithm stops when no such neighbor exists. This approach is 
also known as hill climbing (when searching for a maximum) or descent (when 
searching for a minimum). In the case of a complicated objective function with 
multiple modes, this type of search is likely to reach a local optimal point, so one 
common variation is to repeat the procedure with a number of starting positions and 
take the best result. 

The disadvantages of the Greedy Search method are as follows (ref). 

• By definition, Greedy Search algorithm terminates in a local optimum and 
there is generally no information as to the amount by wliich this local 
optimum deviates from a global optimum. 

• The obtained local optimum depends on the initial configuration or starting 
point, for the choice of which generally no guidelines are available. 



• In general it is not possible to give an upper bound for the computation time of 
greedy seaarch. 

However, this algorithm has the advantage of being broadly applicable: 
Configurations, a cost function and a generation mechanism are usually easy to 
doline. llcsidcs, though upper bounds for computation times are missing, a single run 
of an greedy algorithm (for one initial configuration) can on the average be executed 
in a small amount of computation time. 

fo avoid some of the aforementioned disadvantages, on : might think of a number of 
common modifications to the basic greedy search. 

1 . Execution of the algorithm for a large number of initial configurations, say N 
(e.g. uniformly distributed over the set of configurations R) at the cost of an 
increase in computation time; for N oo, such an algorithm finds a global 
optimum willi a probability of 1, ifonly for the fact that a global minimum is 
encountered as an initial configuration with probability 1 as N->oo; 

2. (Refinement of 1.) Use of information gained from previous runs of the 
algorithm to improve the choice of an initial configuration for the next run 
(this information relates to the structure of the set of configurations); 

3. Introduction of a more complex solution generation mechanism (or, 
equivalently, enlargement of the neighborhoods), in order to be able to ‘jump 
out’ of the local optima corresponding to the simple generation mechanism. 

4. Acceptance of the non-improving transitions in the cost function in a limited 
way. 

The second and third approaches aue problem-dependent and hence their applicability 
is limited. The first approach is the traditional way to solve combinatorial 
optimization problems approximately by greedy search. 

4.2 Initial Solution Generation Mechanism: 

To tackle the satellite clash resolution problem, we have repeated the greedy search 
procedure with many different starting solutions. In this problem context, if we look 
at the search space more closely, it would resemble many hills distributed at various 
locations in the whole search space. The peculiarity of this search space is that these 
various hills have no continuity between them, and hence while searehing, one has to 
sometimes jump from one hill to another in order to come out of local optimum. 



Furthermore, even if we consider a single hill, it may also have more than one peak 
and hence search becomes even more complex. 

Thus in this problem the challenge would be to select suffieiently many different 
starting solutions so that no part of the search space would remain unexplored. We 
took the advantage of the fact that, a satellite, if supported, has to be given at least a 
minimum time of support. Otherwise the satellite is said to be unsupported and if 
unsupported, the support time is zero. Therefore we generated 2"-l different binary 
combinations of possible support scenarios where n = number of satellites clashing 
over a station, fhe combination.s generated are simply in accordance with the binary 
numbers where a ‘0’ represents that the particular satellite is “not supported" and M" 
represents that the satellite is “supported” with minimum support allocation. This is 
explained below with an example. 

Suppose ?i satellites arc clashing over a station. Let’s call them satl, sat2, sat3. Now 
we would generate 2‘’-l =7 different binary combinations as shown below. 

001, 010, on, 100, 101, no, in. 

Note that 000 ( No support to any satellite ) is not an interesting case, hence it is 
(unitted. Now the first binary number, that is 001, indicates that, in the first initial 
solution generated, satl, sat2 will not be supported as there is ‘0’ in the first two 
places of the binary number, and only sat3 will be supported as the binary number has 
n ' at the third place. Similarly, the second initial solution generated will have sat2 
supported and ‘no support’ for satl and sat3. Similarly 7 different initial solutions will 
be generated with final solution supporting all the three satellites. 

One way to further reduce the number of initial solutions is to consider only those 
initial solutions, in which satellites having maximum weight factors clash. This is 
because, after experimentally running the algorithm for all the initial solutions we 
found that, most of the time, the best solution obtained contains the support for the 
satellites with maximum value of the weigh factors {vj}. Hence if we consider only 
those initial solutions in which the satellites having maximum weight factor clash, we 
would be able to reduce the computation time without likelihood of loosing the best 
solution existing. 

As already noted, there is one major problem with the greedy algorithm. Even though 
we would generate one initial solution for one hill in the solution space, some hills 
have multiple peaks and the greedy algorithm may end up in a local optima. Hence we 



may have to consider more initial solutions than 2" . We will consider this case with 
an example. 

Suppose satl , sat2 and sat3 are supported in one of the initial solutions(case of ’ 1 1 1’). 
'I'lie situation is like this. 

• • Satl 

vsl vel 

• • Sat2 

vs2 ve2 

• • ^ ^ 

Sat3 

vs3 ve3. 

Time 

► 


Where, 

vsi =visibility start time of i‘'’ satellite. 

til 

vei = visibility end time of i satellite 

In this situation we have to generate two initial solutions. The first one will have sat2 
supported from vs2 side, that is, the minimum support allocation is done from vs2 
side and hence. 


Allocation start time of sat2 = max( ( allocation end time of satl+station 
reconfiguration time),(visibility start time of sat 2)) (4.2.1 ) 


Allocation end time of sat2 - Allocation start time of sat2+ minimum support time... 




Non-Allocated Window 

wmmmmmmmmmmmm Allocated Window. 

R - Station reconfiguration time. 

M - Minimum support time. 


The second situation will be like this- Satl and Sat3 will be allocated windows as 
above and sat2 will be allocated a window from ve2 side. Now the equations become, 

Allocation end time of sat2 = min( ( allocation start time of sat3-station 
reconfiguration time), (visibility end time of sat 2)) (4.2.3) 


Allocation start time of sat2 = Allocation end time of sat 2- minimum support time. 

(4.2.4) 



The figure is shown next. 



Non allocated window. 
Allocated window. 


R - Station rcconliguration time. 

M - Minimum support time. 

'I'hese two initial solutions will both correspond for the ‘ 1 1 1’ situation. 

Mence the modified formula for number of initial solutions is. 

Number of initial solutions = 2" + Z!i<=n {("Ci * 2^'‘^V{"Ci)} (4.2.5) 

The algorithm proceeds as follows: 

Generate initial solutions. 

Repeat 

Take first initial solution. 

Repeat 

Make this initial solution as current solution. 

Find all feasible .solutions in the neighbourhood of the current solution. 
Evaluate all of them. 

Find the best solution amongst them. 

Make this as current solution. 

Until (no improved solution is found or the best solution remains same for a 
number of iterations) 

Until (all initial solutions are checked.) 

Store the best value obtained. 



Greedy Search Flow Chart 



generate initial 
solution 



select the best solution from 
the feasible neighbourhood 
such that this solution is 
better than curent solution. 


mark the selected 
solution as current 
solution 







4.3 Implementation and results obtained: 


We implemented this algorithm using JAVA. The code is given in Appendix. We 
tested the algorithm on various clash situations. Here is result: 

Four satellite visibilities are clashing and their visibility start timings and visibility 
end timings are as shown follows- 

Vis_start_sall =6520 
Vis_start_sat2 = 7244 
Vis_start_sat3 = 7712 
Vis start_sat4 = 8027 

Case A: Linear objective function 
F = Zovcralli{(ei-Si)*Vi} 

Case B: Non Linear objective function. 

F = lover all i{Arf(l-e-<‘^"'’'>*'’‘’‘''' * Vi)/.00958} 

Ai= minimum support time * vj 
Cj = end of support of satellite i 
Si = start of support of satellite i 

Vj - value contributed to I’ per unit time by supporting satellite i 


Vis_end_satl =7115 
Vis_end_sat2 = 8318 
Vis_end_^sat3 = 8650 
Vis end sat4 = 9]96 



Table 4.1 Results obtained using Greedy Search 


Problem 

Satellite 

Vj 

Case 

A: Linear 

objective 

Case 

B: 

Non-linear 

Scenario 

Number, 


function 



objective function. 



i 


Si 

Si 

Xi= 

Function 

Si 

ei 


Function 






ei- 

Si 

Value 



1 


1 

1 

1 

6520 

6520 

0 







2 

2 

7244 

7613 

369 


7244 

7884 

640 



3 

3 

8213 

8213 

0 


7712 

7712 

0 



4 

4 

8213 

9196 

983 


8484 

9196 

712 


' 

2 

1 

4 

6520 

7115 

595 

4467 

6520 

71,15 

595 

3855.55 


2 

3 

7715 

8318 

603 


7715 

8318 

603 



3 

2 

7712 

7712 

0 


7712 

7712 

0 



4 

1 

8918 

9196 

278 


8027 

8027 

0 



3 

1 

1 

6520 

6520 

0 

4670 

6520 

6520 

0 

3415.99 


2 

4 

7244 

8227 

983 


7244 

7956 

712 



3 

3 

7712 

7712 

0 


7712 

7712 

0 



4 

2 

8827 

9196 

369 


8556 

9196 

640 



4 

1 

1 

6520 

7112 

592 

4344 

6520 

7115 

595 

2881.89 


2 

3 

7712 

7712 

0 


7712 

7712 

0 



3 

4 

7712 

8650 

938 


7715 

8650 

935 



4 

2 

8027 

8027 

0 


8027 

8027 

0 




4.4 Graphs: 

4.4.1 Conversion Graph of Greedy Algorithm 

Conversion Graph for Greedy Search 
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4.4.2 Satellite Support seenario Graphs 


Linear Ohjecl ive function: 

Problem case 1 : vl=l v2=2 v3=3 v4=4 



443 



Problem case 2: vl=4 v2= 3 v3= 2 v4=l 










Problem case 4: vl=l 


v2=3 


v3=4 


v4=2 



6520 7020 7520 8020 8520 9020 9520 

■ ; \ . 'I : ’ 


Non-Linear Objective Function. 


Problem case 1: vl =1 v2=2 v3 = 3 v4=4 

Satellite Support Scenario with Non-linear 
objective function, Case 1 



6520 7020 7520 8020 8520 9020 9520 

Time 



Satellites I: Satellites 


Problem case 2: vh'4 v2“ 3 v3=2 v4=l 

Satellite Support Scenario with Non-linear 
objective function, Case 2 
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Time 


cm case 3; vl = l v2=4 v3=3 v4=2 


Satellite Support Scenario with Non-linear 
objective function, Case 3 
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Problem case 4; vl=l 


v2= 3 v3=4 


v4=2 


Satellite Support Scenario with Non-linear 
objective function, Case 4 
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Chapter 5 

Approach using Tabu Search 

5.1 Tabu Search 

Tabu search is an iterative procedure for solving discrete combinatorial optimisation 
problems. It was first suggested by (Jlovcr |7) and since then has become increasingly 
used. It has been successfully applied to obtain optimal or sub-optimal solutions to 
such problems as scheduling, timetabling, travelling salesperson and layout 
optimisation. 

The basic idea of the method, described by Glover [7], is to explore the search space 
of all feasible solutions by a sequence of moves. A move from one solution to another 
is the best available. However, to escape from locally optimal but not globally optimal 
solutions and to prevent cycling, some moves, at one particular iteration, are classified 
as forbidden or (ahu (or taboo). Tabu moves are based on the short-term and long- 
term history of the sequence of moves. A simple implementation, for example, might 
classify a move as tabu if the reverse move has been made recently or frequently. 
Sometimes, when it is deemed favourable, a tabu move can be overridden. Such 
aspiration criteria might include the case that, by forgetting that a move is tabu, leads 
to a solution which is the best obtained so far. 

I’ormally, tabu search is applied to an optimisation problem as Ibllows. Suppose F is 
the real-valued objective function on a search space S and it is required to find a ceS 
such that F( c) has maximal value. For combinatorial hard (NP-complete) problems, 
this requirement needs to be relaxed to finding ceS such that F(c) is close to the 
maximal value (a snh-optimal value). This is because any known algorithm to 
determine the maximal solution requires time that is exponentially increasing with the 
problem si/e ( 1 7]. Sub-optimal problems may be solved by halting when a certain 
threshold for an acceptable solution has been achieved or when a certain number of 
iterations have been completed. 

A characterisation of the search space S for which tabu search can be applied is that 
there is a set of k moves M={ml,..mk} and the application of the moves to a feasible 



the variables changed during a move is taken as the attribute and used for enforcing 
tabu status. 

The neighbourhood structure we are using is very simple. We add and subtract step 
size value to the current solution and get the neighbourhood solution. This we call as 
local neighbourhood. In our global neighbourhood we try to make a big jump in the 
search space through divcrsillcation. Wc use problem specil'ie knowledge to define a 
number of diversified solutions. 

We diversify our search to a new region when we find no improving solution in the 
current region for a number of moves. After jumping to the new search area we 
intensify our search in that area. The main advantage in this particular problem is that 
we already know which are the most promising areas and can limit our search to that 
only. 

5.2.1 Results obtained 


We tested our algorithm on various clash situations. One of them is as shown below. 
Four satellites are clashing and their visibility start timings and visibility end timings 
are as shown follows- 

Vis_start_satl = 6520 Vis_end_satl =7115 

Vis_start_sat2 = 7244 Vis_end_sat2 = 8318 

Vis_starl_sat3 = 7712 Vis_cnd_sat3 = 8650 

Vis start sut4 8027 Vis end sat4 = 9I90 


Case A: Linear objective function 

F “■ Sover all i {(6i”Sj)*V|} 

Case B: Non Linear objective function. 

f == lover all i { Ai+( 1 * Vi)^^ 

Aj= minimum support time* j 
C| end of support of satellite i 
s, = start of support of satellite i 
Vj = priority weight of satellite i 










Table 5.1 Results obtained using Tabu Search: 


Problem 

Scenario 

Satellite 

Number, 

i 

Wi 

Case A: Linear objective 

function 

Case B: Non-linear objective 

function. 


Si 

Cj 
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ei- 

Si 

Function 

Value 

Si 

Ci 

xr 
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Si 
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Value 

1 

1 

1 

6520 

6520 

0 

4670 

6520 
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3481.79 

2 

2 

7244 

7613 

369 

7600 
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3 

3 

8213 
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0 

7712 
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0 

4 

4 

8213 
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983 

8680 

9196 
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2 

1 
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7115 

595 

4467 

6520 

7033 

513 

3962.04 

2 

3 

3 

2 

7715 

7712 

8318 

7712 

603 

0 

7633 

7712 

8116 

7712 
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4 

1 

8918 
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8716 

9196 

480 
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1 

1 

6520 
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0 

4670 

6520 
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3481.79 

2 

4 

7244 

8227 

983 

7600 
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3 

3 

7712 

7712 

0 

7712 

7712 

0 


4 
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8827 

9196 
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8716 

9196 
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4 

1 

2 

1 

3 

6520 

7712 

7112 

7712 

592 

0 

4344 

6520 

7600 

7000 
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2971.34 

3 

4 

7712 

8650 

938 


7712 

7712 

0 


4 

2 

8027 
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0 


8716 
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480 






5.3 Graphs 

5.3.1 Conversion Graph for Tahu Search: 
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5.3.2 Satellite Support Scenario Graphs: 

Linear Objective function: 


Problem case 1 : vl=l v2=2 v3=3 v4=4 





Linear Objective Function, Case 1 
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Problem case 2: vl - 4 


v2= 3 
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Problem case 4; vl=l 
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Non Linear Objective Function: 

Problem case l:vl=l v2=2 v3=3 v4=4 

Satellite Support Scenario with Non-linear 
objective function 
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Chapter 6 

Approach using Simulated Annealing 


6.1 Simulated Annealing: 

Siiiuilatcd aiiiicaliiig is a numerical optimisation tccliniquc based on the principles of 
thermodynamics. Annealing refers to the process in which a solid material is first 
melted and then allowed to cool by slowly reducing the temperature. The particles of 
the material attempt to arrange themselves in a low energy state during the cooling 
process. The collective energy states of the ensemble of particles can be considered 
the "configuration" of the material. The probability that a particle is at any energy 
level can be calculated by use of the Boltzmann distribution. As the temperature of 
the material decreases, the Boltzmann distribution tends toward the particle 
configuration that has the lowest energy. Metropolis et al [ll] first realized that the 
thermal equilibrium process could be simulated for a fixed temperature by Monte 
Carlo methods to generate sequences of energy states [1 1], The s}'.stem is perturbed to 
yield a new configuration of the particles. The energy level before perturbation (Fs) 
and the energy level alter perturbation (Ei) are compared. If Es is greater than Ei (i.e., 
E > 0), the new perturbed system is accepted as the new configuration of particles. If 
\i > 0, the probability of accepting the perturbed system Ibllows the Metropolis 
criterion shown in following equation 


p = exp(-E/kT) 


where k is the Boltzmann constant and T is a fixed temperature. Using this criterion, 
the material will eventually reach its equilibrium configuration. 

Ibis basic concept can be applied to numerical optimisation problems. Simulated 
annealing (SA) applies the Metropolis criterion to a series of variable settings 
(configurations) for the system being optimised. The new variable settings are 
obtained by perturbing the current configuration and can' be thought of as steps or 
movements on the response surface. For numerical optimisation, the response 
(objcelivc) I'unction score (R) replaces the energy terms. The concept of temperature 
is retained. However, it is now combined with k and used as an important control 



factor. The probability of accepting a detrimental step (i.e., R > 0, assuming 
minimization) is governed by following equation 

p ■ cxp(-R/ 1) 

A random number, P, is drawn from a uniform random distribution on the interval 
[0,1]. If P > p, the detrimental step is rejected and a new step taken from the current 
position. If p< P, the detrimental step is accepted and the new configuration replaces 
the old one. A new step is then taken relative to this configuration. This criterion 
allows the possibility of a new variable being accepted as the new configuration even 
when it has a worse response function score than the current configuration. Moves 
that are very poor (i.e., large R) are less likely to be accepted than moves that are not 
poor. This feature allows the algorithm to "walk" from local optima. New steps are 
taken until some termination eriterion is reached. 

Analogous to the physical process, T is slowly reduced causing the probability of 
accepting a poor move to decrease with time. The schedule by which T is reduced is 
called the cooling schedule and is critical to the success of SA. Two important 
parameters governing the cooling schedule are the step size for perturbation and T. 
The iiaraineter settings suggested by most researchers arc step sizes and f values that 
allow approximately 80% of the poor moves to be accepted [12], 

The primary advantage of SA is the ability to move from local optima. Thus, the 
ability to lind the global optimum is not related to the initial conditions (i.e., the 
starting point). SA is also very simple to implement. The primary disadvantages to SA 
are the subjective nature of choosing the SA configuration parameters (c.g., f and 
step size) and that it typically requires more response or objective function 
evaluations than other optimisation approaches. SA has been termed a "biased random 
walk". Unlike other optimisation approaches that attempt to make intelligent moves 
on the response surface, the steps in SA are taken randomly. Thus, SA is classified as 
having a weak search heuristic. The fact that SA employs no knowledge of the 
response surface can be an advantage or a disadvantage depending on the application. 



Chart 




6.2 SA on this problem 


The basic algorithm remains same. As we have do deal with discontinuous nature of 
search space, we have to have some mechanism to make a jump in the search space. 
One way is to define the neighbourhood structure in such a way that the selected 
solution from neighbourhood automatically jumps to the unexplored region. This 
makes the number of neighbourhood solutions very large that is of the order of 
laisc to ?.''rai.so to 2'’...upto (nl 1)), which i,s very clifflciill to construct and not 
practical. 

In our case the problem structure is such that we can jump in the search space if we 
define some initial solutions already; its equivalent to starting algorithm with different 
initial solutions in greedy search. And here the advantage is that we can have only 2" 
different solutions at hand from where we can start. 

As mentioned earlier, in this problem, the individual hills in the search space have 
different peaks, but SA can skip away from local optima and find the global optima in 
individual hill. 


6.3 Implementation and Parameter optimisation. 

The critical parameters in simulated annealing are, the initial high temperature and the 
fraction by which initial temperature is brought down to freezing point. This is called 
cooling schedule. We tried to find out best value of T (Initial Temperature) and a 
(fraction by which f is reduced). 

Results: 

We considered the following sample problem; 


Vis_start_satl = 6520 
Vis_start_sat2 = 7244 
Vis_start_sat3 = 7712 
Vis_start_sat4 = 8027 


Vis_end_satl =7115 
Vis_end_sat2 = 8318 
Vis_end_sat3 = 8650 
Vis_end_sat4 = 9196 


Case A: Linear objective function 

f = Sover ;illi{(Ci-Si)*Vi} 

Case B: Non Linear ohiective function. 

1' love, all i{ Ai * Vi)/.00958} 



where, 

Ap initial support time * Wj 
Ci = end of support of satellite i 
Si = start of support of satellite i 
Vj ==- priority weight of satellite i 

Table shows the average and standard deviation of the best value obtained after five 
different seeds given to random number generator. 

Table 6.1 Average of the Best function value Generated for different 
values of T & a using different seeds for random number generator. 

/. Linear objective function. 


a t 

► 

10 

15 

20 

25 

30 

0.5 

Avg. 

4344 

3602 

4344 

4344 

4344 


SI) 

0 

0 

0 

0 

0 

0.7 

Avg. 

4344 

4344 

4344 

4344 

4344 


SD 

0 

0 

0 

0 

0 

0.9 

Avg. 

4344 

4344 

4344 

4344 

4344 


SD 

0 

0 

0 

0 

0 

0.99 

Avg. 

4344 

4344 

4344 

4344 

4344 


SD 

0 

0 

0 

0 
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2. Non-linear objective function. 




10 

15 

20 

25 

30 


Avg. 
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2764.77 

. 



2762.56 




4.633 

4.633 

4.805 

5.085 

0.7 

Avg. 

2764.77 

2769.89 

2764.77 

2769.89 

2769.89 


SD 


4.536 

4.633 

4.536 . 

4.536 

0.9 

Avg. 

2769.89 

2769.89 

2769.89 

2971. .34 

2971.34 


SD 

4.536 

2.609 

2.609 

0 

0 

0.99 

Avg. 

2769.89 

2971.34 

2769.89 

2971.34 



SD 

2.609 

0 


0 

0 













Though there is no effect of seed on best value obtained when using linear objective 
function, above table ( with non-linear objective function values) clearly indicate that 
higher llie initial leinperalurc and finer the fraction better is the result. Hence wc 
select the initial temperature value T as 25 and fraction a = .99. 

I'ollowing fable shows the best objective function value obtained by setting T=25 
and a == .99. 

Table 6.2 Results obtained using Simulated Annealing 
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6.4 The Graphs: 

6.4.1 Conversion Graph for SA. 


Conversion Plot for SA 
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6.4.2 Satellite Support Scenario graphs 

Linear Objective function: 

Problem case I : vl=l v2=2 v3=3 v4=4 





Problem case 2: vl = 4 v2= 3 v3= 2 
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Problem case 3: vl=l v2=4 v3=3 v4=2 ' 
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Non Linear Objective Function: 

Problem case 1; vl = l v2=2 v3=3 v4=4 


Satellite Support Scenario with Non-linear 
objective function 
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Chapter 7 

Approach using Genetic Algorithm 

7.1 Genetic Algorithm 

Genetic algorithms are search algorithms based on the mechanics of natural selection 
and natural genetics. I hey combine survival ol the littest among string structures with 
a structured yet randomised information exchange to form a search algorithm with 
some of the innovative (lair of human search. In every generation, a new set of 
artificial creatures (strings) is created using bits and pieces of the fittest of the old; an 
occasional new part is tried for good measure. While randomised, genetic algorithms 
are no simple random walk. They efficiently exploit historical information to 
speculate on new search points with expected improved performance. [ 8 ] 

7.2 GA on this Problem 

Kumar and liagchi (2000) solved this visibility clash resolution problem using GA. 
We describe here the approach given by them. 

7.2.1 GA Chromosome Construction 

'I'hey use binary coding to code the chromosome because in this problem the decision 
variables {xi} are discrete numerical quantities. A single chromosome would contain 
all information for the full duration of support for each satellite involved in a clashed 
set of visibilities at a station. 

For example: 

At a station 4 satellite visibilities are clashing, then the chromosome would contain 
values for Xi, Xj, X 3 , and X 4 . 

A typical chromosome would look as below. 


dlOOO 10000110 10100 0101010 0011110 


X],/ [-bit 


X 2 ,/ 2 -bit 


X 3 ,/ 3 -bit 


Xn,/n-bit 


7.2.2 Deciding the Length of the Chromosome 


In GA lonmilalion using binary coded variables, the length of a substring representing 
real a variable depends on desired accuracy of that variable. If one do f^-bit coding 
for a variable, the obtainable aecuracy in that variable is approximately 

In ISRO's case the upper limit on operation timing is 16 Minutes 23 Second. Thus the 
attainable accuracy with a 10- bit substring for all Xj’s will be; 

( 1 6*60+23 V2 =0.959 seconds 

this level of |irecision will always keep the error to be less then one second if one 
selects =10 bits. 

7.2.3 Feasibility Assurance 

The feasibility of a species is maintained by incorporating constraints inside 
chromosome coding. This is done by taking modified upper and lower limits of 
|X)ssible values of time. 

I'or example, consider decision variable x# . 

The limits for x/ are [ai , bi] (Fig 3.2) 

And limits for X 2 will be from 

[(Greater of ei+r and aa) and b 2 J 
the above approach would ensure feasibility of all chromosomes. 

7.2.4 GA Operators 

7.2.4. 1 Crossover 

In the crossover operation, new strings are created by exchanging information among 
strings of the mating pool . In this application they u.sed single point cros.sover. 



The Crossover used here is one-cut point and it exchanges the right part of two 
parents to generate offspring . The single point crossover operator is performed by 
randomly choosing a crossing site along the string and by exchanging all bits on right 
side of crossing site between the two parents . As applied to this problem, a typical 
crossover is shown as below. 


Crossover site 

A 

Parent l-> 1 0 1 0 1 0 1 OOOOOOO 111101010110 
l’aren(2-> 11000010101000 111010101001 

Chi 1 d 1 1 0 1 0 1 0 1 0000000 111010101001 
Child2^ llOOOOlOlOlOOOl 11 1010101 10 


7. 2. 4. 2 Mutation 

Mutation alters one or more genes with a probability equal to the mutation rate . So, if 
the probability of mutation is set at 0.01, so we expect that, on average, 1% of total bit 
of the population of solution would undergo mutation. 


7.2.5 Fitness Function 

For calculation of fitness function (Value function), first the decoded values of 
chromosome string are calculated by a linear mapping rule. The rule used is: 


y Un _ (/.) 

X; -f— * — — "-i — decodedvalue{s i) 


-1 


Mere decoded value {Xi} is the integer conversion of binary string [sj}. 
1 he decoded value of binary string is calculated as: 



f-l 

Decoded value (8,)= X2'.v, Where Si 6(0,1) 

I 0 

Al ter all Xi’s are calculated, total value function of the chromosome may be calculated 
by Eq.3.4.1 as: 

P = lover all i{Vi*Xi} 

Or by Eq. 3.3.5 as: 

P = lover all * Vi)/.00958} 

7.2.6 Selection 

Selection strategy used in the their work is Elitist selection. Elitist selection ensures 
that the best chromosome is passed into the new generation. In this strategy all 
individual fitness values are scanned and best member is selected. Number of 
selections made is such that the size of population remains constant. 

7.3 Implementation and Results 

They tested the algorithm on a real life situation problem of clashing visibilities. The 
visibilities are shown below. 

Vis_start_satl = 6520 Vis_end_satl =7115 

Vis_start_sat2 = 7244 Vis_end_sat2 = 8318 

Vis_start_sat3 = 7712 Vis_end_sat3 = 8650 

Vis_start_sat4 = 8027 Vis_end_sat4 = 9196 

Case A: Linear objective function 
Maximize, 

F = Zovcralli{(ei-Si)*Vi} 

Case B: Non Linear objective function. 

Maximize, 

F = lover all i{Ai+(l-e-‘‘=‘-^‘^' “‘'-"' * Vi)/.00958} 

Ai= minimum support time * Vj 
Ci = end of support of satellite i 



Decoded value (Si)= ^ 2' s, Where Sj e (0,1 ) 

I 0 


Alter all Xj’s are calculated, total value lunction of the chromosome may be calculated 
by Eq. 3. 4.1 as: 

P = lover all i{ Vi* Xi} 


Or by Eq. 3.3.5 as: 

P = Eoveraii * Vi)/.00958} 

7.2.6 Selection 


Selection strategy used in the their work is Elitist selection. Elitist selection ensures 
that the best chromosome is passed into the new generation. In this strategy all 
individual fitness values are scanned and best member is selected. Number of 
selections made is such that the size of population remains constant. 

7.3 Implementation and Results 


1 hey tested the algorithm on a real life situation problem of clashing visibilities. The 
visibilities are shown below. 


Vis_start_satl =6520 
Vis_start_sat2 = 7244 
Vis__start_sat3 = 7712 
Vis_start_sat4 = 8027 


Vis_end_satl =7115 
Vis_end_sat2 = 8318 
Vis_end_sat3 = 8650 
Vis_end_sat4 = 9196 


Case A: Linear objective funciion 
Maximize, 


F = Eoveralli{(ei-Si)*Vi} 

Case B: Non Linear objective function. 
Maximize, 

F = lover all i{Ai+(l-e-<'=‘-'’'>* "‘’‘'-"' * Vi)/.00958} 


Aj- minimum support time * Vj 
Ci = end of support of satellite i 



Si = start of support of satellite i 

Vi = value contributed to F per unit time by supporting satellite i 


Table 7.1 


Problem 

Scenario 

Satellite 

Number, 

i 

Wi 

Case A: Linear objective 

function 

Case B: Non-linear objective 

function. 


Si 

Cj 

xr 

Ci- 

Si 

Function 

Value 

Si 

Cj 

Xr 

Ci-Si 

Function 

Value 

1 

1 

1 

6520 

6520 

0 

4670 

6520 

7000.01 ■ 

483.4 

3478.67 

2 

2 

7244 

7613 

369 

7600.01 

8082.5 

480.4 

3 

3 

8213 

8213 

0 ' 

7712 

7712 

0 

4 

4 

8213 

9196 

983 

8682.5 

9196 

512.2 

2 

1 

4 

3 

6520 

'7715 

7115 

595 

4467 

6520 

7023.^5 

5015 

3958.48 

2 

8318 

603 

7623.55 

8115.62 

492.1 

'3 

"2 

7712 

7712 

0 

7712 

7712 

0 

4 



8918 

9196 

278 

8715.62 

9196 

480.3 


3 

1 

1 

6520 

6520 

0 

4670 

6520 

7000.1 

480.1 

3481.01 

2 

4 

7244 

8227 

983 

7600.1 

8115.8 

515.7 

3 

4 

3 

¥ 

7712 

8827 

7712 

0 

369 

7712 

7712 

0 

9196 

8715.8 

9196 

480.1 


4 

1 

1 

6520 

7112 

592 


6520 



2971.01 

2 

3 

7712 

7712 

0 


7600.07 

8115.01 

516 


3 

"T” 

2 " 

7712 

8650 

938 


7712 

7712 

0 


4 

8027 

8027 

0 


8715.01 

9196 

480 



flic parameters used for (lA in above case are: 


I’m™ Probability of Mutation = O.l 


Pc Probability of Crossover “ 0.95 










Population Size = 50 
Generations = 5000 

7.4 Graphs 

7 . 4.1 Conversion Graph for GA : 


Conversion Graph of GA 
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7.4.2 Conversion Graph for different values of seed: 


Conversion of GA for different values of seed 
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7.4.3 Satellite Visibility support Graphs for GA: 

Linear ( )hjeciive function: 

Problem case 1 : vl-1 v2~2 v3=3 v4=4 


Linear Objective Function, Case 1 
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Problem case 3 : vl==l 


Linear Objective Function, case 3 
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J^on Linear Objective Function: 

Problem case l:vl=l v2=2 v3=3 v4= 

I Satellite Support Scenario with Non-linear 

i objective function 
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Problem case 2: V 1=4 v2=3 v3=2 v4=l 

Satellite Support Scenario with Non-linear 
objective function 
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Problem case 3:vl = l 
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Chapter 8 
Comparison 

this chapter we compare the results obtained by applying previously discussed 
algorithms such as Greedy Search, Tabu Search. Simulated Annealing and the GA by 
Kumar and Bagchi [3], We take here a sample problem with four scenarios, 
(each obtained by changing the weight factor of satellites that are involved in 
visibility clash.) 

8.1 The Results of comparison: 

The problem we consider here is a real life clash situation in which four satellite 
visibilities clash. The visibility start times and visibility end times are as follows. 

vis_slart I = 6520 vis_end 1 = 71 15 

vis slarl 2 7244 vis ciul 2 831S 

vis._start 3 77 1 2 vis end 3 8650 

vis_start 4 == 8027 vis__end 4 = 9 1 96 

I he cases described below are obtained by changing the priority factor of different 
satellites. These (actors determine the total value generated from the supports 
piovided. 1 lie dilferent priority (actor combinations considered are as follows. 


Case 1 : vl =1 

v2=2 

v3=3 

v4-4 

Case2: vl =4 

v2=3 

v3=2 

v4=l 

Case3: vl =1 

v2--4 

v3=3 

v4=2 

Case 4: vl =1 

v2=3 

v3=4 

v4=2 


We considered two different value generating objective functions. The first is linear 

(applied (or Payload operatiotis) and the other is non-linear (applied for TIC type 
satellite operations). 

iMear Objective function: 

Maximise 

1^ ~ 2]((vcr nil i{(e|“Si)*Vi} 


.. ( for PL operations) 



Non-Unear Objective function : 


Maximise 

F = lover all i { A| +( 1 * v,)/.00958} ..(for TTC operations) 


where 

Ai= initial support time * Vj 
Ci = end of support of satellite i 
Si = start of support of satellite i 
Vj ■“ priority weight of satellite i 

Table 8.1 


Best Solutions obtained for the four different algorithms: 


Problem 

Linear objective function 

Non-Linear objective function 

Case 

GA 

Greedy 

Tabu 

Simulated 

GA 

Greed).' 

Tabu 

Simulated 

no. 


Search 

Search 

Annealing 


Search 

Search 

Annealing 

1 

4670 

4670 

4670 

4670 

3478.67 

3415.99 

3481.79 

3481.79 

2 

4467 

4467 

4467 

4467 

3958.48 

3855.55 

3962.04 

3962.04 

3 

4670 

4670 

4670 

4670 

3481.01 

3415.99 

3481779" 

3481.79 

4 

' 4 ^ 4 "^ 

4344 ”’'"'^”' 

4344 ~ 

4344 

2969.49 

2881.89 

2971.344 

2971.344 


I able 8.2 

Computer time ( in milliseconds) taken by four different algorithms: 

The algorithms were run on a Pentium III , 800 MHz processor machine with 128 MB 


RAM and 20 GB HDD. 


Problem 

case no 

Linear objective function 

Non-Linear objective function 

GA 

Greedy 

Search 

Tabu 

Search 

Simulated 

Annealing 

GA 

Greedy 

Search 

Tabu 

Search 

Simulated 

Annealing 

1 

i 

540 

390 

390 

4730 

7460 

490 

490 

4780 

2 

450 

330 

390 

5540 

7523 

550 

600 

5760 

3 

480 

330 

390 

4510 

7765 

440 

500 

4780 

4 

560 

330 

330 

4510 

7612 

440 

600 

4720 



Table 8.3 Support Scenario generated by four different algorithms: 


Problem Case: vl =1, v2=2, v3=3, v4=4 


"Sgorithm 

Support Start Time 

Si 

Support End Time 

Ci 

Total Support 

given: ej - s\ 

GA 

6520 

7000.01 

480.01 

7600.01 

8082.5 

482.489 

7712 

7712 

0 

8482.5 

9196 

513.5 

Greedy Search 

6520 

6520 

0 

7244 

7884 

640 

7712 

11 n 

0 

8484 

9196 

712 

Tabu Search 

6520 

7000 

480 

7600 

8080 

480 

7712 

7712 

0 

8680 

9196 

516 

Simulated Annealing 

6520 

7000 

480 

76()() 

8080 

480 

■ 

7712 

0 

8680 

9196 

516 


8.2 Graphs 

8.2.1 Graph for best solution obtained with different algorithms: 


Best Solution Obtained with different algorithms 
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8.2.2 Graph for Computer time taken by different algorithms: 


Computer time for different algorithms 
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(liraplis for Satellite Support Scenario obtained with different algorithms: 


Satellite Support Scenario with Non-linear 
objective function (GA) 
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Satellite Support Scenario with Non-linear 
objective function (Tabu Search) 
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Satellite Support Scenario with Non-linear 
objective function (Simulated Annealing) 
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8.3 Analysis: 

8.3.1 Best Results Obtained: 

The results obtained .show that there is no dilTcrcncc in the best solution obtained by 
four algorithms when we use the linear objective function for evaluation (i.e. Payload 
visibility su|)port). We lind sonic minor dilTcrcncc in the bc.st value obtained when wc 
use the non-linear objective function, (i.c. for T'l'C visibility support). First we 
consider the Greedy Algorithm. As mentioned earlier, if number of initial solutions is 
increased, it increases the likelihood of reaching to global optima. Here though we are 
using a large number of initial solutions (Chapter 4, Eq. 4.2.5), the shape of the search 
space is very complex and hence the Greedy Algorithm terminates in a local optima 
(which is very near to the global optima). But most of the times we find that this 
algorithm also works fine to get the global optima for this problem. 

The other three algorithms, GA, Tabu Search and Simulated Annealing are well 
known to find a global optimum solution. The critical thing is setting of parameters 
for these algorithms. We, as de.scribed before, optimised the parameters using Design 
ot Experiments. As shown above all these three algorithms are able to find global 
optimum solution. 



8.3.2 Computer Time: 


We tested the algorithms on a IBM machine having Pentium III processor with 
SOOMhz speed, 128 RAM, 20 HDD. The computer time taken by Tabu Search and 
(Irccdy Scarcli is almost a second or less than that for resolving a visibility clash of 
(bur satellites. However, the time taken by Simulated Annealing and GA is 
comparatively more. This can be explained as follows. 

The Simulated Annealing algorithm needs a number of iterations at each temperature 
value and the algorithm proceeds with a small decrement in temperature at each step, 
(slow cooling). 'I'hus, to obtain final result one has to make comparatively a large 
number of iterations in SA. Hence, the computer time required is more as compared to 
Tabu Search or Greedy Search algorithms. 

In case of GA, we have to run the algorithm for as much as 2000 generations to get to 
the best solution obtained by Tabu Search and Simulated Annealing. This is because 
the shape of the search space is very complex. Hence GA also takes comparatively 
more time to reach the global optimum solution. 

The interesting thing to note here is that Tabu Search is able to reach the global 
optima in a comparatively less amount of time because of the diversification criteria 
used. We select different combinations of initial solutions ( as explained earlier) and 
run the algorithm with these solutions as starting solutions, hence the global optimum 
solution is obtained comparatively quickly. 



Chapter 9 
Conclusions 


Tlie problem of sclieduling satellite operations is a complex process due to a sizable 
luiinher and type of constraints imposed from operations, satellites, stations and 
controller, and non-linear objectives and constraints. Also, it is experienced that this 
problem of satellite scheduling is difficult to be mapped to existing solution 
paradigms in the literature, such as linear or integer programming. Thus, conventional 
methods cannot be used to produce a satellite operations support schedule. 

One of the critical things in TTC operations scheduling is the designing of objective 
function. Because of many qualitative decisions involved, it is very difficult to map 
all the factors in the objective function. We have developed a exponentially decaying 
function which takes into consideration most of the factors discussed with ISRO and 
which also tackles the issue of dynamic priority. 

In this study we tested four different algorithms on the same problem. Though all the 
three metaheuristic algorithms implemented here, viz, GA, Tabu Search and 
Simulated Annealing and also Greedy Search reached to the optimal or near to 
optimal solution, we could see that Tabu Search reached the global optimum in 
comparatively less amount of time. If we look at the results of the three algorithms 
more closely we find the following things. 

In case of Greedy Search, the algorithm performed very well when we used the linear 
objective function, but many times the algorithm was not able to reach the global 
optimum solution when we used non-lincar objective function. This shows that this 
algorithm is mostly applicable for the clash resolution of the payload type operations 
support where the objective function is linear. 

Genetic Algorithm performed very well on both type of objective functions, i.e linear 
and non-linear. In case of non-linear objective function, GA took more computer time 
to reach to tlic global optimum .solution than other two metalieuristic algorithms ( 1 abu 
Search and Simulated Annealing) tested. 



Tabu Search turned out ,o be the best algorithm among the other two aa it reached 
global optimum solution in case of both the linear and the non-linear objective 

function. Furthermore, it took comparatively much leas time to reach ttte global 
Optimum solution than other two algorithms. 


in case of Simulated Annealing, we find that the algorithm reached the global optima 
m both the cases of hnear and non-linear objective function. The algorithm took 
computer time that is more than Tabu Search but less than GA. 


This work establishes that Tabu Search works more effectively than Genetic 

Algorithm and Simulated Annealing in resolving clashed visibilities for TTC type of 
operations, the key difficulty in satellite scheduling. 



Bibliography 


[1 ] Agnese J C and Brousse P(1 998) “Scheduling Techniques for a 

constellation of Visibilities,” Spaceflight Dynamics; Proc AAS/ 
GSFC International Symposium, Greenbelt, MD, May. 

[2] Bagchi T P ( 1 999), Multiohjeclive scheduling by Genetic Algorithms, 
Kltiwer . 

[3] Bagchi T P, Kumar S (2000), “LEO Satellite Visibility Clash 
Resolution by Genetic Algorithms”, OR-2000 XXXIII Annual 
Convention of the OR Society of India, Ahmedabad. Operations 
ManagcmoiU Conference 4, Dec 22-24, 2000. Chennai. 

[4J Bottchcr A, Jahn A, Lutz E and Werner M (1993),” Analysis of basic 

system parameters of communication networks based on low earth 
orbit .satellites.” 

(5 1 Deb K ( 1 995) Optimisation for Engineering Design: Algorithms and 

Examples New Delhi, Prenticc-llall. 

[6] (ialibcr I' and Dimitrov S A (1998)” A Satellite Mission Planning 

Framework.” 

I?] Glover F, Laguna M (1997) Tabu Search, Klnwer Publication 

[8 1 Cioldberg 1) E ( 1 989), Genetic Algorithms in .search, Optimization and 

Machine Learning. Addition-Wislcy, Reading Mass, USA. 

19| GRI'-AS, (ieneric Resource, F.vcnt and Activity Scheduler. 

http://wvvvv.vcridian.com/Products/grweb2.asp 

[ 1 OJ .lohnston M D and Miller G D, “SPIKE: Intelligent Scheduling of 
Hubble Space 'Felescope Observations” 

http://www,.stsci.edu/apsb/doc/papers-and-meetings/93-lntelligent- 

Sche diiling/spikc/spikc-chapterd.html 

[ 1 1 J Kennedy D J, Parnell G S and Rowell W F (1988) “An Expert System 
for Scheduling Satellite Supports”, Interfaces ,\%\6 Nov-Dec 1988, pp. 
28-34 

[12] Kirkpatrick S, Gelatt C D, Vecchi M P Jr (1983), “Optimization by 
Simulated Annealing”, Science, 13 May 1983, Volume 220, Number 


4598. 



Laarhoven P J M and E.H.L. Arts (1987), Simulated Annealing: 
Theory and Applications. 

Learning Optimisation from Nature, Simulated Annealing, 
ChemoWeb, Chemometrics Research Group, Naval Research 
Laboratory, Washington DC. 

http://chemdiv-www.nrl.navv.mi1/6 1 1 Q/sensors/chemnmetrics/optsa.html 
Mostert S and Milne G W (2000), “Sunsat: Solutions for Remote 
Sensing”. 

http://sunsat.ee. sun. ac.7a/sspapers/sats.95/sat.s.html 
Mukherjee A(1998), “Skirmishes in the sky”, Computers Today, 
March 1998. 

Rao J D, Soma P, Padmashree G S (1998), “Multisatellite scheduling 
system for Leo Satellite Operations.” 
h(ip://www.nasda.go.ip/SpaccOps/paper98/track2/2b002.pdf 
Thiel S IJ, “Tabu Search, Overview of Modem Heuristic Methods.” 
hitp://egLM-ia.cs.cf.ac.iik/lJser/S.U.Thiel/ra/subsection3_7_2.html 
Wolfe W J and Soresen S E (2000), “Three scheduling Algorithms 
Applied to the Earth Observing Systems Domain,” Management 
Science, Vol. 46, No 1, January 2000 pp. 148-168. 



APPENDIX 1 
COMMON PROGRAMS 


1 .Evaluation . java 

import java.io.*; 
import java.util.*; 
import java. lang. *; 
class Variable 
( 

int [ ] S, E, s, e, s w, s n; 


Var iabl e ( ) 

{ } 

Variable ( int [ ) vis^start , int [ ] vis^end, int [ ] allocate_start, int [ ] 
allocate end, int [) sat_wt,int[] sat_no) 

{ 

S = n€^w int [vis_start . length] ; 
s = new int [allocate^start . length] ; 

E - new int [vis_end. length] ; 
e new ;i nt [allocate_end. length] ; 
s w - new int [sat_wt . length] ; 

£3 n new int [sat_no. length] ; 

i 0 r ( 1 r X t .t ' 0 ; 1, < S * 1 e n g t h ; i + 1 ) 

( 

S [ i ] -vis_start [i] ; 
s [ i ] -allocate_start [i] ; 

E [i] =vis_end[i] ; 
e [i]=allocate_end[i] / 
s _w [ i ] = s at_wt [ i ] ; 
s_n [ i ] " s a t_no [ i ] ; 


class Eva 1 nation 


Eva I uat i on ( ) 

{ ) 

public static double evaluation {Variable vb) 
throws Exception 
( 

Variable vrb - new VariableO; 

vrb-vb; 

double result; 

int []w; 

double [)A; 

int min_support = 480; 

w " - r \ p w i n t [ v r 1) . s _w . 1 e n g t h ) ; 

A ^ new d o u b 1 e [ v r b . s _ w , 1 e n g t h ] ; 



for (int i=0; i<vrb. s_w. length; i++) 

( 

w|i| ( V( b. s _ w I 1 I ) ; 

A[i|- (min_support*w[i] ) ; 

} 

result^O; 
int flag; 

for (int i=0; i<vrb. s_w. length; i++) 

{ 

if ( (vrb.e[i] -vrb.s[i] )>=inin support) 

flag = 1; 
else 

flag = 0; 

result = result+(A[i] + { ( (l-(Math.pbw(2. 718281, (- 

( (vrb. e [ i] -vrb. s [ i] - 

min support) * ( .00958) ) ) ) ) ) *w [ i ])/. 00958 ) )*flaq; 

) 

return rosu It ; 

} 


2. Feasibility . java 

class Feasibility 
{ 

Feasibility { ) 

( ) 

public static boolean feasibility (Variable vr) 

{ 

double min = 480, max = 983, recon =600; 

Variable vrb = new Variable (); 
vrb =vr; 

boolean flag = false; 

for (int i=0; i<vrb. S . length; i++) 

( 

f lag=false; 

if ( (vrb.s [i] <=vrb.e[i] ) && (vrb.e[i]<=vrb.E[i] ) 

& & ( vrb . S [ i ] <=vrb . s [ i ] ) ) 

{ 

if ( ( ( (vrb. e [i] -vrb.s [i] ) >=min) | | ( (vrb.e [i] - 
vrb. s [i] ) ==0) ) && ( (vrb.e [i] -vrb. s [i] ) <=max) ) 
{ 

if ( (vrb.e [i] -vrb.s [i] ) >=min) 

( 

if (i ! = (vrb . S . length-1) ) 

{ 

if ( (vrb . s [i+1 ] -vrb. e [i] ) >=recon) 
flag=true; 

else 

brealc; , 

} 

O 1 

flag-true; 



} 

el se 
I 

i f ( 1 ! { vrb . S . length-1) ) 

{ 

if ( (vrb.s[i+l]-vrb.e[i] )>= 0 ) 

f lag=true; 

else 

break; 

} 

else 

f lag=true; 

} 


else 

break; 

} 

else 

break; 


return fla 


g; 


3 , FineTuning . j ava 

import j ava . 1 ang . ; 
import java. io. * ; 
import j ava . u t i 1 . * ; 

class FineTuning 
( 

FineTuning ( ) 

{} 

public static Variable f ineTuning (Variable__array v^array) 
throws Exception 
{ 

Vector v; 

Vector vl; 

Vector v2; 

Variable va « new VariableO; 

Variable_array varray = new Variable_array ( ) ; 
double optimum=0; 
varray = v_array; 

f o r; ( i n t m= 0 ; m< va r r a y . v r ba . 1 e ng t h ; m+ + ) 

{ 

double 

f inalmax-0, previous=0, next=0, index__previous=0, index__n 

ext-0; 

boolean f lag=true, flagl, flag2=true; 
double answer-0, answer 1; 
double ( ) t:emp; 

int 1 ter=10000,prob,countl=0; 



Random rm new Random ( (long) 19. 00) ; 

Variable vrb= new Variable(); 

Variable vrb_final = new Variable! ); 

Feasibility fb = new Feasibility!)/ 

Evaluation ev « new Evaluation!); 

Neighborhood! nh = new Neighborhood! () ; 

Variable_array vr_array = new Variable__array(); 

vrb=varray .vrba[m] ; 
for(int j=0;j<iter;j++) 

{ 

vr_array = nh. neighborhood! (vrb) ; 

V = new Vector! ) ; 
vl « new Vector!); 

f lag fb. tea55ibility! vrb); 
if !hag) 

! 

answer=ev. evaluation (vrb) ; 
if (answer>f inalmax) 

{ 

f i n a 1 ma x=ans we r ; 
vrb_f inal=vrb; 

} 

previous^next ; 
next finalmax; 
if !previous==n.ext) 

{ 

if !f]ag2) 

( 

f 1 ag2--f alse; 
index_previous = j; 
index_next=j ; 

} 

COuntll+; 

index_next=j ; 

} 

if ! (countl==20 ) && ( (index_next- 
index_previous ) ==19) ) 
break; 

n 1 ‘If* 

( 

if (countl=-20 ) 

f lag2=true; 
countl=0; 

} 

for (int i=0; i<vr_array. vrba. Length; i++) 

{ 

flagl = fb.feasibility(vr_array.vrba[i] ) ; 
if(flagl) 

answer], . oval uat i on (vr array . vrba { i} ); 



answGr) ) ! 


v.addElement(new Double(answerl- 
vl.addE lament {new Integer (i) ) ; 

} 

Enumeration vEnum = v. elements () ; 
double max=-12000; 

int element=0, index= 0 ,indexl= 0 ; 

v2 = new Vector ( ) ; 

temp = new double [v. size ()] ; 
while (vEnum. hasMoreElements ( ) ) 

{ 

temp [element] = 

Double. parseDouble (vEnum. nextElement ( ) .toString 

0 ) ; 

11(1 1 1 ' 1 ('nK'nt 1 > max) 

( 

if (temp [element] >max) 

{ 

max=temp [element] ; 
index = element; 

) 

else 

{ 

max=temp [element 1 ; 

v2 . nddE.l omont (new Integer (element ) ) ; 

} 

} 

element++; 

} 

i f (vl . isEmpty ( ) ==false) 

{ 

if ( v2 . isEmpty ( ) ==true) 

{ 

vrb = 

vr array . vrba [ Integer .parseint ( vl . elementAt (ind 
ex) . toString { ) ) ] ; 

) 

else 

{ 

float interm = ( (rm.nextFloat { ) ) * (v2 . size ( ) +1) ) ; 

prob= (int) interm; 

if (prob==v2.size() ) 

indexl=index; 

else 

indexl=^Integer .parseint (v2 . elementAt (prob) .toSt 

ringO); 

vrb = 

vr array . vrba [Integer. pa rselnt(vl. elementAt ( ind 
exl) .toStringO ) ] ; 

} 

} 

else 

break; 

} 

also 



break; 


} 

i f ( flag) 

{ 

i f ( f i nalmax>~optimum) 

( 

op t i mum= f i na Iruax ; 
va=vrb_f inal ; 

) 

} 

else { } 

} 

return (va) ; 

) 

} 


4 . initialSolution . java 

import java.lo.^'; 
import java. util . 
import java, lanq, * ; 

class I n i t i a 1 S o u t i o n 

( 

s t: cl t i c:: i n t:. [ ] v s s t a r 1 0 ; 

s t a t i. c i n t [ ] vis ^ o n ci 0 ; 

static ini [) saL_wtO; 

static int [] sat_noO; 

static i nt" mi n 4 90; 

Sta lie i f 1 1 t ' - r'en ('0 0 ; 

steit i e int, max OR ^ ; ‘ 

static ini [|[] binary; 

publ 1 c Sint i c Vfrr i abl c array i nit ial solution ( InputClass i class ) 
I h r'f>wn PKCf 'pi I f )!! 

{ 

Variable [] vra; 

Variable array varray; 

vis^^ startO - new int [( iclass . S . length) ] ; 

vis_endO = new int ( (iclass. E. length) ] ; 
sat_wtO =“ new int [ (Iclass. s_w. length) ] ; 
sat noO new int [ (iclass, s_n. length) ] ; 

for ( int i~0; Kiel ass , S , length; it+) 

{ 

vis_startO [i] ' == iclass.S[i]; 
vis_endO(i] « iclass. E[i]; 
sat_wtO[i] « iclass .s_w[i] ; 
sat'*no0[i] » iclass . s_n[i) ; 

} 

int. n - ■" V i n 1: a r t; 0 . 1 e n q t: h ; 
int b 

double p ™ Math . pow (2 . 0, n) ; 



vra = new Variable [ (int) p] ; 

int [][ lbinary = new int [ (int) p] [n] ; 

for ( int i~=0; i< ( int) p; i + +) 

{ 

b=c ; 

for (int j=0; j<n; j++) 
i 

binary [i] [ j]=b%2; 
b=b/2; 

} 

c + ; 

} 

for (int i=0; i<p; i++) 

{ 

int count=0; 
for(int j=0; j<n; j++) 

{ 

if (binary [i] [j]==l) 

( 

count++ ; 

) 


} 

int []vls_start = new int[count]; 
int []vis_end = new int [count]; 
int []all__start = new int [count]; 
int []all__end = new int [count]; 
int [ ] sat_wt = new int [count]; 
int ( ] sat__no new int [count] ; 

i, n t s 0 ; 

for (int j="0;j<n; j + t) 

( 

if (binary[i] [ j]==l) 

{ 

V I nlartl.s) vi n nliartO [ j ] ; 
vis end(sl "■ vis_GndO[j); 
sat_wt[s] = sat_wtO[j]; 
sat no[s] = sat_noO[j]; 


int )c=count-l; 

for (int j=n-l;j>=0;j — ) 


If (l<--0) 

{ 

if (binary [i] [ j]==l) 

[ 

all__start [ )<] =vis_start [)c] ; 
al 1 end [)<] =a Upstart [)c) +min; 

k__: 

} 



} 

el se 
{ 

if (k=» (count-1) ) 

{ 

if (binary[i) [ j]==l) 

{ 

all_end[k]=vis_end[k] ; 
all_sta,rt[k]=all end[kl -min; 
k— ; 

} 

} 

else 

{ 

if (binary [i] [ j 1==1) 

( 

a 1 1 _end [ k ] =Ma th . min ( ( al l_s t a r t [ k+ 1 ] - 

recon) , vis_end[k] ) ; 

all start [k] =all_end[k] -min; 

k--"7 

} 

) 

} 

vra[i] " new 

Variable (vis_start , vis^end, all_start , all_end, sat_wt, sat_no) ; 

} 


varray “-new Vari able^array ( vra) ; 
return (varray) ; 

} 

} 


5 . InitialSolutionl , java 

import java.io.*; 
import java. vit:il . *■ ; 
import java * lanq . S' 

//this class is useful for greedy heuristic., 
//this genierates more than 2"" n initial solutions, 
//the formula is given below 

//formula - (sigma i=3 to n) (nCi*2'" (i-2) -nCi) 

class InitialSolutionl 
{ 

I n i t i a 1 S o 1 u 1 1. on X ( ) 

( } 

static int [] vis^startO; 
s t: a 1 1 c 1 n t [ ] v i s_e nd 0 ; 
static int [] sat_wtO; 
static int [] sat^noO; 
s t a t: i,. c i n t m 1 n ^ 4 8 0 ; 
static int recon«600; 



static irit max=983; 


public static Variable . 

throws Exception (inputClass 


Variable [] vra; 

Vari able^array varray; 
double p7d=0; 
int n,b-0,c-0; 


vi s star to 
vis_endO 
sat_wtO = 
sat noO = 


new '.nt:( (Iclass.S. length) ] ; 
new int[ (iclass.E. length)]; 
new int[ (iclass.s_w. length) ] • 
new int [ (iclass . s_n. length)]; 


for ( int i=0; Kiclass.S. length; i++) 

[i ] = iclass.S[i]- 
vis_endO[i] = iclass.Efi]; 
•'’nt_wtO[i] = iclass.s wfi]- 
sat noO[.i] =■= lclass.s>[i]; 


n vis .startO. length; 
p ' Math.pow( 2 . 0 ,n) ; 

i f (n<3) 
d’-'p; 

elseZ/formula = 2''n+(sigma i=3 to n) (nCi*2- (i-2)-nCi) 

i f (n==3) 
d 0; 

i f {n--“4 ) 
d=23; 
if (n==5) 
d=64; 
if (n==6) 
d=186; 

1 f ( n -- 7 ) 
d=551; 
if(n==8) 
d-lG'lS; 

} 

int [][] binary = new int [ (int)p] [n] ; 

int [][] binaryl = new int [ (int)d] [n] ; 
for (int i=0; i< (int)p; i++) 
b--c ; 

for (int i'^O; j<n;i++) 

( 

binary [i] [j]=b% 2 ; 
b=b/ 2 ; 

} 

Cff; 

) 


iclass) 



int K-0; 

double count!; 

for (int i=0; i< (int)p;i++) 

{ 

count 1 =0 ; 

f o r ( i n t j =0 ; j <n ; j 4-4* ) 

{ 

if (binary[i] [ j]=-l) 
count 14- 4-; 

} 

if (count 1<3) 

{ 

f or ( i nt J 0; j<n; j M ) 
binary! [x] [ j ] ==binary[i] [ j] ; 

X4-4' ; 

} 

e .1 s e 

{ 

for (int k==0; k<Math,pow(2, (counti-2 ) ) / k+H- ) 

for(int j==0; j<n; j4-4') 
binary! [x] [j ]^"bjnary(i) [j]; 

X If ; 

} 

} 

} 

b o o 1 e a n f ! a g “ t r u e ; 

int i_4=0; 

int i_5=0; 

int i_6-0; 

int i_7-0; 

int i_8=0; 

vra == new Variable [ (int) d] ; 

int nev;_index=0; 
for (int i=0; Kd; i + +) 

( 

int count=0; 

for (int j=0; j<n; j4-4*) 

( 

if (binary! [i] (j]=”!) 

COUnt + 4 "; 

) 

int []vis_start = new int [count]; 
int []vis_end » new int [count] ; 
int []a!]_start = new int [count]; 
int []al!_end = new ’int [count] ; 
int []sat_wt * new int [count]; 
int []sat_no » new int [count]; 

In! n- 0; 

for (int j=0; j<n; j4-+) 
f 


if (hinaryl [i] [ j]--*!) 



s start [s] = vis_startO( j ] ; 
vis ond[s] = vis_endO[j]; 

;jat_wt[s] '= sat_wtO[j]; 
sat jio[s] = sat__noO[j] ; 
n + )- ; 

1 

} 

switch (count ) 

( 

case 0: break; 
case 1 : 

all_start[ 0 ]=vis_start[ 0 ] ; 
all_end[0]= all_start [0] +min; 
break; 

t ' ti \ 

al l st art ( 0) =vis_start [0] ; 
all end [01= all_start [0] +min; 
a 1 1 f'lid [ 1 I - vi n nntl [ 1 ) ; 
a 1 1 St. art 1 1 ] “al l_end ( 1 ] -min; 
t vrna k; 
c.' a s e 3 * 

i f ( t 1 aq ) 

( 

an__start [ 0] =vis_start [0] ; 

a 1 1 end [ 0 ] ■” a.ll__start [0] +min; 

a.l L_end [2] =vis_end[2] ; 

all__start [2] =all_end[2] -min; 

all_start [l]=Math.max( {all_end[0] +recon) , 

vis__.start[l] ) ; 

all __end [ 1 ] = all_start [1] +min; 

f lag=false; 

break; 

) 

else 

{ 

all_start [ 0] =vis_start [0] ; 
all_end[0]= all_start [0] +min; 
al] _end [2 ] =vis_end [2] ; 
all_start [2] =all_end[2]-min; 

all__end(l]-Math.min( (all_start [2] -recon) ,vis_end[l] ) 
all__start [1] = all_end[l] -min; 
f lag=true; 
break; 

} 

case 4 : 

if (i__4==0) 

{ 

all__start [0] =vis_start [0] ; 

all__end[0]= all_start [ 0] +min; 

all__end[3]=vis_end[3] ; 

all start[3]"all end[3]-min; 

all^tartfl] « Math.max( (all_end[0] +recon) , 

v.i s__start [ 1] ) ; 

all cnd[l]= all__star,t(l]+min; 

al 1 .st art [2] =’ Math.max ( (all_end[ 1 ] +recon) , • 

vis .start(2)); 



all_end[ 2 ]= all_start [2] +min; 

i_4++; “ 

break; 

) 

if (i_4==l) 

{ 

all__start[ 0 ]=vis_start[ 0 ] ; 
all_end[0]= all_start [0] +min; 
all_end[3] =vis_end[3] ; 
all_start [ 3] =all_end[3] -min; 

^ ^ ^ start [ 1 ] —Math , max {(all end [ 0 ] + re con ) 

vi s_ start [ 1 ]) ; ~ 

all__end[l] = all_start [1] +min; 

a 1 1 2 1 —Math . rain ( ( al l_st art [ 3 ] -recon ) , vis end [ 2 ] ) ; 

all start(2]= al 1__ end [2 ] -mi n; ~ 

i I I ; 

break; 

} 

if(i 4““2) 

{ 

ail__start [0] =vis_start [0] ; 
all ond(01"" all start ( 0) +min; 
a I 1 1] • vl.'! oiulUl ; 

all__start [ 3] =all_end[3] -min; 
all__end(l]=vis_end[l] ; 
all__start [1]= all_end(l]-min; 
all__start [2]=Math.max( (all_end[l] +recon) , 
vis__start [2] ) ; 

all__end[2]= all_start [2] +min; 

i_ /I 4. . 

break; 

) 

i ( ( i ^ A) 

{ 

all_start [0] «vis_start [0] ; 
all_end[0]*' all_start[0]+min; 
all_end[3] =vis_end[3] ; 
a 1 1 _ s t a r t [ 3 ] 1 l_end [ 3 1 -min; 

all_end[2] -Math. min ( (all_start [3] -recon) , vis___end[2] ) ; 
all_start [2] - all_end(2] -min; 

all_end[l]®=Math.min( (all_start [2] -recon) ,vis_end[l] ) ; 

all start [1]“ all_end[l] -min; 

i_4=0; 

break; 

case 5: 

if(i_5««0) 

{ 

all_start [0] =vis_start [0] ; 
all_end[0]“ all_start [0] +min; 
all_e^d[4 ] =vis_end[4] ; 
all^start [ 4 ] *«all_end[4 ] -min; 
all‘^start [1] «Math»max( (all_end[0] +recQn) , 
vis start [ 1] ) ; 

all^end[l]- all_^start[l]+min; 

al r’st.art [2]«Math.max( (all^endf 1] +recon) , 

vis .start [2]); 



al 1 endl2]- all__st;art (2) +min; 
a 1 1 _start [ 3 ] =Math . max ( (al l_end [ 2 ] +recon) , 

vis_start [ 3] ) / 

all _end[3]= all_start [3] +min; 

i 5 f -f ; 
break; 

} 

if (i_5==l) 

{ 

all_start [0] =vis_start [0] ; 
all_end[0]= all_start [0] +min; 
all_encl[4] =vis_end[4] ; 
all_start [4] =all_end[4] -min; 
all_start [1] ^Math.max ( (all_end[0] +recon) , 

V i a nt a rt f 1 ) ) ; 

a 1 1 ( 1 ] a 1 1 nl ,ir t 1 1 | i mi ri; 

all_start [2] =Math.max ( (all_end[l] +recon) , 

vls_start 1 2] ) / 

a 1 1 end [ 2 ] a 1 l_s t a r t ( 2 ] +min ; 

al 1 _end ( 3 ] . min ( ( all_end [ 4 ] -recon) , vis_end [ 3 ] ) ; 

all start: [31 - all_end[3] -min; 

i Sft; 

break; 

} 

if(i_5'^^2) 

{ ^ 

all_start [ 0] =vis_start [0] ; 

all end[0]- all_start [0] 4-min; 

a 1 1 ond [ 4 ] v Is^end ( 4 ] ; 

al l_start [ 4 ] -ail_end[ 4 ] -min; 

all_start [ 1 ] »Math .max ( (all_end [0] +recon) / 

vis_start ( 11 ) ; 

all_end [ 1 ] « all_start [1] +min; 
all_end[2] =vis_end[2] ; 
all_start [2] =all_end[2] -min; 
all_start [3] =Math .max ( (all_end[ 2] + recon) / 
vis_start [ 3] ) ; 

all_end[3]= all_start[3]-fmin; 

i 5 + *f ; 
break; 

} 

if(i 5«=3) 


all_start [0] -vis_start [0] ; 

al 1 end [0] » all_start [0] +min; 

all_end [4 ] “vis_end [4 ] ; 

all_start [4 ] *«all_end[4] -min; 

all start [ 1 ] *Math.max ( (all_end[ 0] +recon) r 


vis_start [ 1] )* ; 

ali_end(l]« all_start [l]+min; . 

all end[3]=Math.min( (all_start[4]-recon) ,vis_enai jj 


all 

all 


^start[3]®=‘ all_^end [ 3] -min, , vis endr21)* 

end[2]=Math.min{{all_start[3]-recon) , vis_enai/j ) , 


all~start[2]= all_end[2] -min; 


i 5 4’ + ; 


break; 




if(i 5==4) 

all__start [ 0] =vis_start [0] ; 

a]I_end[0]= all_start[0]+min; 

all__end[4] =vis_end[4] ; 

all _start[4]=all_end[4]-inin; 

al ] __end [1] =vis_end[l] ; 

all_start [ 1 ] =vis_end [ 1 ] -min; 

all_start[2]=Math.max( {all_end[l]+recon) , 

vis_ start [ 2 ] ) ; 

al l_end (2 1 = all_start [2] +min; 

all_start[3]=Math.max( (all_end[ 2 ]+recon) , 
vis__start [3] ) ; 

all_end(3]= all_start [3] +min; 

i 5 + + ; 

break; 

1 

lf(l 5-=5) 

{ 

a 1 1 start [ 0] =vis_start [0] ; 

all r'nfi[01= all atari [0] +min; 

a 1 1 ^aiid [ 4 ] “vis__end [4 ] ; 

al 1 start [ 4 ] =all_end[ 4] -min; 

al 1 end[ 1 ] =vis_end[l] ; 

all start [ 1 ] =vis__end[ 1] -min; 

al 1 start [ 2 ] =Math .max ( (all_end[l]+recon) , 

vi.s_start [2 ] ) ; 

all (•■nd(21 all_ start [ 2 ] -tmin; 

a L 1 ond ( 3 ] - Math .min ( (all_start [4] -recon) , vis_end[ 3] ) 
all start [3] = all_end [ 3] -min; 

l !) i t ; 
tiraak; 

} 

if (i_5— 6) 

{ 

all_start [0] =vis_start [0] ; 
all_end[0]= all_start[0]+min; 
all_end[4] =vis_end[4] ; 
all^start [4]=all__end[4]-min; 
all start [ 3 J =vis_start [3] ; 
all__end[3]= all_start(3)+min; 

all^end[2] =Math.min ( (all_start [3] -recon) , vis_end[2] ) 
all__start [2] = all_end[2]-min; 

all~end[l]=Math.min( (all_start [2] -recon) , vis_endt 1] ) 
all_start[l]= all_end[l] -min; 
i_ 5++ ; 
break; 

} 

if(i 5==7) 

{ 

all_start[0]=vis_start[0] ; 
all_end[0]= all_start [0] +min; 
all_endC4]=vis_end[4] ; 
all start [ 4 ] =all_end(4] -min; 

al l'’'end [ 3 ] =Math . min ( ( al l_start 1 4 ] -recon ) , vis_end [ 3 ] ) 
al restart [3] = all end[3]-min; 

all an, 1(21 Mai h.inlriK.-ill .nt .'U t ( V] - r nrion) , v I r. '-nri ( 2 1 ) 



all__start [2] = all_end[2] -min; 

al I end [ 1 ] =Math . min ( (all_start [2] -recon) , vis end [ 1] ) 

a.l l start. [1]- all_end(l] -min; 

i_5=0; 

brealc; 

} 

case 6: 

if ■. 

{ 

all_start [0] =vis_start [0] ; 
all_end[0]= all_start [0] +min; 
all_end[5] =vis_end[5] ; 
all_start [5] =all_end[5] -min; 
all_start[l]=Math.max( (all_end[ 0 ]+recon) , 

vis___start [ 1] ) ; 

all__end [ 1 ] = all_start [ 1 ] +min; 

al 1 start [2] =Math.max{ {all_end[l]+recon) , 

vl s start 121); 

ai l_end [2 ] = all_start [2] +min; 

a.l l__start [3] =Matti.max( (all_end[2] +recon) , 

vis___start [3] ) ; 

an__end[31" all_start [ 3] +min; 
all__start [4] *Math,max{ (all_end[3) +recon) , 
vis_start [ 4 ] ) ; 

all_end[4]= all_start[4]+min; 

i 6 + + ; 

brealc; 

) 

if(i 6==1) 

{ 

all start [0] =vis_start [0] ; 
a I 1 f'tui f 0 1 ■ all .start [0] +mi n ; 

() 1 1 i !) 1 vl a_ tnul [ !)] ; 
ail_start [ 5] =all_end[5] -min; 
all__start [1] =Math.max( (all_end[0]+recon) , 
vis_start [ 1] ) ; 

a 1 1 _end [ 1] = all_start [1] +min; 
all__start [2] =Math.max( (all_end [ 1] +recon) , 
vis_start [2] ) ; 

all_end(2]= all_start [2] +min; 
all__start [3] =Math.max ( (all_end [2] +recon) , 
vis _start [3] ) ; 

all__end[3]= all_start [3] +min; 

all~end [4 ] =Math7min ( (all_start [5] -recon) , vis_end[4] ) 

all _ start [4] = all___end[4] -min; 

i_67+; 

brealc; 

} 

if (i_6==2) 

{ 

all_start [0] =vis__start[0] ; 
all_end[0]= all_start [0] +min; 
all_end[5]=vis_end[5] ; 
all__start {5] =all_endt53-min; 
all_start[l]=Math.max( (all_end[0]+recon) , 

vis startlH); 

all"end [ 1 ] = all__start [1] +itiin; 



st32rt [2] Msth . ru3x ( ( 31 1 6 nci [ 1 ] +x 6 con ) 

vi s start [2] ) ; ~ 

n.l ,1 _nnd [ 2 ] -r all_start [ 2 ] +min; 

a^ .l_f5)nd[3] =vis_end[3) ; 

a U pstart [3] = all_end[3]-min; 

all_^start [4] =Math.max( (all_enit 3 ] +recon) , 

vis_start [ 4 ] ) ; 

all_end[4]= all_start[4]+iiiin; 

i_6++; 

break; 

} 

if (i_6==3) 

{ 

all__start [ 0] =vis_start [0] ; 

al l__Gnd [0] * all_start [0] +min; 

all_ en(l[5] =vis_end[5] ; 

all start [ 5] =all_end[5] -min; 

a I l_st:art [ 1] =Math.max ( (all_end[0] +recon) , 

vis__start [ 1 ] ) ; 

all__end[l] = all_start[l]+min; 
all_start [2] =Math.max( (all_end[l] +recon) , 
vis_start [2] ) ; 

ari__end[2] = all_start[2]+min; 

al l__end[4] =Math.min( (all_start [5] -recon) ,vis_end[4] ) 
all_start [ 4 ] = all_end[4]-min; 

al,l_end[3]»=Math.min( (all_start [4 ] -recon) ,vis_end[3] ) 
all _start [3] = all_end[3] -min; 
i 6 1 ' i ; 
break; 

} 

i f ( i 6--4 ) 

( 

a I I .‘)l ,ir I. ( 0 I vl 1. (0) ; 

all_end [0 ] = all_3tart [0] +min; 

all_end[5] =vis_end[5] ; 

all_start [5] =all_end[5] -min; 

all_start [ 1] =Math.max{ (all_end[0] +recon) , 

vis_start [1] ) ; 

all_end[l]= all_start[l]+min; 

all_end [2] =vis_endt2] ; 

all_start [2] = all_end[2] -min; 

all_ start [3] =Math.max( (all_end[2) +recon) , 

vis_start [ 3] ) ; 

all__end[3]= all_start [3] +min; 

all start [4] =Math.max( (all_end[3]+recon ) , 

vis_start 1 4 ] ) ; 

all_end[4]= all_start [4] +min; 

i_6++; 

break; 

) 

if {i_6==5) 

{ 

all start [ 0] =vis_start [0] ; 
all’"nnd[0]= all_start [0] +min; 
all end[5] =vis_end[5] ; 
all^start [5]=all_end[5]-min; 

all start[l]=Math.max( (all_end[0]+recbn), 



vis_start [ 1] ) ; 

all end[l]= all_start [ 1] +min; 

al 1 ond|?] -vis end (2) ; 

all .■!ta!-t(21= all end[21-min; 

al l^ntart: [3] =Math.max ( {all_end[2] +recon) , 

vis start [3] ) ; ’ 

a 1 l__end [ 3 ] = all_start[3]+min; 

[ 4 ] =Math .min ( (all_start [5] -recon) , vis end [4] ) ; 
all_start[4]= all_end[4] -min; 

i 6't+ ; 

break; 

} 

1f(i 6 -f.) 

{ 

all_start [ 0] =vis_start [0] ; 
all__end[0]= all_start [0] +min; 
all_end[5] =vis_end[5] ; 
all_start [5] =all_end[5]-miji; 
all_start [1] =Math.max( (all_end[0] +recon) , 
vis__start [ 1] ) ; 

all__end[l] - all_start [ 1) +min; 

all_ end [2 ] =vis_end[2] ; 

a 1 1 start [2] = all_end[2] -min; 

a 1 1 _end [3] =vis_end[ 3] ; 

all start [3] = all_end[3] -min; 

a 1 l_start [4] =Math.max ( (aXl_end[3] +recon) , 

vis_start [4] ) ; 

all_end [ 4 ] = all_start [4 ] +min; 

i__6 + + ; 

break; 

} 

if {i_6-=7) 

{ 

all__start [0] =vis_start (0) ; 

al ,1 end [0] = all_start [0) +min; 

all end ( 5 1 -vi s_ end [ 5) ; 

all_start [5] =ail_end(5] -min; 

all_start [1] =Math.max { (all_end[0] +recon) , 

vis_start [1] ) ; 

all_end[l]= all_start [1] +min; 

all_end[4]=Math.min( (all_start [5] -recon) ,vis_end[4] ) ; 
all^start [ 4 ] = all_end[4] -min; 

al r"end( 3] '=Math .min ( (all_start (4 ] -recon) , vis_end[31 ) ; 
ail_start [3] = all_end[3] -min; 

all2end[2]=Math.min( (all_start [3] -recon) ,vis_end[2] ) ; 

all^start [2] = all_end[2] -min; 

i_6++; 

break; 

} 

if(i_6==8) 

{ 

all .start [0] =vis_start [0] ; 
a 1 1 "end ( 0 ] - a 1 1 _start [ 0 ] +min ; 
a 1 1 end [ 5 ) - v I s_ end ( 5 ] ; 
all__start [ 5] =ail_end[5] -min; 
all end[l]=vis_end[l] ; 
all”start [ 1 ] =all_end[l]-min; 



f ^ ^ Msth . max { ( al 1 end[l] +r6con) 

vis_start [ 2] ) ; ”* ’ ' 

a 1 1 _end [ 2 ] = a 1 l^start [ 2 ] +min ; 

^ ^ tart, [ 3 ] -Math, max ((all e nd [ 2 ] + ra con ) , 

vis St art [3]); “ 

all [ 3] - all^start [3] +min; 

ci 1 l^s t a r t [ 4 ] —Math . max ( (all and [ 3 ] ■f-rscon ) , 
vis_start [ 4] ) / 

all_end[4]= all_start [4] +min; 

i_6++ ; 

break; 

} 

if(i 6—9) 

( 

all_start [0] =vis_start [0] ; 
all_end[0]« all_start [0] +min; 
all_end[5] «vis_end[5] ; 
al l_start [5] -all_end[5] -min; 
all_end[l] =vis_end[l] ; 
all_start [ 1] =all_end[l] -min; 
all_start [2]=Math.max( (all_end[l]+recon) , 
vis_start [ 2] ) ; 

a 1 1 end [21- all^ start [2] +min; 
a 1 1 si art [ 3] " Math .max ( (all_end[2] + recon) , 
vis_start ( 3] ) ; 

al l_end [ 3 ] = all_start [3] +min; 

al l_end [ 4 ] «Math . min ( ( all_start [ 5 ] -recon) , 

vis_end[4 ] ) / 

all_start [ 4 ] - all_end[4] -min; 

i_6ii ; 
break; 

} 

if {i_6=»10) 

{ 

all_start [0]=vis_start [0] ; 

all_end [ 0 ] = a Upstart [0] +min; 

all_end [5] -vis_end[5] ; 

all^start [5] =ail_end[5]-min; 

all_end[l] “vis_end[l] ; 

a 1 l_s t a rt [ 1 ] =al l_end [ 1 ] -min; 

a Upstart [2] =Math,max ( (all_end[l] Brecon) , 

vis_start [2] ) ; ^ 

all_end ( 2 ] - a Upstart [2] +min; 

a 1 l_end [ 3 ] =vi s_end ( 3 ] ; 

all_start [3] « all_end[3] -min; 

all_start [ 4 ] =Math.max ( (all_end[3] +recon) r 

vis_start ( 4 ] ) ; 

all_end[4]- all_start [4] +min; 

1^6++; 

break; 

) 

( 

rj 1 1 St art ( 0) vis si nr l [0 1 ; 
all jnu:i[Q]- all^start [ 0] +min; 
all^end [5] =^vis_end'[5] ; . 
all start [5] =all___end[5]-min; 



all__end[l] =vis_end[l] ; 

al 1 start [1 1 =all_ end [ 1 ] -min; 

all_star:l [ 2) ==Mat'h.max ( (all_end[l] +recon) , 

vis___start [2] ) ; 

all ond[21 == aU_start [2] +min; 

dl 1 _end [ 4 ] -Math . min ( (all_start [5] -recon ) , vis end [ 4 ] ) ; 
all_start[4]= all_end[ 4] -min; 

all_snd [ 3 ] =Math . min ( (all_start [ 4 ] -recon) , vis end [ 3] ) ; 
all_start [3] = all_end[3]-min; 
i 6 ■» + ; 
break; 

} 

if (i_6==l2) 

{ 

all__start [0] =vis_start [0] ; 
all__end[0]= all_start[0]+min; 
all_end[5]=vis_endt5) ; 
all__start [5] =all_end [5] -min; 
all_end[l] =vis_end[l] ; 
all_start [ 1] =all_end[ 1] -min; 
all__end[2] =vis_end(2] ; 
al l__start [2] = all_end[2] -min; 
all_start [3] =Math.max( (all_end[2] +recon) , 
vis__start [3] ) ; 

all_end[3]= all_start[3]+min; 
all_start[4]=Math.max( (all_end[3] +recon) , 
vis_ start [ 4 ] ) ; 

all end[4]- all start [ 4 ] +min; 
i_6 in- 
break; 

) 

i I ( i (. I '. ) 

{ 

all_start [ 0] =vis_start [0] ; 

all_end[0]= all_start [0] +min; 

all_end[5] =vis_end[5] ; 

all_start [5] =all_end[5]-min; 

all_end(l]=vis_end[l]; 

all_start(l]=all_end[l]-min; 

all_end[2] =vis_end[2] ; 

all__start [2] = all_end[2]-min; 

all'^start [3] =Math.max( (all_end[2) trecon) , 

vis_start [3] ) ; 

all_end[3]= all_start[3]+min; 

all”end[4 ) =^Math.min,{ (all_start [5] -recon) , vis_end[4] ) ; 

all_start [4] = all_end[4]-min; 

i_6++; ^ 

break; 

} 

if(i 6==14) 

{ 

all__start[0]»vis_start[0] ; 
all__end[0]= all_start [0] +min; 
all end[5] =vis_end[5] ; 
all_^start [5]=aTl_end[5] -min; 
al l_cnd [ 1 ] ='vis_end[l] ; 
all start [l]=all_end[l]-rain; 



all__end[2] =vis_end[ 2 ] ; 
all^start [2]= all__end[2] -min; 
a 1 1 __nnd ( 3 ] =vi s_end [ 3 ] ; 
all__start [3]= all_end[3]-min; 

all__start[4]=Math.max( (all_end[3]+recon) , 

vis__start [ 4 ] ) ; 

all_end[4]= all_start[4]+min; 

i_6++; 

break; 

} 

if (i_6==15) 

( 

all__start [0] =vis_start [0] ; 
all__end[0]= all_start[01+min; 
a] l_end[5]=vis_end[5] ; 
a 1 ,1 j', t, a r t [ 5 ] =al l_end [ 5 ] -min ; 

al 1_ ond[4 ] =Math.rain( (all_start [5] -recon) , vis_end[4] ) ; 
a.l ,l__start [4]* all_end[4]-min; 

all._end[3]==Math.min( (all_start[4]-recon) ,vis_endt3] ) ; 

aUj5tart[3]= all_end[3] -min; 

nil end [2] *Math.min( (all_start [3] -recon) , vis end[2] ) ; 
nl l^stnrt [2] = all_end[2] -min; 

al l __ond[l ]=Math.min( (all_start [2] -recon) ,vis_end[l] ) ; 
nil, start [1]= all end[l]-min; 

break; 

} 

('ii.'ii' 7: 

11(1/0) 

( 

all Kt art [ 0] -vis_start [0] ; > 

all oMui[OI^ all_start [0] +min; 
al l.__end [ 6] =vis_end[ 6] ; 
all__start [6] =all_end[6]-min; 
all__start [1] =Math.max( (all_end[0] +recon) , 
vis_start [ 1] ) ; 

a 1 1__ end [ 1 ] = a 1 1 __st a rt [ 1 ] +min ; 

nil fst.art [2] ■■■Math.max( (all_end[ll irecon) , 

vis start [2] ) ; 

nil _ end[2]= all__start [2] +min; 

al 1 _start [3] =Math.max ( (all_end[2] +recon) , 

vis__start [ 3] ) ; 

all_end[3]=' all_start [3]+min; 
all~start [4] “Math.max ( (all_end[3] +recon) , 
vis__start [ 4 ] ) ; 

all_end(4]= all_start[4]+min; 

al 1 _start [5] »=Math.max( (all_er.d[4] +recon) , 

vis _start [ 5] ) ; 

al 1 __end { 5] = all_start [5] +min; 

i / 1 t ; 

break; 

1 

if(i 7~=!l) 

{ "■ 

all start [0] =vis_start[0] ; 
al l'~end[0]’^ all_start[0]+min; 
all end [ 6] =vis_end [ 6] ; 



all__start [6]=all_end[61-min; 

all start [ 1] =Math.r[iax ( (all end[0] +recon) 
vis slartd! ) ; 

all_end[l]= all_start [ 1] H-min; 
n 1 1 _ s t a r t [ 2 ] =Math . max ( (all end ( 1 ) +recon ) , 
vi s^stnrt [2] ) ; 

all_end[2]= all_start[2]+min; 

all_start[3]=Math.max( (all_end[2]+recon), 

vis_start [3] ) ; 

all_end[3]= all_start [3] +min; 
all_start [4}=Math.max( (all_end[3]+recon) , 

vis_start [4] ) ; 

a 1 1 end [ 4 ] =• al l_st art [4]+min; 

'^^■^■_®f''^l5]=Math.min( (all_start [6] -recon) ,vis end(5] ) ; 
all start[5]= all_end[5] -min; 
i _7 ++ ; 
break; 

} 

if (i_7-=2) 

{ 

an_ start [01 "Vis_start[0] ; 

al l^ondlOJ-^ all_start(0] +min; 

nl l__end[6]=vis_end[6] ; 

nl 1_ start [6]=aTl_end[6]-min; 

a.ll_start [l]=Math.max( (all_end[0]+recon) , 

vi s^start [ 1 ] ) ; 

al l_cnd(l]= all_start [1] +min; 
all__start [ 2 ] “Math. max ( (all_end[l] +recon) , 
vis__.start [2] ) ; 

al l_end [ 2 ] = all_sta,rt [2] +min; 

al 1_ start [3]*Math.max( (all_end[21+recon) , , 

vi.s''.sfart[31); 

all nnd(3]- all_start ( 3] +min; '' 

all__end[4 ] '=vis_end[4] ; 

a 1 1 start [ 4 ] =al l_end[ 4 ] -min; 

all start [')] Ma t h . max ( (a 1 1 end [ 4 ] i recon ) , 

V i .s_.‘;tar t. [ 5J ) ; 

all end[5]~ all .start ( 5] +min; 

17 11; 
break; 

} 

if(i 7=-3) 

{ 

al l__stnrt [0] =vis_start [0] ; 
all"end[0]= all_start [0] +min; * 

al l__end[6]=vis_end[6] ; 
all__Start [6] =all_end[6] -min; 
al l3start[l]=Math.max( (all_end[0J+recon) , 
vis__start [ 1 ] ) ; 

all endll]= all_start[l]+min; 
all'‘start(2]=Math.max( (all_end{l]+recon) , 
vi s__start [ 2] ) ; 

all end[2]= all_start [2]+min; 

al r start [3] *Math .max ( (all_end(2] +recon) , 

vis .st.art[3]); 

al, l2®nd[3]= all_start[.3J+min; 

ail end [ 5 ] “Math .min ( (all_start [ 6] -recon) , vis_end [ 5] ) 



all__start[5]= all_end[5] -min; 
all_end [ 4 ] =Math .min ( (all_start[5]- 
nl.l _start[41= all end[4]-min; 

:1 1 t ( ; 

bfonk; 


recon) ,vis_end[4] ) ; 


} 

if(i 7-=4) 

( 


all_start [0] =vis_start[ 0 ] ; 

all_end[0]= all_start [0] +min; 

all_end[6] =vis_end[6] ; 

all__start [ 6] =all_end[6] -min; 

nil nt art; [1 ] -Math .max ( (all^ end [0] +recon) , 

vi s _stat t [ 1 ) ) ; 

a.ll_end[l]= all_start [1] +min; ^ 

all, .start [2] =Math.max( (all_end[l]+recon) , 
vi.s ;'.tart[2]); 

al i__cnd[2] = all_start [2] +min; 
all_end[3] =vis_end[3) ; 
al l_.start [3] '=all_end[3] -min; 


al 1 .start [4] =Math. max ( (all_end[3] +recon) , 

V i ,>!l art ( 4 ] ) ; 

■ ill end [4]-" all_start [4] +min; 
a 1 1 .start. [ 5] - Math.max{ (all_end[4] +recon) , 
vis .start [ .51 ) ; 

all f'nd[51= all start [5] +min; 

1 J I t; 
break; 

) 

,if(i 7==5) 

{ 

al 1 _start [ 0] *vis_start [0] ; 

all_end[0]= all_start [0] +rain; 

aJ l__end[6] =vis_end[6] ; 

all_start [6] =all_end[ 6] -min; 

all__ start [ 1] =Math .max ( {all_end[0] +recon) , 

vi.s__start ( 1] ) ; 

all_end[l]= all_start [1] +min; 
all__start [2] =Math.max{ {all_end[l] +recon) , 
vis__start [2] ) ; 

all__endt2]= all_start [2]+min; 

■al 1 _ end[3] =vis_end[3] ; 
al l__start ( 3] =all_end[3] -min; 
al l^start [ 4 ] =Math.max ( (all_end[3] +recon) , 
vis__start [4] ) ; 

■all endl4]= all start[4]+min; 

all (-'iid 1 ') 1 Math .min ( (al l_start [6] -recon) , vi3_end[5i ) 
< 1 1 1, _ start ( 5| = ail_end(5] -min; 
i_7 + + ; 
break; 

1 

if{J 7 — 6) 

{ 

A 1 1 fitr'irt [0] =vis_start (0] ; 
air’end[0]- all_start [0] +min; 



all and 1 fi] ■ vi s_end(6] ; 
all ntrirt 16) af]__end[6]-min; 

■‘11 ni ar t 1 1 1 Math. max( (all end[0] +recon) , 

vin 111); 

all an.ldl al, l,_start[l]+inin; 

all ar t 12) Math.max((all end[ll+recon) , 

via at ai I (71 ) ; 

all an. 117) ali_.-5tart[2,]+min; 

all and ( 3 ] ■ vis__end [3] ; 

a] 1 [ 31 .^aTl_end[3] -min; 

a 1 1 _ end [ 4 ] =vi s_end [ 4 ] ; 

an_ start [4]=:aTl_end[4]-min; 

all .mart [51=^Math.max( (all end[4 ] +recon) , 

V i ; nt a rt [ h ] ) ; 

all i'n<1[!)l* all start [5] +min; 

1 /It; 
bri'ak; 

) 

i ( ( i 7 / ) 1 

( 

al l .-datt 1(11- vi.s^start[0); 

all nn.l(oi all .start [0] +min; 

a 1 1 iMii i t (t 1 vis end[ 6] ; 

al 1 M aril (11 all jtnd [ 6 1 -min; 

all .-.t.iit 11 ) Math. max( (all end[0]+recon) , 

vl.n .-dart [ID; 

all nndlD a 1 1 !d art. ( 1 ) imin; 

all .•datli;i| Math.max( (all end[l] +recon) , 

vi.‘i :d art (71 ) ; 

.ill an(il2)'^ all .start[2]+Inin; 

all and IS] Mal.h.min( (all_start[6]-recon), 

V i and 1 ') 1 ) ; 

all :;tari(‘il-' al l__end [ 5] -min; 

all aii.l [ 4 1 ''Math .min ( (all_start [5] -recon) , 

vi.'-. nnd[4]); 

all .•:tart[4]* al l_end ( 4 ] -min; 

all and [ 3 1 ''•Math .min ( (all_start [4] -recon) , 

V i end ( 1 j ) ; 

alt ;;t .art (.31'-- all__end[3) -min; 
i / 1 t ; 
la nak; 

) 

11(17 (1 ) 

( 

all rd .11 t 1 0) •■■vis_start (0) ; 
all ntid(0)"= nll_start[0)+min; 
n 1 1 end [ 6 ) -•’vi .s_end [ 6] ; 
all ?;t .art ( 6] ■»all_end( 6] -min; 

,il I '.st art t) I’^Math-maxl (all_end[0]+recon), 

V i ;; .'!( .art [ D ) ! 

.ill tMidlD* a ll__start [ 1] +min; 

.1 1 1 etuil 7 1 "v Intend [2] ; 

.ill .'d art t2]»all_end(2]-min; 

,ill .‘dart 1 3) '*Math.max( (all_end[2]+recon) , 
V i a ;d ar t I 3) ) ; 

a I 1 and I lb' all .'start [3] +min; 

all .'daitli) Math .max ( (all end(31 +recon) , 



vi s^start [4]); 

all_end[4 1= all_start [ 4 ] +inin; ' 

all__start[5]=Math.max( (all_end[4] +recon) , 

vis start [5]); 

a 1 1 onci ( 5 ] - a 1 1 __st a rt [ 5 ] +mi n ; 
i / 1 I ; 

brsak; 

} 

if (i_^7==9) 

{ 

all start [0] “vis_start [0] ; 

nl 1 on(i[0] • all__start[0) imin; 

a 1 1 _end [ 6 ] i s_end [ 6 ] ; 

all start: [ 6) «all_end[6] -min; 

a 1 I s t a r t [ 1 ] -Math . max ( ( al l_end [ 0 ] +recon) , 

vis St a r t 1 1 ) ) ; 

all fMuil, 1 all_start tl] tmin; 
a 1 L imd [ 2 ] -vi s_end [ 2 ] ; 

.il 1 St art; [2]«all_end[2]-min; 

ril 1 start [3]=Math.max{ (all_end[2] +recon) , 

v i s start [ 3] ) ; 

all tsK,l(31"" all^ start [3] +min; 

al I .St .,n t (r)“^"^Math.max( (all_end[3] +recon) , 

vi s St at; t ( 4 ) ) ; 

a 1 I raid ( ^ ) " all star t [ 4 ] +min ; 

all o 1 1 d I ') I Ma I 1 1 . m i n ( ( a 1 1 s t a r t ( f> ) “ raron } , 

V i .‘5 (Hid ( 1 1 ) ; 

a ] 1 s t. a r t [ 5 ] - a 1 l_end [ 5 ] -min ; 
i 7 f f ; 
lu (sik ; 

} 

ii:(i 7 
{ 

all _start [0] =vis_start [0] ; 

all_end[0] = all_start [0] +min; 

a 1 l_end [ 6 ] -vi s_end [6]; 

a Upstart [6] =all_end[6] -min; 

al 1 _start [ 1] ^Math.max ( (all_end[0] Brecon) , 

V i s_star;t [1] ) ; 

al 1 end [ 1 ) « all_start [1] +min; 

al 1 ond[2]«vis_end[2] ; 

all .start: [ 2 ] **all_end[2] -min; 

al 1 St art [3] -Math.max{ (all_end[2] +recon) , 

vj s^sStart [3] ) ; 

al l^^^end[3] all_start [3] +min; 
a 1 l_end [ 4 ] wvis_end[4 ] ; 
all^start [4] »all_end[4]-min; 
all^start [5] ==Math.max( (all_end[4] trecon) , 

vis_start [5] ) ; 

a 1 1 end [ 5 ] - al l_start [ 5 ] +min ; 

i 7 } f ; 

hroak; 

} 

lf(i 7^^11) 
f 

al I St ,ir t(0) vis si art 10) ; 



all end [01= all_start [ 0] +inin; 

a 1 1 end ( 6] =vis_end['6] ; 

a 1 1 n t art: [ 6 ) =ai l__end [ 6 ] -min ; 

all !5t art. [l]=Math.max( (all end[0] +recon) , 

via nt.art [1] ) ; 

nll_(!nd[l]= all_start [1] +min; 
a 1 1 _end [2 1 =vis_end[2] ; 
nJ l_.sf;art[2]=all_end[2]-min; 
a Upstart [ 3] =Math.max ( (all_end[2] +recon) , 
vls_ start [ 3] ) ; ' 

all _end[3] = all_start [3] +min; 

a],l end[51=Math.min( (all start[6] -recon) , vis end[5]); 
al I :d atl |'i| - ,il 1 ond ( -mi li; 

end[4]=Math.min( (all_start [5] -recon) , vis_end[4 ] ) ; 

al 1_ start [ 4 ] = all_end[4] -min; 
i 7 f I ; 
lit oak; 

) 

if (i' 7==]2) 

( 

a 1 1 St art [ 0] --vl s^start [0] ; 

all end(()]= all_start (0] +min; 

al 1 _ond [ 6] -vis_end[6] ; 

a 1 1 start [ 6] =ail_end[ 6] -min; 

al l^ntart. [ 1 1 =Math.max ( (all_end[0] +recon) , 

V i s_ start [ 1] ) ; 

all c'nd(l]= all start [ 1 ] +min; 

,1 1 I Olid I 1 V I M (nKl[21 ; 
all ;it,a rl, 1 2 | a 1 1 end [ 2 ] -min; 
a ,1 1, ond [ 3 1 ■ vi and ( 3 ] ; 
al! start [ 3]'ai I_ ond(31 -min; 

.ill .'il.irt[4) Matii. max ( (all_erid[ 3] +recon) , 
vis jd.art [4] ) ; 

al 1 _end[4]= all_start [4 ] +min; 

al 1 _start [5]=Math.max( (all_end[4] +recon) , 

vi .s^start [5]); 

all_end[5]= all_start [5] +min; 

i 7 1 1 ; 

broak; 

} 

if(i 7==13) 

( 

al l_start [0]=vis_start [0] ; 

all_end[0]* all_start[0]+min; 

all__end ( 6] “vis_end [6] ; 

all_start [6]*all_end[6]-min; 

al l__start ( 1 ] =Math .max ( (all_end [ 0 ] +recon) , 

vis .5tart[l]); 

all^endd]” all_start(l]+min; 

all end[2]=vis_end[2] ; 

a 1 1 2.S t a r t ( 2 ] =all_end [ 2 ] -min ; 

nil end[3) =vis_end[3] ; 

alf'start [3]=all_endl3]-min; 

al restart [4] “Math. max( (all_end[3] +recon) , 

vl s__start [ 4 ] ) ; 

nil end[4]= all start [4] +min; 

a 1 1 ‘"end [ 5 ] =Math7min ( (all_start [ 6] -recon) , vis_end [ 5] ) 



all_start[5]= all_end[5] -min; 

i 1 + + ; 
break; 

) 

11 ( i 7^—14 ) 

{ 

all_start [0] =vis_start [0] ; 

all_end[0]= all_start [0] +min; 

all_end[6]=vis_end[6] ; 

a] 1 start[6]=all_ end[6]-min; 

al 1 start [ 1 ] =Math.max ( (all_end[0] 4 recon) , 

vis_start [1] ) ; 

al]__end[l]= all_start [ 1] +min; 
a 1 1 _ _end [ 2 ] =vi s_end [ 2 ] ; 

all _start [ 2] =all_end[ 2 ] -min; 

a 1. 1 _nnd [ 3 ] =vi s_end [ 3 ] ; 

^ 1 ~ all_end [3] -min; 

al 1_ rnd(4]=vis_end[4] ; 
al 1 _start [4] = all_end[4]-min; 
al.l__start [51=MathTmax( (all_end[4]+recon) , 
vis_start [5] ) ; 

all end[5]= all start [5] +min; 

1 7 t't ; 
break; 

) 

if(i 7--15) 

{ 

all__start [0] =vis_start[0] ; 

all_end[0]= all_start[0] +min; 

all_end [6] =«vis_end [ 6] ; 

all__start [ 6] -all_end[6] -min; 

all_start [ 1] =Math.max ( (all_end[0] +recon) , 

vis_start [ 1] ) ; 

all_end [1 ] =" all_start [1] +min; 

all_end[5] =Math.min ( (all_start [6] -recon) , vis_end[5] ) ; 
all_start [ 5] = all_end[5]-min; 

all_end[4] =Math.min ( (all_start [5] -recon) , vis_end[4] ) ; 
a Upstart [4] = all_end[4] -min; 

al l__end[ 3] --Math. min ( (all_start [4] -recon) , vis_end[3] ) ; 
all start [3] = all_end[3] -min; 

a 1 J end [ 2 ] -Math .min { (all_start [3] -recon) , vis_end[2] ) ; 
a 1 I ;; t a r; t [ 2 1 a 1 l_end [ 2 ] -min ; 

i / M ; 
br'o.ik; 

} 

if(i_7-=16) 

a 11 .St art [ 0 ] -"Vi .'5_start [ 0 ] ; 

all nnd[01-' all_start [0] +min; 

all end [ 6] -vis__end( 6] ; 

all. .start [ 6] -ai l_end[ 6] -min; 

all firKi[l]^vis_end[l] ; 

al 1 .sfarl: ( 1 ! ■•all __end[ll -min; 

al 1 .start f2b-Math. max ( (all_end[l] +recon) , 

vi.s .start [2]); 

all end[21^- all_start [2] +min; 

all' start [3]=Math.max{ (all_end[2] +recon) , 



vi s_start [ 3] ) ; 

airend[3]= all_start [ 3] +min; 


al l_st:art [4]=Math.max( (all end[3] +recon) 
vis__start[4] ) ; ~ , 

a I l__end [ 4 ] = all_start [ 4 ] +min; 
all.__start [5] =Math.max( (all_end[4] +recon) 
vis__start [ 5] ) ; ' 

all_end[5]= all_start [5] +min; 
i__7 + +; 
break; 

} 

if{i__7==17) 

{ 

al 1 start [0] '^vis___start [0] ; 

all eridlO)" all ntart ( 0] itni ti; 

a 1 1 end 1 6 1 vi ;>_ end [ 6] ; 

a 1 1 start ( 6) =ail_end( 6] -min; 

all end] 1 1 =vis__end[l] ; 

all nl art [ 1 ] =ari_end[l] -min; 

all .start[2]=Math.max( (all_end[l] +recon) , 

vi.-, ,‘;tart[2]); 

all (Midi 2]-= all_start [2] +min; 

.ill .-dart [ 3] -Math.max ( (all_end[2] +recon) , 
vl r. Ml art ( 3) ) ; 

all (Midl'l] al 1 .start [ 3) imin; 
a 1 1 .‘Start ( 4 1 Math. max ( (all_end[3) trecori) , 

V i .‘i is tart 1 4 ] ) ; 

a 1 1 onrt[4 ] -■ all__start [4] +min; 

< 1 1 1 end 1 h) -Math .min ( (all_star.; [ 6] -recon) , vis_end[5] ) ; 
all .‘startlfi]” all_end[5] -rain; 
i _ 7 t r ; 
break; 

} 

if(i 7==18) 

{ 

al l___start [0] =vis_start [0] ; 
all_end[0]= all_start [0] +min; 
all__end[6] “vis_end[6] ; 
al 1 __start [ 6] =all_end[6] -min; 
a 1 l__end [ 1 ] ='vis_end ( 1) ; 
a 1 1__ .start [1] =all_end( 1] -min; 
all__start[2)=Math.max( (all_end[ 1) +recon) , 
vi .s start [2] ) ; 

all ('nd[2]- all_start [2] +min; 

al 1 _.start [ 3] =Math.max ( (all_end[2] +recon) , 

vis start [3]); 

al l''end[3]= all_start [ 3] +min; 
al 1 _cnd[4]=vis_end{4] ; 

,al 1 '.start [4]=all_end[4]-min; 

al 1 "startlSl^Math.maxl (all_end[4]+recon) , 

vis .sta rt 1 .‘j I ) ; 

all fMull5]- all_start[5]+min; 

i 7 I t ; 

bterik; 

1 



if(iJ7==19) 

( 

a 1 i_ start [ 0] =vis_start [0] ; 
al 1 _eiKi[ 0 ]= all_start[0]+min; 

a 1 1 end [ 6 ] =vis_end [ 6 ] ; 
a 11 St a f t [ fi 1 -~ai 1 onrl [ 6 ] -mi ri ; 
a 1 1 (null 1 1 " vis_ond [ 1] ; 

al l^start [ 1] =all_end[ 1 ] -min; 

start [2] —Math . max ( (all and [ 1 ] +r6con ) , 

vis_start [2] ) ; ~ _ ' 

all_end(2]= all_start [2] +min; 
all_start [3] “Math.max( (all end[2] +reconi 
vis__start[3]); 

all_end[3] = all_start [3] +min; 

al l__end[5]=Math.min( (all_start [6] -recon) ,vis end[5) ) 
all St, art, [5]^ all_end[5] -min; 

all ('nci[(|l Math. min((all_start[51-recon), vis end[4]) 
,d 1, rA art [4) = all_endt4] -min; 
i 7 M ; 
ht ('ak; 

} 

lf(i 7-20) 

{ 

all s t a 1 1 1 0 1 "vi s_start [ 0 ] ; 

.ill (nui(01" all__start [0] +min; 

,il 1 mul [ fi 1 -vl .s_end [ 6] ; 

al 1 r.Larl [ 6) =‘-ai l__end[ 6] -min; 

all ('ud ( 1 1 >=vis__end [ 1] ; 

al! fUart [ 1, ] '-■all_end[ 1 ] -min; 

al 1 Mt art ( 2 ] -Math .max ( (all_end[l] +recon) , 

vi n st a rt [ 2 1 ) ; 

fill <nid I 2 I a 1 1 nta r t I 2 1 I ml ii; 

.ill mid I 3 1 - vis__Gnd[3) ; 

a 1 1 .s t a r t [ 3 ] -al l_end [ 3 ] -min ; 

a 1 1 !;t;ar t [ 4 ] “Math .max ( (all_end[3] +recon) , 

V i ; : 1 a r t [ 4 ] ) ; 

,.i 1 1 mid[4]- all_start [ 4 ] +min; 

ri II nt art [ 5] =Math .max ( (all_end [4 ] +recon) , 

vi.n_ f5t .U. L [5J ) ; 

all end [5]= all_start [5] +min; 

.1 _7't t ; 
branik; 

) 

if(.l 7”=21) 

( 

,i 1 1 .start [0] ''Vis_start [0] ; 

.illjnid(0|= all_start[0]+min; 

.ill ('rid ( () ] »vis_end [ 6] ; 

,il 1 .ntart [ 6] ~al l__Gnd( 6] -min; 
all (Mid I 1 1 =vis_Gnd [ 1] ; 

,il 1 [ 1] -^all__end(l]-min; 

a I r St art [2] -Math. max ( (all_end[l] +recon) , 
v i ;; .‘it art (2)); 

,ill muil21= all__start [2)+min; 

,ill mull 31"vls__end[3] ; 

all ;;1 ,irf ( 3] -aTl _ Gnd(3] -min; 

,il 1 st ,it. t [ 4 ) -"Math.max ( (all__end(3] +recon) , 



v:i s_start: [ 4 ] ) ; 


al Ijmci [ 4 ] == all_start [ 4 ] +min; 

al 1 'indl 5| =Math.min( (aU_start[6]- 

.itl iUatt [h]-- a.l 1 end [ 5] -min; 

1 /II; 


recon) ,vis_end[5I ) ; 


lf(i 7 ==22) 

{ 

all_start[ 0 ]=vis_start[ 0 ] ; 
a.ll_ end[0]= al.l_start [0] +min; 
al 1 _end[6]=vis_end[6] ; 
al I nt-nrt [ 61 =al.l_end[ 6) -min; 
all nnd[ 1 l=^vis__end[l] ; 

11 .‘d ar 1 [ 1 1 ~ail_end( 1] -min; 
all .start [2] ==Math.max( (all_end[l] +recon) , 
vis .St at,r [21 ) ; 

ti 1 1 _ (,M\d [2] - al.]_start [2] +min; 

al 1 ond[3] -via__end[3] ; 

al 1, .sinri; [31’^allj'ind[3]-min; 

a 1 1 oiui [ 4 1 "'vi .';__end(4] ; 

all .slatt(41 ai ,l_end[4] -min; 

all ;; f a r t ( 6 1 "•Math . max ( (all_end [ 4 ] +recon) , 

V i .St art (51); 

al! ond[,61« all start [ 5] +min; 

1 1 M ; ~ 

break; 

) 

lid 7»>»23) 
t 

< 1 ! 1 ;•.! art. (0) '"vls_start [0] ; 

<ill endl01= all^start [0] +min; 

all end ( 6 1 »vis__cnd [ 6] ; 

all St art' I 61 -all end[ 61 -min; 

a 1 1 f'lid I 1 I V 1 s end [ 1 ] ; 

all start [1] all_Gnd[ 1 ] -min; 

al 1 St art, [2] -Math.max( (all_end[l] +recor)) , 

V i ;; .s t a 1 1 [ 2 ) ) ; 

all f'nd[2|>^: al l___start [2] +min; 

, 1 1 1 end ( !) ] ’•Math .min ( (all_start [ 6] -recon) , vis_end[5] ) 
all St art 1 ,61 al l_end [ 5] -min; 

all end 1 4 ] •-Math .min ( (all_start [5] -recon) , vis_end[4] ) 
a 1 1 .St art, [ 4 ] » al l_end ( 4 ) -min; 

al 1 enfl [ 3 1 - Math. min ( (all_start [ 4 ] -recon) , vis_end[3] ) 

a 11 .s t a r t 1 3 d all _end [ 3 ] -min ; 

i/M; 

break; 

) 

ltd 7 r-.7A) 

{ 

a ! 1 .St art [ 0 j *vis__start [0] ; 

.ill endlOl'- all__start [0] +min; 

a 1 1 end [61 -v i .<;_ pnd(6] ; 

all s.t at t [ 6 1 -a i 1 end[6] -min; 

all end 1 2 r'vi.s__end[2] ; 

all st.art. [21--all end[2)-min; 



1 1 
al 1 

al ] 

V i ;; 
a 1 1 
a 1 I 

V i r> 
a 

all 
vi s 
all™ 

i y 

hrn 

} 

i f ( 

( 

a ! i 
al 1 
ai 1 
a 1 1 
.1 1 1 
a! 1 
al 1 
a 1 ! 

. 1 1 1 
\/ ! * ! 
al 1 
.tl 1 

V i 
all 
,il I 

V 1 n 
.1 1 1 

i / i 
la f M 


1 I -Math. mint (All oi- 

at art. [ 1 ]** all^end [ 1 J ,vis_end[l] ) 

;;;;;;;|3l-M«h,„ax,,aU_e„dt2Hreco„,, 

'M 1 Math. max( (all endnt + 

an [4] ) ; •"•"-®'^a[3]+recon) , 

_ end [ 4 ] 3a_start [4] +inin; 

,...«f.an[5]=Math.max((all endr4n 

start [5]); ■^aend[4]+recon), 

>dt5]= all__start[5]+min; 

t t , 

ak; 

i 7 2b) 

n ,i( t [O] vi s ntart [ 0 ] ; 

''‘'..stnrtfOJ+min; 

'‘tid((, I vifi endffij; 

•'■’■“"a ■'>' .,nnd(6]-min; 

'‘lull. ’I vi.n ntulf:::]; 
n'.nl l-M all end(2I-mln; 

I ' I Mai li.niin{ (nM start-roi v= 

■''■'■""11 .ill ..,i,iiri-„i„ ' '‘”'°"’'’'»=-“‘‘nii.- 

;;;;;;; I ;|^».i'i...iMMi.iii„;„,|j|,,,„„,^ 

'‘"''Ml ■'! 1 [ Jj 

■> ' I 1 a ^ Mat h . max ( (all__end [3] +recon) , 

'’'"i[4 1 ,11 I Ht artt4J+min; 

;;||;J|;•l^M.dh.min((all 

all end [5] -min; 

1 ; 
k; 


} 

if(i 7 -26) 

I 


a I I 
al 1 
a 1 1 
a 1 1 
.1 1 1 
.1 I I 
.i I 1 
a I 1 
a 1 J 
V } a 
.1 1 1 
al 1 
a 1 j 
a 1 1 
V i ■: 

all 

i / I I 


7 t (()| vi;; ;!lnrt[0]; 

'"* 1 ( 1 ) a 1 1 .Tt.art [0] +min; 

nd| f, I vl •! end [6] ; 

Mad [f,l al 1 end[6J-min; 

•Iid(,.'| Vi;i end (2) ; 

n .u 1 dd -a! 1 end[2]-mln; 

;'"i( 1| Math min((all_start[ 2 ]-recon),vis end[l]) 

M nd ( I J all end [1] -min; 

Mad Ml Mdlh,max( (all end(2] +fecon) , 

M .11 1 f 1) ) ; 

•'I'll a - nil „start(3)+min; 

’lull 4 I vi ;i end (4] ; 

'ad [4 1 -al 1 end[4]-min; 

ladl'd Mat h.max( (all end(4] +recon) , 

I ad I hi ) ; 

"d I 'i ) , 1 1 1 .-.t ,1 1 1 ( h] imin; 



) 

ifd 

I 

. 1 1 I St .1 rl I 0 1 vi s start [0] ; 

a 1 1 I'lu 1 1 0 1 a 1 .1 sta rt [ 0 ] +min; 

all 'Mitilfij viR_end[61; 

al 1 _start [ 6] ==ail_end [ 6] -min; 

.al 1 _and [2] =vis_end[2] ; 
al l..start [2]=all_end[2]-min; 

® ^ 5 “Math . min ( {all_start [2 ] - recon) vis end f Di- 
al l__start [ 1 ] = all_end[l]-min; 
a 1. 1 start [ 3 ] =Math . max ( (all_end [ 2] +recon) , 
vis s tart [3]); 

. 1 1 I^cnd [31= all_start [3] +min; 

all_^end[5]-Math.min( {all_start[6]-recon) ,vis end[5] ); 
all start [51 » all__end[5]-min; 

a! l.,.«'nd[^ l-^Math.min( (all_start[51-recon) ,vis_end[4] ) ; 
al 1 start [41= all_end[4] -min; 
i 7 I 1 ; 
hr 0,1 k; 

1 

il(i 7 =*28) 

[ 

all start |01 vi s^start [0] ; 

all ondlO)- all_start [0] +min; 

all otKi 1 6 1 =vis_end [ 6] ; 

nl 1 st art [ 6) ••ail_end[6] -min; 

all ondl21 vis_ond[2]; 

all St ,irt [ 2 d'al l_ond[2] -min; 

,il I oiidl 1 1 -'Math .min‘( (all_start [2] -recon) , 

V 1 s ond (11); 

.ill si art [ 1 1 “ all _end[l] -min; » 

al! and I 3 1 vis end [ 31 ; 

, 1 1 1 s. I a 1 1 [ 3 ] ' nl l^ find ( 3 ] -min ; 

.ill stall 1 4 ] Math .max ( (allend [3] trecon) , 

V 1 s .st.art [4 1 ) ; 

all end [41= all_start [4] +min; 

al 1 start. [ 5) =''Math.max ( (all__end[4] trecon) , 

vis .start [5]); 

al 1 ond [5]= all._start [5] +min; 
i 7 11; 
hi oak; 

1 

i.f(l 7 =--=2 9) 

( 

a 1 ,1 start: [ 0] =vis_start [0] ; 
al 1 ’’end [01= all_start(0]+min; 

,i 1 1 f'lui [ fi 1 ■'•V 1 .s _ end ( 6) ; 
a 1 1 .St art [ fil =ail_end[6] -min; 

,i 1 1 _ oiul ( 2 1 ■>'vis_end [2] ; 

.ill" start [ 21 =aTl end(2)-min; 

,il 1 '"oiidll )>*Math.min((all_start[2]-recon) ,vis_end[l] ) ; 
,il r St art 1 1 I - al, l_end[l]-min; 

.ill Olid I 31 -'vhe^endlSl ; 

all st.irt 131=ari__end[3]'-min; 

.ill” .St art [4 1 '"•Math. max ( (all_end[3] -trecon) , 

V 1 s st.at t [41); 



} 

if(i 121 ) 

{ 

a 1 1 Ni at I I 0| vi ;; at ari [0] ; 

■>i> 'Mid 10 ] all ;itart[ 01 +min; 

all iMidlf)] vis f'tui(61; 
all ;datt[ftl n.U_end ( 6] -min; 
all and [ 2 1 "vis__ond[2] ; 
all at art- [21 '-al'l__end[2]-min; 

^ I ^ 1 ■ Math .min ( (all_start [2 ] -recon) , vis endf 11 ) • 

all _start[]]= all_endtl] -min; “ 

all s t . n r t [ 3 1 =Ma t h . max ( ( al l_end [ 2 1 + r econ ) , 

V i !' .st art (31); 

al l _(>nd(31i" all__start [3] +min; 

all and['5l*Math.min( (all_start[6]-recon) ,vis end! 51 ) • 
ail .’itartl'll" all__end[5] -min; ~ 

all -nuK^l “Math. mTn((all_start[5]-recon), vis end[ 4 ]); 
all .‘d.att (4 1« aJ,l_end[4] -min; 

lilt; ~ 

ht aak; 

1 

11(1 7 ?.H) 

I 

a 11 .Mt at I ( (1 1 vi .s start [0] ; 

.ill (Midllll al]__.start[0]+min; 

.ill tnidlhl vi;i_end[6]; 

.ill ;:t,itt[fil all end[6)-min; 

. 1 1 1 and (21 V 1 .M ond [21 ; 

.til Mt,irt[71 al 1 _end[2]-min; 

all i'nd[ll MaLh.min'( (all__start[2)-recon) , 

V 1 ;i t'lid I 1 1 ) ; ' 

all fii .It I |1| al 1 end [ 1 1 -min; ^ 

.til .‘nd Ml V 1 F, nrui ( 3 1 ; 
all .St, lit 1 t] .il 1 _ end [ 3] -min; 

.ill .‘vtattl4| M.d ii .max ( {a.1.1 _end [ 31 +recon) , 
u 1 :: ;d, at t [41); 

.ill ('nd[4 1* all _start [4] +min; 

all .St at t. [ 5] -Miath .max ( (all end[4] trecon) , 

vl.M -ttart [h]); 

till (‘nd[h)« all_ start [5] +min; 

1 7 I I ; 
bt aak; 

1 

1((1 7 -=^29) 

{ 

,1 1 1 .Mt ,11 1. [01 «vis__start [0] ; 
all imdlOl” all__start [0] +min; 

,ill <Mid 1 f) M V i s end [ 61 ; 

.ill .start [61 •'■■a.i l_^endt6]-min; 

.ill and [ 2 1 •»vi.'3^end[2] ; 

.ill Ht art, [2 Mall end(2]-min; 

.11 1 and( 1 MM.ath.mint (all_start[2]-recon) , vis_e>^d[l] ) ; 
. 1 1 I .Mt art [ 1 I all end[ll -min; 

.il 1 (Midi 1) ■ vl,s orid[3] ; 

.ill ;;t .iTt [ 3Mai,l_ end(3)-min; 

.ill Mt . 11 1 1 4 1 -^Mat’h.maxt (all_endI3]+recon) , 

V i ; ; ,s t a 1 t [ 4 1 ) ; 



all_end[4]= all_start[4]+min; 

all_end[5]=Math.min({all_start[6]-recon),vis end[5]); 

all_start[5]= all_end[5] -min; 

i 7 + + ; 

break; 

} 

if(i 7==30) 


all^start [ 0] =vis_start [0] ; 
all__end[0]= all^start [0] +min; 
all_end[6] =vis___end[6] ; 
all_start [ 6] =all_end[6] -min; 
all_end[2] =vis_end[2] ; 
all_start [ 2 ]-all__end[ 2 ]-min; 

all^end[l]=Math.min( (all__start[2]-recon) ,vis end[l] ) ; 
all^start [1]= all_end[l] -min; 
all_end[3] =vis_end[3] ; 
al l^^start [ 3 ] = all_end[3] -min; 
all_end[4]=vis_end[4] ; 
all_start [4] = all_end[4]-min; 
a 1 1 _s t a r t [ 5 ] -Ma t h . max ( ( a 1 l_e nd [ 4 ] + r e c on ) , 
vis_start [ 5] ) ; 

all_end[51 = all_start [5] +min; 

i ^^^7 + 4 '; 
break; 

} 

if (i_7=-31) 

{ 

all_start [0]=vis_start [0] ; 
all_end[0]= all_start [0] +min; 
all_end [6] =vis_end[ 6] ; 
all_start [ 6] «all_end[6] -min; 

al 1 _ond [ 5 ] *^Math . min { (all_start [6] -recon) , vis_end [ 5] ) 
all_start [5] - all_end[5] -min; 

all_end [ 4 ] =Math . min ( (all__start [5] -recon), , vis_end [4 ] ) 
all_start [4] = all_end[4] -min; 

all_end [ 3] =Math.min ( (all_start [4] -recon ) , vis_end[3] ) 
all_start [3] « all_end[3] -min; 

a l,l_end [2] -Math. min ( (all_start [3] -recon) , vis_end [2] ) 
all_start [2] = all_end[2] -min; 
all_start [1] =Math.max( (all__end[0] +recon) , 
vis_start [1] ) ; 

all end[l]= all_start[l]-Hmin; 

1 7-0; 
break; 

} 

vra[i] - new 

Variable ( vi s start , vis end, all_start, all_end, sat_wt, sat_no) ; 


} 

varray =new Variable_array (vra) ; 
return (varray) ; 



6 . InputClass . java 

import java.io.*; 

class InputClass 
1 

int[] S, E, s_w, s_n; 

InputClass ( ) 

{} 

InputClass {int[] vis_start, int [ ] vis_end, int [ ] sat_wt,intt] 
sat_no) 

( 

S = new int [vis__start. length]; 

E = new int [vis_end. length] ; 

s w = new int [sat_wt. length] ; 
s n - new int ( sat_no . length] ; 

for ( int i=0; i<S .length; i++) 

( 

S [ 1 ] ^^vis_start [ i ) ; 

F. ( i ] "-vis^end [ i ] ; 
s_w [ i ] - sat_wt [ i ] ; 
s n(i) « sat_no[i]; 

) 

} 

1 


7 .MyTimer • java 

import jciva . util , * ; 

class MyTimer extends GregorianCalendar 
{ 

publ i c long getTimeInMillis ( ) 

{ ' ' 

return super . getTimeInMillis ( ) ; 

public long difference (MyTimer t) 

I 

return this . getTimeInMillis ( ) - t .getTimeInMillis () ; 

} 

} 


8 .Neighborhood. java 

class Variable array 

( 

Variabl e [ ] vrba; 
Var i abl o array ( ) 

{ 



} 

Variable__array (Variable [ ] vrb_a) 

{ 

vrba » new Variable [vrb_a . length] ; 
vrba = vrb_a; 

} 

} 

class Neighborhood 
{ 

Neighborhood { ) 

{ } 

public static Variable_array neighborhood (Variable vb) 

I 

Variable vrb = new Variable (); 

Variable array va; 
vrb ' vb; 

Variable []vrba = new Variable [ (vrb. S . length*4 )] ; 
int step '10; 

for(inL i "O; i<vrb. S . length; i++) 

( 

vrb . s [i ] =vrb. s [i] +step; 
vrba [ i ] -new 

v,n i ab 1 1 ' ( V rb . h , V t b . K, vr b . .s , V rb . e, vrb . n w, vrb. a n) ; 
vrb . s [ i ] =vrb . s [ i ] -Step; 

} 

for (int i=0, j=vrb , S .length; i<vrb. S . length; i++, j++) 

{ 

vrb . s [i] =vrb. s [i] -step; 
vrba [ j ] =new 

Variable (vrb . S, vrb. E, vrb . s, vrb. e, vrb. s_w, vrb. s_n) ; 
vrb . s [i ] =vrb . s [i] +step; 

) 

for (int i=0, j =vrb. S . length*2; i<vrb.S. length; i++, j++) 

vrb . e [ i ] =vrb. e ( i] +step; 
vrba [ j 1 -new 

Variable (vrb . S, vrb . E, vrb . s, vrb. e, vrb. s_w, vrb. s_n) ; 
vrb . e [i ] =vrb . e [i ] -step; 

} 

(,dr ( 1 nt: !■ 0, ) - vrb . B . lonqth *3; 1 <vrb. S . length; i + + , J f f ) 

vrb . e [ i ] =vrb . e [i] -step; 
vrba [ j ] =new 

Variable (vrb . S, vrb . E, vrb. s, vrb.e, vrb. s_w, vrb. s_n) ; 
vrb. e ( i ] -vrb. e [i] +step; 

} 

va = new Variable_array (vrba) ; 
return va; 

} 



appendix 2 

CODE FOR GREEDY SEARCH ALGORITHM 


import java . lang. 
import java.util.^; 
import java.io.^; 

class Greedy 

( 

public static void main (String args[]) 
throws Exception 
{ 

MyTimer start = new MyTimerO; 

GregorianCalendar goal = new GregorianCalendar { ) ; 

System, out . print ( "Time Before: ") ; 

System . out . pri nt (gcal . get (Calendar . HOUR) +" : ») ; 

System , out . print (gcal . get (Calendar . MINUTE) +” : ”) ; 

System, out , pri ntln (gcal . get (Calendar .SECOND) ) ; 

double opt:,] mum ■ 0 ; 

int lonqth of ^ fi ne_tuning =5; 

i rit [ ] vi f-p st, a,rt new int [ 4 ] ; 
j n t: 1 ,1 V :i s (■' n d - \\ e w i, n t: [ 4 ] ; 

i nt [ ] sa t: wt •" new int [ 4 ] ; 
j rit[] sat jio - new int [4]; 

vis^^^start [01 ^:”"^6520; 
vis start f 1 ] ”'7244 ; 
vis"\st art- [2 1 “7712; 
vis_star t [ 3 1 d3 027 ; 

vis end[0J “■7] IS; 
vis^end[l]=8318; 
vis_end[2] ^-8 650; 
vis_end(3]=9196; 

sat__ wt [0] •=! ; 
sat_wt [ 1 ] =4 ; 
sat_wt [2] =^3; 
sat_wt [3] ==^'2; 

sat_no [0] “^1; 
sat_no[l] ^■■2; 

sat no [2 ] "'“3; 

sat rio [ 3 ] -"'4 ; 

V(vd (H v; 

Vector vl; 

Vector v2; 

Variable^array varray « new Variable_array () ; 

Input Cl ass i class « new Input Cl ass (vis_start^ vis^end^ sat__wt , sat_no) ; 
InitialSolutionlModified insol = new InitialSolutionlModified( ) ; 
//TnltialSo.lutionl insol « new Ini tialSolutionl ( ) ; 

FileWriter fw new EMleWriter ( "clash_soiutions" ) ; 

F’i. n e T u n i n q f t “ n e w F i. n e T u n i n g ( ) ; ' 

Fvn 1 u,i t i on < » v 1 rn^w F.va I tia I, i on ( ) ; 



vnrray = i nsol . initialsolutionlmodif ied (iclass) ; 
//varray = insol . initialsolutionl ( iclass) ; 

Variable vrb_optimum= new VariableO; 
Variable_array input; 

Variable []input__a; 
double []input_d; 

V I ' c t o r V 3 • ■ n e w V e c t o r ( ) ; 

Vector v4=new VectorO; 

//varray . vrba . length 

System, out .print In ( "length '="+varray. vrba. length) 

for (int m=l ;m< varray . vrba. length ;m++) 

{ 

double f inalmax=0, previous=0, next= 0 , 
index_previous=0, index_next=0; 
boolean flag=true, flagl, flag2=true; 
double answer=0, answerl; 
double []temp; 

int iter-lOOOO, prob, countl^O; 

Random rm = new Random (( long) 10 . 00) ; 

Variable vrb= new VariableO; 

Variable vrh_.final = new VariableO; 

Feasibility fb = new Feasibility () ; 

Evaluation ev = new Evaluation {); 

Neighborhood nh = new Neighborhood () ; 
Variable_array vr_array = new Variable_array () ; 

vrb=varray . vrba [m] ; 

for (int j=' 0 ; ■) <1 t:or ; j + + ) 

vr_array -- nh . neighborhood (vrb) ; 

v = new Vector ( ) ; 

vl = new VectorO; 

flag = fb.feasibility(vrb); 

if(f.lnq) 

{ 

answer=ev . evaluation (vrb) ; 

//System, out .println { ''answer="+answer) ; 
fw. wr j te ( " "+answer+"\n") ; 
if (answer>finalmax) 

{ 

f inalmax==answer ; 
vrb flnal=vrb; 

} 

previous=next ; 
next = finalmax; 
if (previous==next) 

{ 

if (flag2) 

{ 

fl. aq2 fnl.sc; 
index provinu.'s =- j; 
indiix nnxt,~‘i; 



} 

if ( ( couri 1 1 - 1 0 ) && ( ( index___next-'index previous) == 9 ) ) 

break; ~ 

else 

{ 

if ( count 1=— 10) 

{ 

f 1 nq2'“U: rue ; 

I'euul 1 0/ 

} 

} 

//System, out. println( "Now let’s check the neighborhood solutions") 
for (int i=0; i<vr_array. vrba . length; i++) 

{ 

flagl = fb.feasibility(vr_array.vrba[i]); 

//System, out .println ("Is solution feasible ? (true/false) •-> 
"+flagl); 

if (flagl) 

( 

answer 1-ev. evaluation (vr_array . vrba [i] ) ; 

//System, out .print In ( "answer l="+answerl) ; 

V . addElement (new Double (answerl-answer) ) ; 
vl . addRlomont (new Integer (i) ) ; 

} 

} 

Enumerati on vEnum == v. elements {) ; 

double max=--12000; 

int element»=0, index=0, indexl=0; 

v2 == new Vector ( ) ; 

temp = new double [v. size ()] ; 

while ( vEnum . hasMoreElements ( ) ) 


//System, out . pri nt In ( "vEnum: - "+ vEnum. next Element ( ) ) ; 
temp [element 1 Double . parseDouble (vEnum. nextElement ( ) .toStringO ) ; 
i f (temp [ element] >=max) 


if (temp [element] >max) 

{ 

max=temp [element ] ; 
index = element; 

} 

else 

{ 

max=-tnmp [element] ; 

v2 . addF’.lement (new Integer (element) ) ; 

} 

} 

element++; 

} 

//Sy.s Lem . out . println ( "max = "+max) ; 
if ( vl . isEmpty ( ) ==false) 

{ 

if (v2 . isEmpty {) =;=true) 

vrb = vr_array. vrba [Integer. parseint 
(vl . elementAt (index) . toString { ) ) ] ; 
else 
( 

//System. out . println ( "size : = 


+v2 . size 0 ) ; 



//System. out. printlnC’interm. . "+interm) ; 

prob= ( int ) interm; 

//System. out. println("prob:= " +prob) ; 

if (prob==v2 . size ( ) ) 

indexl=in.dex; 

else 

indexl-IntegGr.parseInt(v2.elementAt(prob) toStrinafn - 
vrb - vr__array.vrba [Integer. parseint ^ ' 

(vl . elementAt (indexl) .toStringO ) ] • 

) 

} 

else 

break; 

} 

el se 
break; 

} 

i f { f 1 aq ) 

i 

for (int. 1 0; i<vrb_final.s. length; i++) 

v3 . aciciMl nment (vrb_ f inal) ; 

v-l . nddl'! 1 nmont: (now Double ( f inalmax) ) ; 

} 

else 
( } 

) 

int'iil ,1 tK'W Vn r i nhl f' I vl ■ n i 7.0 ( ) 1 ; 

Input d now douti I (> [ . .mI /. f- ( ) | ; 

f or ( i nf. i -O; i<input__d. length; i++) 

( 

input d [ 1 ] =DouBle.parseDouble (v4 . elementAt (i) .toStringO ) 

input a[i] " (Variable) v3 .elementAt (i) ; 

} 

double tempi; 

Variable temp2 = new Variable (); 
for ( int i~0; i<input_d. length; i++) 

( 

for (int j=i; j<input_d. length; j++) 

{ 

if (input d [ j ] >input_d [i] ) 

( 

I urtif)! ■ i npiut, d [ j ] ; 
iiifiul- fj( j) - jnput_d[i]; 
i nput fl [ i I t.ompl ; 

I.<‘mp2 -i np>it _a [ j ] ; 
i uf.iul .1 ( i I input__a [i] ; 

,i iif.iut a 1 ,i ] temp2 ; 

} 

} 

) 

Variabiu f]lnput_al = new Variable [length_of_f in 

for, (int, i-d); i < 1 ongth__of_f ine_tuning;i++) 
int.ml, al|il ■■ input__a[i]; 
iufiut now Vari able__array (input_al) ; 
vrt) opt imnm ■" ( t . fineTuning ( input) ; 
optimum evl .evaluation (vrb_optimum) ; 



for ( int opt=0 ; opt<vrb_optimum . s . length; opt++ ) 

System, out . print In (vrb_optimum. s [opt] + ”\t*’+vrb_optimum. e [opt] + ” 
\t ”+vrb_optimum. s__n [opt] ) ; 

System. out .println( "optimum: = "+optimum) ; 

fw. close { ) ; 

GrrMior i anCa 1 ondar qcall - now Groqori anCa 1 ondar ( ) ; 

S ys I can . out . pi i n t ( "Time AI Loi : " ) ; 

System. out .print (gcall.get (Calendar .HOUR) +" : ") ; 

System. out .print (gcall . get (Calendar .MINUTE) +" : ") ; 

System. out .print In (gcall . get (Calendar . SECOND) ) ; 

MyTimer stop = new MyTimerO; 

long minis = stop.difference(start); 

System. out .println ( "Program Took: " + millis + " 
milliseconds" ) ; 

} 


} 



APPENDIX 3 

CODES FOR TABU SEARCH ALGORITHM 


l.Ta±)U.java , 

import java.lang.^; 
irnpo r t j a va . u t i 1 . * ; 
import j a va , i o . * ; 

class Tabu 

( 

public static void main (String args[]) 
throws Exception 
{ 

MyTimer start = new MyTimerO; 

GregorianCalendar gcal == new GregorianCalendar ( ) ; 

System, out .print ( '’Time Before: "); 

System. out .print (gcal ^ get (Calendar . HOUR) +" 

System. out .print ( gcal . get (Calendar .MINUTE) +” : " ) ; 

System. out .println(gcal. get (Calendar .SECOND) ) ; 

do ub 1 e opt imum== 0 / 
int[] vis_start = new int[4]; 
int[] vis^end = new int[4]; 
int[] sat_wt = new int[4]; 
int[] sat__no = new int[4]; 

vis_start [0] =6520; 
vis_start [1] =7244; 
vis_start [2 ] =7712; 
vis_start [3] =8027; 

V i s_e nd [ 0 ] = 7115 ; 
vis_end [ 1 ] =8318 ; 
vis_end[2] =8650; 

V is _end [ 3 ] =9196; 

sat_wt [0]=1; 
sat_wt [1]=2; 
sat_wt[2]=3; 
sat_wt [3] =4 ; 

sat_no [0] =1; 
sat__no [ 1 ] =2; 
sat_no[2]=3; 
sat_no [3] =4 ; 

Vector v; 

Vector vl; 

Vector v2; 

Variable_array varray = new Variable^array ( ) ; 

InputClass iclass = new InputClass ( vis__start , vis_end, sat_wt , sat_ no) ; 

//InitialSolutionModif ied insol = new Init ialSolutionModif ied ( ) ; 

Init ialSolution insol = new Init ialSolution ( ) ; 

FileWriter fw = new FileWriter ( "clash_solutions” ) ; 

FineTuning ft = new FineTuning ( ) ; 

Evaluation evl = new Evaluation () ; 



//varray = insol . initialsolutionmodified(iclass) ; 
varray = insol.initialsolution{iclass); 

//System. out .println ( "I am here”); 

Variable vrb_optimum = new Variable(); 

Variable_array input ; 

Variable [ ] input_a; 

d u b ]. G [ ] i n p u t d ; ^ 

V c t o r V 3 - n g w V (5 c t: o r ( ) ; 

Vector v4=new Vector () ; 

int f inetuning__array_size =3; 
int tabu^tenure = 1; 

//varray . vrba . length 

for (int m=l;m<varray. vrba. length ;m++) 

{ 

double f inalmax=0, previous=0, next=0, index_previous=0, index_next=0 

boolean flag, flagl, flag2=true; 

double answer===0, answer 1/ 

double []temp; 

int iter=150, countl=0,* 

Variable vrb- new Variable (In- 
variable vrb_final = new Variable On- 
Feasibility fb = new Feasibility 0 ; 

Evaluation ev = new Evaluation () ; 

Neighborhood nh == new Neighborhood () ; 

Variable array vr_array - new Variable_array { ) ; 

Variable array vr array feasible ~ new Var i ahle_array ( ) ; 

Select si = new Select (); 

TabuList tl = new TabuList On- 
Aspiration ap = new Aspiration On- 
boolean terminate -false; 

V t b V (^1 r 1 ‘ a y . V r b a [ rn J ; 

flag == fb.feasibility(vrb); 

//System. out . println (flag); 
if (flag) 

{ 

for(int ' j==C)' j'^iter; j++) 

{ 

//System. out . println ( ” iter:= ” + j ) ; 
answer-ev. evaluation (vrb) ; 

/bSys tern, out .println (” answer := ” tanswer) ; 

System. out .println ( ” finalmax:= ” +f inalmax) ; */ 
if (answer>f inalmax) 

{ 

f inalmax=answer; 
vrb_f inal=vrb; 

} 

previous=next ; 
next = finalmax; 
if (previous==next ) 

{ 

if ( f lag2 ) 

{ 

f lag2==f alse; 

1 n c i n X p r e V i o \ i s ™ j ; 
index_next= j ; 

) 

count 1++; 



index_next=j ; 

} 

if ( (countl==20) && ( { index_next-index_previous ) == 19 ) ) 

break; 

else 

{ 

if (count 1==20) 

{ 

flag2=true; 
count 1=0; 

} 

} 

vr_array = nh . neighborhood ( vrb) ; 
vr_array_f easible = si . select (vrb, vr_array) ; 

/ /System . out . print In ( ''length=="+vr_array_f easible . vrba . length) ; 
if ( vr___array__feasible , vrba . length==0) 
break; 
else 
{ 

for (int i=0; i<vr__array_f easible . vrba . length; i++) 

bool nnn fY)nd1t lon « t: 1 . 1 (vrh, vr_array__fna.*^ible . vrba [1 1 ) ; 

//System. out .println ( ”tabu_condition = "+tabu_condition) ; ‘ 
if (tabu_condition) 

{boolean aspiration_condition = 

ap . isAspirant (vr_array_feasible . vrba [i] ) ; 

//System. out .println ( "aspiration^condition = ”+aspiration_condition) 
if (aspiration_condition) 

{ 

String tabu_element= tl . comparision (vrb, vr__array_f easible . vrba [ i ] ) ; 

1 1 . addToTabuList ( tabu_element ) ; 
vrb = vf array_f easible . vrba [ i ] ; 
break; 

} 

} 

else 

{ 

String tabu_element= tl . comparision (vrb, vr__array_f easible . vrba [i] ) ; 
//System. out .println ( "tabu_element = ”-i-tabu_element ) ; 

tl . addToTabuList (tabu_element ) ; 
vrb = vr_array_feasible.vrba[i] ; 
break; 

} 

if ( i== ( vr_array__f easible . vrba . length- 1) ) 
terminate = true; 

} 

if ( j >=tabu_t enure) 

tl . removeFromTabuList ( ) ; 

if (terminate) 

break; 

} 

} 

v3 . addElement ( vrb__f inal ) ; 

v4 . addElement (new Double ( f inalmax) ) ; 

} 

c' .1 s e 

{ } 



} 

System. out .print In ( "size=="+v3 . size { ) ) ; 
input_a = new Variable [v3 . size {)] ; 
input_d = new double [v4 . size ()] ; 

//System, out .println ( "well, I am here") ; 
for (int i=0; i<input__d. length; i++) 

{ 

input_d [i] =Double.parseDouble {v4 .elementAt (i) . toString ( ) ) ; 
input_a [i] = (Variable) v3 . elementAt (i) ; 

} 

double tempi; 

Variable temp2 = new Variable {); 
for (int i=0; i<input__d. length; i++) 

{ 

for (int j=i; j<input_d. length; j+-f ) 

{ 

i f ( input_d [ j ] >input_d [ i ] ) 

{ 

tempi = input_d[j]; 
input_d [ j ] = input_d [ i ] ; 
input_d [ i ] =templ ; 
temp2 “ input_a [ j ] ; 
input_a [ j ] = input_a [ i ] ; 
input_a [i] =temp2; 

} 

} 

} 

Variable []input_al == new Variable [ finetuning_array size]; 

for (int i-0; i<f inetuning__array_size; i + t) 

input_al[i] = input_a[i]; 

input = new Variable_array (input_al) ; 

vrb^optirnum - f t . fineTuning ( input ) ; 

optimum = evl . evaluation (vrb_optimum) ; 

for ( int opt=0 ; opt<vrb___optimum. s . length; opt++ ) 

System, out . println ( vrb_optimum. s [opt] +"\t "+vrb_optimum. e [opt ] +"\t"+vr 
b_opt imum . s_n [ opt ] ) ; 

System. out . println ( "optimum: = "+optimum) ; 

GregorianCalendar gcall = new GregorianCalendar ( ) ; 

System. out .print ( "Time After: "); 

System. out .print (gcall . get (Calendar . HOUR) +" : ") ; 

System. out .print (gcall .get (Calendar .MINUTE) +” : ") ; 

System. out .println (gcall . get (Calendar . SECOND) ) ; 

MyTimer stop = new MyTimerO; 

long minis == stop.difference(start); 

System. out .println ( "Program Took: " + millis + " milliseconds"); 

} 

} 



2 . Aspiration . java 


class Aspiration 

{ 

static double previous_best; 

Aspiration 0 

{ 

previous_best = 0; 

) 

public static boolean isAspirant (Variable vr) 
throws Exception 
{ 

Variable vrb = new Variable () ; 
vrb=vr ; 

Evaluation ev = new Evaluation () ; 

double this_solution = ev. evaluation (vrb) ; 

/*System. out . println ( "this_sol : ="+this_solution) ; 
System. out .println ( "pre_best : ="+previous best ) ; */ 


if (this_solution>previous_best) 

{ 

previous_best=this_solution; 
return (true) ; 

} 

else 

return (false) ; 

) 


} 


3 .Select. java 

import java .util . 
import j ava . io . * ; 

class Select 

{ 

Select 0 

{ } 

public static Variable__array select (Variable vr , Variable__array vrr) 
throws Exception 
{ 

Variable_array varray = new Variable_array ( ) ; 

Variable_array varray__return; 

Variable vrb = new Variable (); 

Feasibility fb = new Feasibility () ; 
vrb = vr; 
varray == vrr; 

Variable [ ] vrb__temp; 
double [ ] f easible__array; 

Evaluation ev = . new Evaluation () ; 
doubl e answer ev . evaluation (vrb) ; 



Vector V = new Vector (); 

Vector vl = new Vector (); 

for (int i=0; i<varray. vrba . length; i++) 

{ 

boolean flagl = fb. feasibility (varray. vrba [i] ) ; 
if (flagl ) 

{ 

double answerl=ev. evaluation {varray. vrba{i] ) ; 
v.addElement (new Double (answerl-answer) ) ; 
vl , addElement (varray . vrba [i] ) ; 

} 

) 

Enumeration vEnum = v. elements () ; 

Enumeration vlEnum = vl . elements () ; 

f easible_array = new double [v. size ()] ; 
vrb__temp = new Variable [vl . size ()] ; 

int element = 0; 

while ( vEnum. hasMoreElements ( ) ) 

{ 

feasible_array [element] = 

Double . parseDouble (vEnum. next Element { ) . toString ( ) ) ; 
vrb_temp [element] = (Variable) vlEnum. nextElement { ) ; 
element++; 

} 

int u=0; 
double tempi; 

Variable temp2 =new Variable (); 

for (int i=0; i<feasible__array. length-1; i++) 

{ 

for (int j=i+l; j<feasible_array. length; j++) 

{ 

if ( feasible_array [ j ] >feasible_array [i] ) 

{ 

tempi * feasible_array [ j ] ; 

feasible_array [ j ] = f easible_array [i] ; 

feasible_array [i] =templ; 

temp2 = vrb_t emp [ j ] ; 

vrb_t emp [ j ] ” vrb_t emp [ i ] ; 

vrb__temp [i] =temp2; 

} 

} 

} 

varray _return - new Variable_ar ray ( vrb_temp) ; 
return (varray_return) ; » 

} 


4 .TabuList . java 

import java.io.’^; 
import j a va . lang . * ; 
import java. util.*; 



class TabuList 

( 

TabuList ( ) 

{} 

static Vector v = new Vector (); 

public static boolean isTabu (Variable cur_sol, Variable new sol) 

( 

Variable current_solution = new Variable!); 

Variable new^solution = new Variable (); 
boolean ret; 

current___solution = cur__sol; 
new^solution = new_sol; 

String index = comparision {current_solution, new_solution) ; 
if (v. contains (index) ==true) 

{ ret = true; } 
else 

{ ret = false; } 
return(ret); 

} 

public static String comparision (Variable cur_sol, Variable new_sol) 

{ 

Variable current_solution = new VariableO; 

Variable new_solution = new VariableO; 
current_solution = cur_sol; 
new^solution = new__sol; 

String sr =" ”; 

for (int i=0; i<current_solution. e . length;i++) 

{ 

if (current_solution . s [i] ! =new__solution. s [i] ) 

{ 

s r~ ” s ” + i ; 
break; 

} 

if (current_solution. e [i] ! =new_solution . e [ i] ) 

{ 

sr=”e''+i ; 
break; 

} 

} 

return (sr) ; 

} 

public static void addToTabuList (String s) 

{ 

String str = s; 

V. addElement (str) ; 

) 

public static void removeFromTabuList ( ) 

{ 

String removed = (String) v. remove (0); 

//System. out . println ( "yes removed" ) ; 

} 

} 



APPENDIX 4 

CODE FOR SIMULATED ANNEALING 


1. Simulated Annealing, java 

import j ava . lang . ^ ; 
import j ava . ut il . * ; 
import java.io.*; 

class SA 

{ 

public static void main (String args[]) 
throws Exception 
{ 

MyTimer start = new MyTimerO; 

GregorianCalendar goal = new GregorianCalendar ( ) ; 

System. out .print ( "Time Before: "); 

System. out .print (gcal.get (Calendar . HOUR) +" : ") ; 

System. out .print (gcal . get (Calendar .MINUTE) +" : ”) ; 

System. out .print In (gcal.get (Calendar .SECOND) ) ; 

int[] vis_start = new int[4]; 
int[] vis_end = new int[4]; 
int[] sat_wt = new int[4]; 
int[] sat_no - new int[4]; 

vis_start [0] -6520; 
vis___start [1]=7244; 
vis_start [2]=7712; 
vis__start [3] =8027 ; 

vis_end[0]=7115; 

vis_end[l]=8318; 

vis_end[2]=8650; ' 

vis_end[3]=9196; 

sat_wt [0] =1; 
sat__wt [1]=3; 
sat_wt[2]=4; 
sat_wt[3]=2; 

sat_no [0] =1; 
sat__no [ 1] =2; 
sat_no [ 2 ] =3 ; 

.s a t n o [ 3 ] - 4 ; 

Variable_array varray = new Variable_array ( ) ; 

InputClass iclass = new InputClass ( vis_start, vis__end, sat_wt , sat 
//InitialSolutionModif ied insol = new InitialSolutionModif ied ( ) 
InitialSolution insol = new InitialSolution ( ) ; 

FineTuning ft = new FineTuning ( ) ; 

Evaluation evl = new Evaluation () ; 

//varray = insol . initialsolutionmodif ied (iclass) ; 
varray = insol . initialsolution (iclass) ; 

Variable vrb__optimum= new Variable (); 

Variable vrb_final = new Variable (); 

Variable vrb_initial = new Variable (); 

Variable vrb = new Variable (); 

Evaluation ev = new Evaluat ion ( ) ; 



Variable_array input; 

double T, epsilon, r; 

int iter-O, iterl=0/ 

double solution=0, solution_optimum=0; 

int f ine_tune_length 5; 

Variable []input_a; 
double [ ] input___d; 

Vector v3=new Vector (); 

Vector v4=new Vector (); 

FileWriter fw = new FileWriter ( "clash_solutions”) ; 

I t ) r ( in t, nr 0 ; tn< va r r a y . v ! ha . I nriq t. h / rtu f ) 

( 

boolean flag=true; 

T=25; 

epsilon =1; 
r=. 99; 

long seedl = 7; 
long seed2 == 10; 

vrb-varray . vrba [m] ; 
vrb_initial=varray . vrba [m] ; 

Feasibility fb = new Feasibility {) ; 

Neighborhood nh = new Neighborhood () ; 

Variable_array vr__array = new Variable_array(); 
Variable_array vr_array_f easible; 

Select sit « new Select {); 

Random rm = new Random (seedl ) ; 

Random rm_new = new Random (seed2) ; 

while (T>=epsilon) 

{ 

iter=0; 

while (iter<50) 

{ 

iter++; 

vr_array = nh . neighborhood (vrb) ; 
flag = fb.feasibility(vrb); 

//System, out . print In ( ’*Is solution feasible ? (true/false) :-> 
”-i flag) ; 
if ( flag) 

{ 

int f easibility_count=0 ; 

for ( int f=0; f <vr_array . vrba . length; f++) 

{ 

if ( fb . feasibility (vr_ar ray . vrba [f ] ) ) 
f easibility__count++; 

} 

Variable []vrb__temp = new Variable [feasibility_count ] ; 
int f^count = 0; 

for (int f=0; f<vr_array . vrba . length; f+4-) 

{ 

if (fb. feasibility (vr__array. vrba [f] ) ) 

{ 

vrb temp [ f_count ] = vr__array . vrba [ f ] ; 
f_count+f ; 

} 

} 



vr_array__feasible = new Variable_array (vrb__temp) ; 

vrb = sit .select (vrb,vr_array_f easible,T,rm. nextDoubleO , 
rm_new. next Double ( ) ) ; 
double tempsol = ev. evaluation (vrb) ; 
if (solution<=tempsol) 

{ 

solution =tempsol; 
vrb_final = vrb; 

} 

} 

else 

break; 

} 

fw . write ( "tsolution+”\n'* ) ; 

T=T-^r; 

iterl++; 

) 

v3 . addElement ( vrb__f inal ) ; 

v4 , addElement (new Double (solution) ) ; 

} 


input_a = new Variable [v3. size ()] ; 
input_d = new double [v4 . size ()] ; 
for (int i=0; i<input_d. length; i++) 

input_d [i] =Double .parseDouble (v4 . elementAt (i) . toString { ) ) ; 
input_a [i] = (Variable) v3 . elementAt (i) ; 

} 

double tempi; 

Variable temp2 = new Variable () ; 

for (int i«0; i<input_d. length; i++} 

{ 

for (int j=i; j<input_d. length; j++) 

{ 

if (input_d[ j ] >input_d[i] ) 

{ 

tempi = input_d[j]; 
input_d[j]= input_d[i]; 
input_d[i] =templ; 
temp2 = input_a[j]; 
input_a[j]= input__a[i]; 
input__a [i] ==temp2; 

} 

} 

} 

Variable []input_al = new Variable [ f ine__tune_length] ; 
for ( int i = 0 ; i<f ine__tune_length; i + + ) 
input_al[i] = input_a[i]; 

input = new Variable_array ( input__al ) ; 

vrb__optimum = f t . fineTuning ( input ) ; 
solution_optimum = evl . evaluation (vrb_optimum) ; 
for (int i=0; i<vrb_optimum. s . length; i4-+) 

{ 

System, out .print In ( vrb_optimum. s [i] + ”\t”+vrb_optimum. e [i] + "\t”+vrb_op 
timum. s_n [i] ) ; 

} 



System, out .print In ( "sol_opt:="+solution_optimum) ; 

GregorianCalendar gcall = new GregorianCalendar(); 

System, out . print ( "'Time After: "); 

System, out .print (gcall . get (Calendar . HOUR) +" :”); 

System, out .print (gcall. get (Calendar .MINUTE) +” :”) ; 

System. out .println(gcall. get (Calendar . SECOND) ) / 

MyTimer stop = new MyTimerO; 

long minis = stop.difference(start); 

System. out .println("Program Took: ” + millis + ” milliseconds” ) ; 

} 

} 


2 . Select. java 

import j ava . ut il . * ; 
import java.io.*; 

class Select 

{ 

Select ( ) 

{ } 

public! si.nl \ r. Vnr i able so 1 ofU- ( Vn n' nhl e vr, Vnr i able a r ray vrr,doiablo 
t, double rml, double rm2) 
throws Exception 

Variable_array varray = new Variable_array ( ) ; 

Variable vrb = new VariableO; 

Variable vrb_temp = new VariableO; 
double T =t , probl, randoml=rml, random:2=rm2; 
int prob; 

vrb = vr; 
varray = vrr; 

Evaluation ev = new Evaluation () ; 

double interm = ( (randoml ) * (varray . vrba . length) ) ; 

prob= ( int ) interm; 

i f (prob==varray . vrba . length) 

vrb_temp=varray. vrba [prob- 1] ; 

else 

vrb___t emp=var ray . vrba [prob] ; 
double difference ; 

difference = (ev. evaluation {vrb___temp) ) - (ev. evaluation (vrb) ) ; 

if (difference >= 0) 

return (vrb__temp) ; 

else 

{ 

double exponent = dif f erence/t ; 

dcniblc! (u<fior\out._va lue ■ Math . pow (2 . 7 1B2H 1 , exponent ) ; 
if (random2<=exponent__value) 
return (vrb^temp) ; 
else 

return (vrb) ; 

} 

1 

} 



APPENDIX 5 

CODE FOR RESOLVING CLASHES IN THE ENTIRE VISIBILITY FILE 
WITH ONE OF THE ABOVE ALGORITHM 


ClashResolve . java 

import j ava . lang . * ; 
import java.io.*; 
import j ava . util . ; 

//THIS PROGRAM PICKS UP THE CLASH SITUATIOlsf FROM THE 
// FILE CLASDID1.TXT AND SENDS IT FOR CLASH RESOLUTION 
//AFTER RESOLVING CLASHES, IT REPLACES THE VISIBILITY START 
//AND VISIBILITY END TIMINGS WITH ALLOCATION START AND 
//ALLOCATION END TIMINGS. 

class TestArrayll 

{ 

String []s = new String[19]; 

TestArrayll ( ) 

{} 

TestArrayll {String[] str) 

{ 

f o r ( i n t i “ 0 ; i <19; i 4 4 ) 
s [i] ==str [i] ; 

} 

} 

class ClashResolve 

{ 

public static void main(String args[]) 
throws Exception 

MyTimer start = new MyTimerO; ^ 

GregorianCalendar gcal = new GregorianCalendar { ) / 

System, out .print ( '"Time Before: ”); 

System. out .print (gcal . get (Calendar . HOUR) 

System, out . print (gcal . get (Calendar . MINUTE) 4-" : " ) ; 

System . out . println (gcal . get (Calendar . SECOND) ) ; 

String svl , sv2 , name; 

StringTokenizer st k, stkl , st k2 ; 

TestArrayll oavl; 

int no_of__ttc_stations = 12; 

String [] station_name = new String [no__of_tt Gestations ] ; 

File fl “ new Fi le ( "clashIDl ” ) ; 

File f__ttc = new File { " ttcstation . txt ” ) ; 

FileReader frl = new FileReader { f 1) ; 

FileReader fr__ttc = new FileReader ( f^ttc) ; 

Buf f eredReader brl = new Buf f eredReader ( f r 1 ) ; 

Buf f eredReader br^ttc = new Buf f eredReader ( fr^ttc) ; 
FileWriter fwl = new FileWriter ( "clashID2” ) ; 

Vector V = new Vector (); 

while ( (svl=brl . readLine ( ) ) 1 =nuil) 

{ 



stkl = new StringTokenizer (svl ) ; 

String []savl = new String [stkl . countTokens ()] ; 

int tnv=0; 

while (stkl .hasMoreTokens {) ) 

{ 

savl [tnv]=stkl.nextToken() ; 
tnv++; 

} 

oavl = new TestArrayll(savl); 

V . addElement (oavl ) ; 

} 

V . trimToSize ( ) ; 
int count =0; 

while ( (name=br__ttc . readLine ( ) ) ! =null) 

{ 

station__name [count ] = name; 
count 'I- 4- ; 

} 

for (int k=0; k< (v. size ( ) ) ; k+4-) 

{ 

boolean isttc = false; 

for (int i=0; i<no_of__ttc__stations; i++) 

{ 

if ( (station_name [i] ) . equals ( ( (TestArrayll ) v. element At ( k) ) . s [5] . 
toString ( ) ) ) 

{ 

isttc = true; 

//System. out . print in { "if condition is satisfied" ) ; 
break; 

} 

} 

//System. out .println ( "I am here"+isttc+k) ; 
if (isttc) 

{ 

i f ( { ( ( Tost: Ar rayl 1 ) v . o 1 ornent At. { k ) ) . s [ 1 5 ] ) . oqua } s ( " ) ) 

f 

Vector vl new Vector (); 

vl . addElement ( v . elementAt ( k) ) ; 

for (int i=k+l ; i<k+100; i++) 

{ 

if (i<v. size ( ) ) 

{ 

1 f ( ( (Tes t Arrayl 1 ) V . elementAt ( i ) ) . s [ 5 ] . equals ( ( (TestArrayl 1 ) v . el 
eiuentAt ( k) ) , s [ 5 ] ) ) 

{ 

if ( ( (TestArrayll) v. elementAt (i) ) .s [0] .equals { ( (TestArrayll ) v. el 
ementAt {k)).s[0])) 

{ 

vl . addElement (v. elementAt (i) ) ; 

} 

else 

break; 

} 

if (vl . size ( ) >8) 
break; 

} 

else 

break; 



} 

if ( vl . size 0 >1 ) 

I 

iiitt ) vis_start « new int [ vl . size ( ) ] ; 
int[] vis_end = new int [ vl . size { ) ] ; 
int[] s__w = new int [vl . size ()] ; 
int[] s_n = new int [vl . size ()] ; 
long temp date ”0; 
for ( int i=0; i<vl . size ( ) ; i++) 

{ 

int date = Integer . parseint {( (TestArrayll) 

vl . elementAt ( i ) ) . s [3] . toString { ) ) ; 

i n t mon t h I n t: ege r . pa rse I n t ( ( ( Tes t Ar ray 1 ] ) 

vl . elementAt (i ) ) . s [2 ] . toString ( ) ) ; 

int year = Integer .parseint (( (TestArrayll) 

vl . elementAt ( i ) ) . s [ 1 ] . toString ( ) ) ; 

int vis_start_hr =Integer . parseint (( (Test Arrayl 1 ) 

v 1 . olomontAt: ( i ) ) . s ( 8 ) . toSt r.i ng ( ) ) ; 

int vis_start_mt ~ Integer . parseint (( (TestArrayl i ) 

vl. elementAt (i) ) .s [9] .toStringO ) ; 

int vis_start_sec = Integer .parseint (( (TestArrayll) 

vl. element At (i) ) .s [10] . toString () ) ; 

int vis_end_hr ===Integer .parseint ( ( (TestArrayll) 

vl . elementAt ( i ) ) . s f 11 ] . toString ( ) ) ; 

int vis_end_mt = Integer .parseint {{ (TestArrayll ) 

vl .elementAt (i) ) .s [12] . toString () ) ; 

jnt: vis^end^sec - Integer .parseint ( { (TestArrayll) 

vl . elementAt ( i ) ) . s [ 13] . toString ( ) ) ; 

Calendar cal = Calendar . getinstance () ; 

Calendar call = Calendar . getinstance () ; 

cal . set (year, month, date, vis_start_hr, vis_start__mt , vis_start sec 


call . set (year, month, date, vis__end_hr, vis__end_mt, vis_end_sec) ; 


if ( ( vis_start_hr==23 ) && ( vis_end_hr==0 ) ) 
call . add (Calendar . DATE, 1) ; 

Date before = cal . getTime ( ) ; 
long bef = before . getTime () ; 
i f ( i .^~^-: 0) 

I (mip date bef; 

Date after ^ cal 1 . getTime () ; 
long aft - after.getTime(); 

/ /System . out . print In ( "bef is: "+bef); 

/ /Sys t em . out . pr i nt In ( ’’ aft is:*’ + aft); 

V i s s t a r t [ i ] = ( i n t ) ( ( be f - 1 emp_da t e ) / 1 0 00 ) ; 
vis_end[i] - (int ) ( (af t-temp_date) /lOOO) ; 
s_n [i] = Integer .parseint ( ( (TestArrayll) 
vl .elementAt (i) ) .s [18] ) ; 
s w[i] ^ Integer . parseint (( (TestArrayll ) 
vl . elementAt ( i ) ) . s [ 17 ] ) ; 

} 

//pass this objet to clash resolution and get an variale object 
back.. //get them converted into 
//date and time format. 

//compared variable and input class objects... 

InputClass iclass = new InputClass (vis__start , vis_end, s__w, s_n) ; 
//one can the below line to tabu or simulated annealing. 



Greedyl greedyl = new GreedyK); 

Variable vrb = greedyl . greedyl (iclass) ; 

for (int i— 0^ i^^iclass . s n . length / 14 - 4 ') 

{ 

boolean match = false; 

for ( 1 nt i ^0 ; j < vrb . s . length ; j + 4 ) 

I 

i i ( ( i c 1 .1 .‘,5 s . s_n [ 1 J ) ( V i b . I J J ) ) 

{ 

match=true; 

break; 

} 

} 

if (match==f alse) 

{ 

Gtr ing (] toinp new String [19]; 
for(int j=0; j<15; j4-+) 

temp[j] = ( { (TestArrayll) v.elementAt {iclass.s_n[i] ) } .s[j] ) 
temp [15] = "Clash-Resolved/NS” ; 
for (int j = 16; j<19; j+4-) 

temp[ j] := ( ( (TestArrayll) v.elementAt ( iclass . s__n [i ] ) ) .s[j] ) 
TestArrayll tsl = new TestArrayll (temp) ; 

V . set Element At (tsl, iclass. s__n [i] ) ; 

} 

} 

int [ ] new_start_hour = new int [vrb, s . length] ; 
int [ ] new_start_minute = new int [vrb. s . length] ; 
int [] newest art_second = new int [vrb. s . length] ; 
int [ ] new_end_hour = new int [vrb. s . length] ; 
int [ ] new___end_minute = new int [vrb. s . length] ; 
int [ ] new_end_second = new int [vrb . s . length] ; 

f or ( int i“0 ; i<vrb . s . length; i4-4- ) 

Calendar cal_new = Calendar . getinstance () ; 

Date dt = new Date (( long) ( (vrb . s [i ]* 1000) 4'temp_date) ) ; 
ca i^now . so tT ime ( dt ) ; 

newest art_hour [i] = cal__new. get (Calendar .HOUR__OF_DAy) ; 
new_start_minute [i] = cal_new. get (Calendar . MINUTE) ; 
new_start_second [i] = cal_new. get (Calendar . SECOND) ; 

Calendar cal_newl = Calendar . getinstance (} ; 

Date dtl = new Date ( (long) ( (vrb. e [ i] *1000) +temp_date) ) ; 
cal_newl . set Time (dtl) ; 

new_end_hour [i] = cal_newl.get (Calendar . HOUR_OF_DAY) ; 
new_end_minute [i] = cal_newl . get (Calendar .MINUTE) ; 
new_end_second[i] = cal_newl. get (Calendar. SECOND) ; 

} 

for (int i=0; i<vrb . s . length; i4- + ) 

{ 

String []temp = new String[19]; 
for (int j =0 ; j <8 ; j 4-+ ) 
temp[j] = ((TestArrayll) 

V . elernentAt (vrb . s_n [ i ] ) ) . s [ j ] .-toSt ring ( ) ; 
temp [8] = ( ” "4-new_start_hour [ i ] ) ; 

temp[9] = (""4new___start__minute[i]); 
temp [10] = (”” + new_start__second[i]); 
temp [11] = ( ” ’'4-new_end_hour [ i ] ) ; 
temp [12] = ( "”+new__end_minute [i] ) ; 
temp [13] = (””4-new__end__second[i] ) ; 



temp [14] = ( (TestArra¥ll)^ p 

v.elementAt {vrb. s___n [if) f .toStringO; 

temp [15] - ”CR TTC"; 

for (int j==16; j<19; j4 1) 

temp[j] *= ( (TestArrayll) 

V.elementAt (vrb.s__n[i] ) ) .s[j] .toStringO ; 

TestArrayll tsl = new TestArrayll (temp) ; 

V , set FM ernen t At (tsl, vrb . s n [ i ] ) ; 

) 

}//if v.size end 
}//if free end 
) //if isttc end 
}//for loop end 

for(int i=0 ; i<v . size ( ) ; i+4) 

for(int j=0; j<19; j++) 

fwl, write ( ( (TestArrayll) v, elementAt (i) ) .s [ j ] +”\t”) ; 
fwl . write ( ) ; 

} 


GregorianCalendar gcall = new GregorianCaiendarO; 

System. out .print ("Time After: ”); 

System. out . print (gcall . get (Calendar . HOUR) +” : ; 

System, out,. print (gcall . get (Calendar .MINUTE) +” : ” ) ; 
System. out .println (gcall .get (Calendar. SECOND) ) ; 

MyTimer stop = new MyTimerO; 

long minis = stop.difference(start); 

System. out .println ( "Program Took: " + millis + ” milliseconds”); 

fwl . close ( ) ; 
f rl . close ( ) ; 


} 

} 



